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 649 records · Page 36

Measurements of Aerodynamic Heat Transfer and Boundary-Layer Transition on a 15 deg. Cone in Free Flight at Supersonic Mach Numbers up to 5.2

Measurements of aerodynamic heat transfer have been made at several stations on the 15 deg total-angle conical nose of a rocket-propelled model in free flight at Mach numbers up to 5.2. Data are presented for a range of local Mach number just outside the boundary layer from 1.40 to 4.65 and a range of local Reynolds number from 3.8 x 10(exp 6) to 46.5 x 10(exp 6), based on length from the nose tip to a measurement station. Laminar, transitional, and turbulent heat-transfer coefficients were measured. The laminar data were in agreement with laminar theory for cones, and the turbulent data agreed well with turbulent theory for cones using Reynolds number based on length from the nose tip. At a nearly constant ratio of wall to local static temperature of 1.2 the Reynolds number of transition increased from 14 x 10(exp 6) to 30 x 10(exp 6) as Mach number increased from 1.4 to 2.9 and then decreased to 17 x 10(exp 6) as Mach number increased to 3.7. At Mach numbers near 3.5, transition Reynolds numbers appeared to be independent of skin temperature at skin temperatures very cold with respect to adiabatic wall temperature. The transition Reynolds number was 17.7 x 10(exp 6) at a condition of Mach number and ratio of wall to local static temperature near that for which three-dimensional disturbance theory has been evaluated and has predicted laminar boundary-layer stability to very high Reynolds numbers (approximately 10(exp 12)).

Rumsey, Charles B.↗

New Hubble Space Telescope Observations of High-Velocity Ly(alpha) and H(alpha) in SNR 1987A

We describe and model high-velocity (approximately 15,000 km S(exp -1)) Ly Alpha and H Alpha emission from the supernova remnant SNR 1987A seen in 1997 September and October with the Space Telescope Imaging Spectrograph. Part of this emission comes from a reverse shock located at approximately 75% of the radius of the inner boundary ofthe innercircumstellar ring and confined within + or - 30 degrees of the equatorial plane. Departure from axisymmetry in the Ly Alpha and H Alpha emission correlates with that seen in nonthermal radio emission and reveals an asymmetry in the circumstellar gas distribution. We also see diffuse high-velocity Ly-Alpha emission from supernova debris inside the reverse shock that may be due to excitation by nonthermal particles accelerated by the shock.

Michael, Eli↗

Regulatory and Technical Issues Concerning the Detection and Treatment of NDMA-Contaminated Groundwater at NASA WSTF

The National Aeronautics and Space Administration (NASA) White Sands Test Facility (WSTF) was established in 1963 primarily to provide rocket engine testing services for several NASA programs. The groundwater underlying the site has been contaminated as a result of historical operations. Groundwater contaminants include several volatile organic compounds (VOCs) and two semi-volatile compounds: N-nitrosodimethylamine (NDMA) and N-nitrodimethylamine (DMN). This paper discusses some of the technical, analytical, regulatory, and health risk issues associated with the contaminant plume. The plume has moved approximately 2.5 miles downgradient of the facility industrial boundary, with evidence of continued migration. As a result, NASA has proposed a pump and treat system using air strippers and ultraviolet (UV) oxidation to stabilize future movement of the contaminant plume. The system has been designed to treat 1,076 gallons (4,073 liters) per minute, with provisions for future expansion. The UV oxidation process was selected to treat NDMA-contaminated groundwater based on successes at other NDMA-contaminated sites. Bench- and pilot-scale testing of WSTF groundwater confirmed the ability of UV oxidation to destroy NDMA and generated sufficient data to design the proposed full-scale treatment system. NDMA is acutely toxic and is a probable human carcinogen. EPA-recommended health risk criteria for the residential consumption of NDMA/DMN-contaminated groundwater was used to determine that a 1.0 x 10(exp -6) excess cancer risk corresponds to 1.7 parts per trillion (ppt). EPA analytical methods are unable to detect NDMA and DMN in the low ppt range. EPA's current Appendix IX analytical method used to screen for NDMA, Method 8270, can detect NDMA only at levels that are orders of magnitude greater than the recommended health risk level. Additionally, EPA Method 607, the most sensitive EPA approved method, has a detection limit of 150 ppt. This corresponds to an excess cancer risk of 9.0 x 10(exp -5), which exceeds the State of New Mexico's water quality standard of a cancer risk less than 1 x 10(exp -5). The treatment system has been engineered to treat contaminated groundwater to levels significantly below the New Mexico standard. However, the inability of EPA-approved analytical methods to detect NDMA and DMN at low ppt levels, and to provide verification of compliance with the 1 x 10(exp -5) cancer risk, introduces a notable risk to the long-term operation of the system. WSTF has been working with Southwest Research Institute to develop a non-EPA analytical method that can achieve a reporting limit of 1 ppt, which corresponds to an excess cancer risk of 7.6 x 10(exp -7). WSTF is currently developing a proposal to obtain approval from the New Mexico Environment Department (NMED) of this non-EPA method.

Wiebe, D. T.↗

A Parallel Multiblock Mesh Movement Scheme For Complex Aeroelastic Applications

A scheme has been developed for the movement of multiblock, structured grids due to surface deformation arising from aeroelastics, control surface movement, or design optimization. Elements of the method include a blending of a surface spline approximation and nearest surface point movement for block boundaries. Transfinite interpolation is employed for volume grid deformation. The scheme is demonstrated on a range of simple and complex aeroelastic aircraft applications using Navier-Stokes computational fluid dynamics and modal structural analyses on parallel processors. Results are robust and accurate, requiring only minimal user input specification.

Potsdam, Mark A.↗

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗

The Gale Crater Mound in a Regional Geologic Setting: Mapping and Probing Surrounding Outcrops for Areas Akin to the Central Mound at Gale

There are several hypotheses on the origin of Gale Crater s central mound. These include ground water upwelling [1], aeolian, ice, volcanic [1-3], lacustrine [1-3], hydrothermal [1-3], and polar deposits [2]. The Mars Science Laboratory rover, Curiosity, landed in Gale Crater on August 6, 2012. It is currently analyzing samples along its traverse towards a channel and layered deposits that will provide insight into the sedimentary history of the crater [4]. Located at 5S, 138E, Gale is a 155km diameter, Late Noachian/Early Hesperian crater. It is situated along the southern highlands/northern lowlands dichotomy boundary and contains a central mound that rises approximately 5km from the crater floor [1]. The highest parts of Mt. Sharp are higher than the northern rim, but are roughly the same height as the southern rim. Mt. Sharp is divided into an upper mound and a lower mound, which are separated by an erosional unconformity [2]. The lower mound s sequences span the Late Noachian/Early Hesperian Epoch [1], while the upper mound s age is poorly constrained. The lower mound s sequences feature parallel beds of varying thickness, albedo, texture, and dip angle that are eroded into channels and yardangs [2]. The upper mound has finer layers at higher angles [1] with yardangs, serrated erosional patterns, and lobate features [3]. The lower mound also exhibits an upward progression of phyllosilicate to sulfate rich sediments, contrasting the upper mound s lack of hydrated minerals [4].

Korn, Lisa↗

Changes in Flat Plate Wake Characteristics Obtained With Decreasing Plate Thickness

The near and very near wake of a flat plate with a circular trailing edge is investigated with data from direct numerical simulations. Computations were performed for four different Reynolds numbers based on plate thickness (D) and at constant plate length. The value of θ/D varies by a factor of approximately 20 in the computations (θ being the boundary layer momentum thickness at the trailing edge). The separating boundary layers are turbulent in all the cases. One objective of the study is to understand the changes in wake characteristics as the plate thickness is reduced (increasing θ/D). Vortex shedding is vigorous in the low θ/D cases with a substantial decrease in shedding intensity in the largest θ/D case (for all practical purposes shedding becomes almost intermittent). Other characteristics that are significantly altered with increasing θ/D are the roll-up of the detached shear layers and the magnitude of fluctuations in shedding period. These effects are explored in depth. The effects of changing θ/D on the distributions of the time-averaged, near-wake velocity statistics are discussed.

Wakes↗

Comparison of Staggered Grid Finite Difference Schemes for Ultrasound Simulation in Curving Composites

The optimization of ultrasonic nondestructive evaluation (NDE) simulation tools for composites has the potential to reduce both individual part inspection time and overall certification time for composite parts and structures. Inspection guidance based on simulation provides increased confidence in the veracity of inspection results in addition to time reductions. This paper outlines ongoing work targeted to advance this objective through the use of finite difference (FD) simulation techniques formulated for composite structures with realistic geometries. Two staggered grid explicit FD schemes which show promise for this purpose are assessed: the Lebedev FD scheme and the rotated staggered grid (RSG) FD scheme. Algorithmic points which provide challenges for complex geometries are addressed, in particular handling of traction free surfaces and bi-material interfaces present at lamina boundaries. Code execution time estimates are performed as well to guide feasible domain sizes relative to algorithm choice and available hardware. Three test cases are simulated: a delaminated plate, a cylinder, and a triclinic lamina. These tests demonstrate that the Lebedev FD scheme needs additional work to handle inter-laminar interfaces and traction free boundaries in the presence of stair-stepping approximations. In contrast, the simple structure of the RSG unit cell makes it more straightforward to construct a 3D simulation technique for curved composite laminates.

Frankforter, Erik↗

Aerosol Size Distribution Properties Associated With Cold-Air Outbreaks in the Norwegian Arctic

The aerosol particles serving as cloud condensation and ice nuclei contribute to key cloud processes associated with cold-air outbreak (CAO) events but are poorly constrained in climate models due to sparse observations. Here we retrieve aerosol number size distribution modes from measurements at Andenes, Norway, during the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) and at Zeppelin Observatory, approximately 1000 km upwind from Andenes at Svalbard. During CAO events at Andenes, the sea-spray-mode number concentration is correlated with strong over-ocean winds with a mean of 8±4 cm−3 that is 71 % higher than during non-CAO conditions. Additionally, during CAO events at Andenes, the mean Hoppel minimum diameter is 6 nm smaller than during non-CAO conditions, though the estimated supersaturation is lower, and the mean number concentration of particles that likely activated in-cloud is 109±61 cm−3 with no statistically significant difference from the non-CAO mean of 99±66 cm−3. For CAO trajectories between Zeppelin Observatory and Andenes, the upwind-to-downwind change in number concentration is the largest for the accumulation mode with a mean decrease of 93±95 cm−3, likely attributable primarily to precipitation scavenging. These characteristic properties of aerosol number size distributions during CAO events provide guidance for evaluating CAO aerosol–cloud interaction processes in models.

accumulation mode particle scavenging↗

Design Considerations of an Ascent Abort Monitor Algorithm for Use During Service Module Aborts

In support of human rating the Artemis missions, NASA's Orion program requires continuous abort coverage from liftoff through mission destination. During a portion of the ascent trajectory, the currently achievable abort mode is determined by an Orion algorithm using the onboard navigated vehicle state. This ascent abort monitor determines achievability for Orion's Mode 2 abort, Untargeted Abort Splashdown (UAS), by propagating the current vehicle state through ascent abort events to determine sufficient timing to perform the abort and to a ballistic touchdown point to approximate landing location relative to prescribed keep out boundaries. The algorithm was updated for Artemis 2 to allow the capability to limit loads for the majority of ascent. Performance of the algorithm has been demonstrated and verified through dispersed trajectory analysis with emulated flight software, unit testing, and hardware in the loop testing.

Esteban Guzman↗

Transient-forced convection film boiling on an isothermal flat plate.

Development of a new approach for the solution of transient-forced convection film boiling on an isothermal flat plate using the boundary layer model. The similarity variables are used to convert the governing partial differential equations to ordinary ones. The results of numerical solutions of these ordinary equations indicate that the transient process can be classified as one-dimensional conduction, intermediate, and the steady-state regions. The time required for the one-dimensional conduction and the time necessary to attain a steady-state condition are obtained. The use of local similarity approximations for the intermediate regime facilitates prediction of complete boundary layer growth. Using the ratio of time at any instant to the steady-state time as abscissa, the curves representing the boundary layer growth can be merged into a single mean curve within 5%. Further, the analysis shows that the average rate of heat transfer during transient is 50 to 100% higher than those at steady state. The average rate of vapor convected away is 10 to 15% lower than at steady state while the average rate of accumulation to form the vapor layer is 1 to 14 times larger.

Nagendra, H. R.↗

A Model for the Sources of the Slow Solar Wind

Models for the origin of the slow solar wind must account for two seemingly contradictory observations: The slow wind has the composition of the closed-field corona, implying that it originates from the continuous opening and closing of flux at the boundary between open and closed field. On the other hand, the slow wind has large angular width, up to approximately 60 degrees, suggesting that its source extends far from the open-closed boundary. We propose a model that can explain both observations. The key idea is that the source of the slow wind at the Sun is a network of narrow (possibly singular) open-field corridors that map to a web of separatrices and quasi-separatrix layers in the heliosphere. We compute analytically the topology of an open-field corridor and show that it produces a quasi-separatrix layer in the heliosphere that extends to angles far front the heliospheric current sheet. We then use an MHD code and MIDI/SOHO observations of the photospheric magnetic field to calculate numerically, with high spatial resolution, the quasi-steady solar wind and magnetic field for a time period preceding the August 1, 2008 total solar eclipse. Our numerical results imply that, at least for this time period, a web of separatrices (which we term an S-web) forms with sufficient density and extent in the heliosphere to account for the observed properties of the slow wind. We discuss the implications of our S-web model for the structure and dynamics of the corona and heliosphere, and propose further tests of the model.

Antiochos, Spiro K.↗

Melting of iron by significant structure theory

From Eyring's method of significant structures, a partition function is derived for liquid iron. The solid at high temperature is described in the Einstein approximation. Magnetic and electronic contributions to the thermodynamic properties of both the liquid and solid phases are considered. The model is compatible with properties (thermal expansion, compressibility, heat capacity, entropy of melting, and volume change on melting) at one atm. The melting temperature at high pressure is found by satisfying the requirement that the Gibbs free energies of the liquid and solid phases are equal at the melting temperature. Under conditions at the earth's core-mantle boundary, the melting temperature of iron is greater than approximately 5000 K, and under inner-outer core conditions the melting temperature is greater than approximately 7000 K. These estimates are consistent with the Lindemann melting law, but not with the Kraut-Kennedy melting law.

Leppaluoto, D. A.↗

Boundary Layer Effect on the Surface Pressure of an Infinite Cone in Supersonic Flow

The theory of Taylor and Maccoll (Ref,1) gives the surface pressure on an infinite cone in supersonic flow as a function of the cone vertex angle and the free stream Mach number and static pressure for a gas of vanishing viscosity. When a slender conical probe is used together with an impact pressure probe to determine the static pressure and Mach number in a low density gas stream, it is desirable to have some theoretical estimate of the effect of viscous boundary layer on the probe readings. Theoretical and experimental results with respect to impact probes have been presented in Refs. 5 and 6. A simple approximation for a conical probe based on linearized supersonic flow and compressible boundary layer theory is presented here.

Lin, T. C.↗

Forward marching procedure for separated boundary-layer flows

A forward-marching procedure for separated boundary-layer flows which permits the rapid and accurate solution of flows of limited extent is presented. The streamwise convection of vorticity in the reversed flow region is neglected, and this approximation is incorporated into a previously developed (Carter, 1974) inverse boundary-layer procedure. The equations are solved by the Crank-Nicolson finite-difference scheme in which column iteration is carried out at each streamwise station. Instabilities encountered in the column iterations are removed by introducing timelike terms in the finite-difference equations. This provides both unconditional diagonal dominance and a column iterative scheme, found to be stable using the von Neumann stability analysis.

Carter, J. E.↗

Inhomogeneous conditions at open boundaries for wave propagation problems

Absorbing boundary conditions contain differential operators even for first-order systems. There is a fundamental difficulty with this, since the conditions are applied on the ingoing variables, and the approximations necessarily become weakly unstable. This difficulty is more pronounced for inhomogeneous boundary conditions which occur if there is a source outside the computational domain D, or if the initial data are nonzero outside D. In this paper this is further investigated and it is shown that reasonable estimates can still be obtained if the solution is smooth. However, it is demonstrated that the approximations are less robust. A cure for this is proposed and an implementation of high-order conditions for first-order systems is described.

Gustafsson, Bertil↗

A comparative study of an ABC and an artificial absorber for truncating finite element meshes

The type of mesh termination used in the context of finite element formulations plays a major role on the efficiency and accuracy of the field solution. The performance of an absorbing boundary condition (ABC) and an artificial absorber (a new concept) for terminating the finite element mesh was evaluated. This analysis is done in connection with the problem of scattering by a finite slot array in a thick ground plane. The two approximate mesh truncation schemes are compared with the exact finite element-boundary integral (FEM-BI) method in terms of accuracy and efficiency. It is demonstrated that both approximate truncation schemes yield reasonably accurate results even when the mesh is extended only 0.3 wavelengths away from the array aperture. However, the artificial absorber termination method leads to a substantially more efficient solution. Moreover, it is shown that the FEM-BI method remains quite competitive with the FEM-artificial absorber method when the FFT is used for computing the matrix-vector products in the iterative solution algorithm. These conclusions are indeed surprising and of major importance in electromagnetic simulations based on the finite element method.

Oezdemir, T.↗

Validation of Blockage Interference Corrections in the National Transonic Facility

A validation test has recently been constructed for wall interference methods as applied to the National Transonic Facility (NTF). The goal of this study was to begin to address the uncertainty of wall-induced-blockage interference corrections, which will make it possible to address the overall quality of data generated by the facility. The validation test itself is not specific to any particular modeling. For this present effort, the Transonic Wall Interference Correction System (TWICS) as implemented at the NTF is the mathematical model being tested. TWICS uses linear, potential boundary conditions that must first be calibrated. These boundary conditions include three different classical, linear. homogeneous forms that have been historically used to approximate the physical behavior of longitudinally slotted test section walls. Results of the application of the calibrated wall boundary conditions are discussed in the context of the validation test.

Walker, Eric L.↗