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 955 records · Page 53

Pseudospectral collocation methods for fourth order differential equations

Collocation schemes are presented for solving linear fourth order differential equations in one and two dimensions. The variational formulation of the model fourth order problem is discretized by approximating the integrals by a Gaussian quadrature rule generalized to include the values of the derivative of the integrand at the boundary points. Collocation schemes are derived which are equivalent to this discrete variational problem. An efficient preconditioner based on a low-order finite difference approximation to the same differential operator is presented. The corresponding multidomain problem is also considered and interface conditions are derived. Pseudospectral approximations which are C1 continuous at the interfaces are used in each subdomain to approximate the solution. The approximations are also shown to be C3 continuous at the interfaces asymptotically. A complete analysis of the collocation scheme for the multidomain problem is provided. The extension of the method to the biharmonic equation in two dimensions is discussed and results are presented for a problem defined in a nonrectangular domain.

Malek, Alaeddin↗

Effects of Adverse Pressure Gradient on the Incompressible Reattaching Flow Over a Rearward-Facing Step

The turbulent, incompressible reattaching flow over a rearward-facing step has been studied by many researchers over the years. One of the principal quantities determined in these experiments has been the distance from the step to the point (or region) where the separated shear layer reattaches to the surface (x(r)). The values for x(r)/h, where h is the step height, have covered a wider range than can reasonably be attributed to experimental technique or inaccuracy. Often the reason for a largely different value of x(r)/h can be attributed to an incompletely developed turbulent layer, or a transitional or laminar boundary layer. However, for the majority of experiments where the boundary layer is believed to be fully developed and turbulent, x(r)/h still varies several step heights; generally, 5 1/2 approximately < x(r)/h approximately < 7 1/2. This observed variation has usually been attributed to such variables as l/h (step length to height, h/delta (step height to initial boundary-layer thickness), R(e)(theta)), or the experimental technique for determining reattachment location. However, there are so many different combinations of variables in the previous experiments that it was not possible to sort out the effects of particular conditions on the location of reattachment. In the present experiment velocity profiles have been measured in and around the region of separated flow. Results show a large influence of adverse pressure gradient on the reattaching flow over a rearward-facing step that has not been reported previously. Further, the many previous experiments for fully developed, turbulent flow in parallel-walled channels have shown a range of reattachment location that has not been explained by differences in initial flow conditions. Although these initial flow conditions might contribute to the observed variation of reattachment location, it appears that the pressure gradient effect can explain most of that variation.

Kuehn, Donald M.↗

A method for calculating turbulent boundary layers and losses in the flow channels of turbomachines

An interactive inviscid core flow-boundary layer method is presented for the calculation of turbomachine channel flows. For this method, a one-dimensional inviscid core flow is assumed. The end-wall and blade surface boundary layers are calculated using an integral entrainment method. The boundary layers are assumed to be collateral and thus are two-dimensional. The boundary layer equations are written in a streamline coordinate system. The streamwise velocity profiles are approximated by power law profiles. Compressibility is accounted for in the streamwise direction but not in the normal direction. Equations are derived for the special cases of conical and two-dimensional rectangular diffusers. For these cases, the assumptions of a one-dimensional core flow and collateral boundary layers are valid. Results using the method are compared with experiment and good quantitative agreement is obtained.

Schumann, Lawrence F.↗

A method for calculating turbulent boundary layers and losses in the flow channels of turbomachines

An interactive inviscid core flow-boundary layer method is presented for the calculation of turbomachine channel flows. For this method, a one-dimensional inviscid core flow is assumed. The end-wall and blade surface boundary layers are calculated using an integral entrainment method. The boundary layers are assumed to be collateral and thus are two-dimensional. The boundary layer equations are written in a streamline coordinate system. The streamwise velocity profiles are approximated by power law profiles. Compressibility is accounted for in the streamwise direction but not in the normal direction. Equations are derived for the special cases of conical and two-dimensional rectangular diffusers. For these cases, the assumptions of a one-dimensional core flow and collateral boundary layers are valid. Results using the method are compared with experiment and good quantitative agreement is obtained.

Schumann, Lawrence F.↗

Boundary Layer Transition in the Leading Edge Region of a Swept Cylinder in High Speed Flow

Experiments were conducted on a 76 degree swept cylinder to establish the behavior of the attachment line transition process in a low-disturbance level, Mach number 1.6 flow. For a near adiabatic wall condition, the attachment-line boundary layer remained laminar up to the highest attainable Reynolds number. The attachment-line boundary layer transition under the influence of trip wires depended on wind tunnel disturbance level, and a transition onset condition for this flow is established. Internal heating raised the surface temperature of the attachment line to induce boundary layer instabilities. This was demonstrated experimentally for the first time and the frequencies of the most amplified disturbances were determined over a range of temperature settings. Results were in excellent agreement to those predicted by a linear stability code, and provide the first experimental verification of theory. Transition onset along the heated attachment line at an R-bar of 800 under quiet tunnel conditions was found to correlate with an N factor of 13.2. Increased tunnel disturbance levels caused the transition onset to occur at lower cylinder surface temperatures and was found to correlate with an approximate N factor of 1 1.9, so demonstrating that the attachment-line boundary layer is receptive to increases in the tunnel disturbance level.

Coleman, Colin P.↗

Coronal magnetic structure at a solar sector boundary

A persistent large-scale coronal magnetic structure associated with a sector boundary is shown to consist probably of a magnetic arcade loop structure extending from one solar polar region to the other in approximately the north-south direction. This structure was inferred from computed coronal magnetic field maps for days on which a stable magnetic sector boundary was near the central meridian, based on an interplanetary sector boundary observed to recur during much of 1968 and 1969.

Wilcox, J. M.↗

Combining Comparison Functions and Finite Element Approximations in CFD

In a variety of potential flow applications, the modal element method has been shown to significantly reduce the numerical grid, employ a more precise grid termination boundary condition, and give theoretical insight to the flow physics. The method employs eigenfunctions to replace the numerical grid over significant portions of the flow field. Generally, a numerical grid is employed around obstacles with complex geometry while eigenfunctions are applied to regions in the flow field where the boundary conditions can easily be satisfied. To handle a wider class of computational fluid dynamics (CFD) problems, the present paper extends the modal element to include function approximations which do not satisfy the governing differential equation. To accomplish this task, a double modal series approximation and weighted residual constraints are developed to force the comparison functions to satisfy the governing differential equation and to interface properly with the finite element solution. As an example, the method is applied to the problem of potential flow in a channel with two-dimensional cylindrical like obstacles. The calculated flow fields are in excellent agreement with exact analytical solutions.

Baumeister, Kenneth J.↗

A High-Lift Building Block Flow: Turbulent Boundary Layer Relaminarization

A working wind tunnel test facility has been constructed at the University of Notre Dame's Hessert Center. The relaminarization test facility has been constructed in the 1.5m x 1.5m (5ft x 5 ft) atmospheric wind tunnel and generates a Re(theta)=4694 turbulent boundary layer in nominally zero-pressure gradient before it is exposed to the Case #1 pressure gradient (K approximately equal to 4.2 x 10(exp -6), which is believed to be sufficient to achieve relaminarization. Future work to be conducted will include measuring the response of the turbulent boundary layer to the favorable pressure gradients created in the test facility and documenting this response in order to understand the underlying flow physics responsible for relaminarization. It is the goal of this research to have a better understanding of accelerated turbulent boundary layers which will aid in the development of future flow diagnostic utilities to be implemented in applied aerodynamic research.

Bourassa, Corey↗

Substorm-injected protons and electrons and the injection boundary model

Analysis of observations of substorm-associated enhancements of proton and electron fluxes: (1) makes a strong case for the existence of a boundary limiting the regional particle injection associated with substorms, (2) supports the hypothesis that the injection process is almost instantaneous (less than approximately 5 min), and (3) indicates that the injection takes place within a large region extending at least several earth radii tailward of the injection boundary. The injection boundary model is superior to others in that it simultaneously explains: (1) the drift and energy dispersion of substorm-injected protons between 1 and 30 keV, (2) the relative behavior of protons with 81 deg and 27 deg pitch angles, (3) the absence of observed electrons below 30 keV, and (4) the time dispersion of impulsively injected electrons seen outside the plasmapause.

Konradi, A.↗

Approximation methods for control of acoustic/structure models with piezoceramic actuators

The active control of acoustic pressure in a 2-D cavity with a flexible boundary (a beam) is considered. Specifically, this control is implemented via piezoceramic patches on the beam which produces pure bending moments. The incorporation of the feedback control in this manner leads to a system with an unbounded input term. Approximation methods in this manner leads to a system with an unbounded input term. Approximation methods in the context of linear quadratic regulator (LQR) state space control formulation are discussed and numerical results demonstrating the effectiveness of this approach in computing feedback controls for noise reduction are presented.

Banks, H. T.↗

Analytical formulas for conditions on blunt wedges in hypersonic flow

An analytical approximation is obtained to the solution of the equation describing the combined effects of bluntness and boundary-layer displacement on the hypersonic flow over a wedge for the case where the constant parameter in the governing equation, which is proportional to the angle of attack, is greater than zero. It is shown that the approximation has good physical accuracy over the entire range of conditions to which the theory is applicable for positive angle of attack.

Bade, W. L.↗

Eigenvalue Sensitivity Computations for Linear Stability Theory

To realize the drag reduction benefit of boundary-layer transition control strategies, it is crucial to integrate transition prediction into the vehicle design through an optimization process. The integration of transition prediction based on linear stability analysis into adjoint d design optimization requires coupling an adjoint enabled computational fluid dynamics (CFD) solver with an adjoint enabled linear stability code. In particular, the boundary-layer transition location is often predicted using the N-factor method based on linear stability theory (LST). Thus, sensitivity of the linear-stability eigenvalues constitute an essential building block for optimizing the laminar flow performance. The present paper describes an implementation of LST eigenvalue sensitivity analysis that can be easily coupled with a CFD solver. Specifically, we describe a discrete adjoint formulation for the transition location prediction based on the N-factor method. The verification of this formulation is carried out by comparing the adjoint-based sensitivity of the local growth rate of a given instability mode with respect to the disturbance frequency, and the adjoint-based sensitivity of the transition location with respect to spanwise wavenumber with those sensitivities computed using a finite-difference approximation. Finally, the adjoint LST formulation is applied to flat-plate boundary-layer flows at transonic, supersonic, and hypersonic conditions, to determine the behavior and sensitivities of the transition location with respect to a range of disturbance spanwise wavenumbers.

Boundary Layer Transition↗

A two-point difference scheme for computing steady-state solutions to the conservative one-dimensional Euler equations

An implicit finite-difference method is presented for obtaining steady-state solutions to the time-dependent, conservative Euler equations for flows containing shocks. The method uses a two-point central-difference scheme for the flux derivatives with dissipation added at supersonic points via the retarded density concept. Application of the method to 1-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions show the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. For 1-dimensional Euler calculations, it is shown that the scheme offers two advantages over the more widely-used three-point schemes. The first is in regard to application of boundary conditions, and the second relates to the fact that the two-point algorithm is well-conditioned for large time steps.

Wornom, S. F.↗

Fomalhaut's Debris Disk and Planet: Constraining the Mass of Formalhaut B from Disk Morphology

Following the optical imaging of exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhaut's debris disk is gravitationally shaped by a single interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass M(sub pl) < 3M(sub J), an orbital semimajor axis a(sub pl) > 101.5AU, and an orbital eccentricity e(sub pl) = 0.11 - 0.13. These conclusions are independent of Fom b's photometry. To not disrupt the disk, a greater mass for Fom b demands a smaller orbit farther removed from the disk; thus, future astrometric measurement of Fom b's orbit, combined with our model of planet-disk interaction, can be used to determine the mass more precisely. The inner edge of the debris disk at a approximately equals 133AU lies at the periphery of Fom b's chaotic zone, and the mean disk eccentricity of e approximately equals 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planet's chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of approximately 100 Myr, and model them separately from their dust grain progeny; the latter's orbits are strongly affected by radiation pressure and their lifetimes are limited to approximately 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fomalhaut b's nominal space velocity does not bear this out, but the astrometric uncertainties are difficult to quantify. Even if the apsidal misalignment proves real, our calculated upper mass limit of 3 M(sub J) still holds. Parent bodies are evacuated from mean-motion resonances with Fom b; these empty resonances are akin to the Kirkwood gaps opened by Jupiter. The belt contains at least 3M(sub Earth) of solids that are grinding down to dust, their velocity dispersions stirred so strongly by Fom b that collisions are destructive. Such a large mass in solids is consistent with Fom b having formed in situ.

Chiang, E.↗

Stability of the laminar boundary layer in a streamwise corner

The stability of viscous, incompressible flow along a streamwise corner, often called the corner boundary layer problem is examined. The semi-infinite boundary value problem satisfied by small amplitude disturbances in the "bending boundary layer' region is obtained. The mean secondary flow induced by the corner exhibits a flow reversal in this region. Uniformly valid "first approximations' to solutions of the governing differential equations are derived. Uniformity at infinity is achieved by a suitable choice of the large parameter and use of an approximate Langer variable. Approximations to solutions of balanced type have a phase shift across the critical layer which is associated with instabilities in the case of two dimensional boundary layer profiles.

Lakin, W. D.↗

Stability of the laminar boundary layer in a streamwise corner

The stability of viscous, incompressible flow along a streamwise corner, often called the corner boundary layer problem is examined. The semi-infinite boundary value problem satisfied by small amplitude disturbances in the 'bending boundary layer' region is obtained. The mean secondary flow induced by the corner exhibits a flow reversal in this region. Uniformly valid 'first approximations' to solutions of the governing differential equations are derived. Uniformity at infinity is achieved by a suitable choice of the large parameter and use of an approximate Langer variable. Approximations to solutions of balanced type have a phase shift across the critical layer which is associated with instabilities in the case of two dimensional boundary layer profiles. Previously announced in STAR as N84-17532

Lakin, W. D.↗

Formation mechanisms of Sn-rich δ phase and its role in strengthening Cu-10Sn manufactured by laser powder bed fusion

Cu-Sn alloys produced via laser powder bed fusion (L-PBF) additive manufacturing (AM) have gained significant attention because they combine the advantages of AM relevant to intricate component design with outstanding combinations of strength, ductility, and resistance to wear and corrosion. However, a detailed understanding of the microstructure that contributes to the enhancement of the mechanical properties of L-PBF Cu-10Sn alloys remains unclear. In particular, there is a lack of understanding of the formation mechanisms of the Sn-rich δ phase commonly observed in Cu-10Sn. This study reveals two distinct variants of the δ phase possessing unique morphological characteristics. These characteristics are attributed to the local solidification conditions inherent to the melt pool boundaries versus those at the interiors of melt pools. A phase transformation pathway that elucidates the origin of the morphological variants of the δ phase from the Sn-rich metastable phases during the cyclic heating of the AM process is proposed. We report superior mechanical properties in L-PBF Cu-10Sn compared to those of conventionally manufactured counterparts due to the synergistic contributions from grain boundaries, dislocations, and the δ phase. Notably, the δ phase alone contributes approximately 22 % to the overall strength observed in the L-PBF Cu-10Sn alloy. The discovery of two types of distinct Sn-rich δ phase offers key insights into precise microstructural control in AM Cu-Sn alloys to enhance mechanical properties, providing practical strategies for improving material performance for diverse applications in automotive, aerospace, and machinery industries.

36 MATERIALS SCIENCE↗