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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Calculation of convective heating on proposed aeroassist flight experiment vehicle

Convective heating calculations of a NASA aeroassist flight experiment (AFE) configuration consisting of a blunted, elliptic cone raked off at the base, and fitted with a skirt-type afterbody with a generous corner radius, are presented. The HALIS inviscid flowfield code is used to supply edge properties to an approximate three-dimensional boundary-layer method to calculate the surface heating. Calculations at Mach number 9.86 and approximate Reynolds number of 500,000, for angles of attack from -10 to 10 degrees, are in general good agreement with Langley 31-inch Mach 10 Tunnel experimental data. The results validate the application of the present overall AFE configuration heating pattern heating contours to aeroassist orbital transfer vehicle heating distribution determinations.

Hamilton, H. H., II↗
NASA NTRS1991

The asymptotic spectra of banded Toeplitz and quasi-Toeplitz matrices

Toeplitz matrices occur in many mathematical, as well as, scientific and engineering investigations. This paper considers the spectra of banded Toeplitz and quasi-Toeplitz matrices with emphasis on non-normal matrices of arbitrarily large order and relatively small bandwidth. These are the type of matrices that appear in the investigation of stability and convergence of difference approximations to partial differential equations. Quasi-Toeplitz matrices are the result of non-Dirichlet boundary conditions for the difference approximations. The eigenvalue problem for a banded Toeplitz or quasi-Toeplitz matrix of large order is, in general, analytically intractable and (for non-normal matrices) numerically unreliable. An asymptotic (matrix order approaches infinity) approach partitions the eigenvalue analysis of a quasi-Toeplitz matrix into two parts, namely the analysis for the boundary condition independent spectrum and the analysis for the boundary condition dependent spectrum. The boundary condition independent spectrum is the same as the pure Toeplitz matrix spectrum. Algorithms for computing both parts of the spectrum are presented. Examples are used to demonstrate the utility of the algorithms, to present some interesting spectra, and to point out some of the numerical difficulties encountered when conventional matrix eigenvalue routines are employed for non-normal matrices of large order. The analysis for the Toeplitz spectrum also leads to a diagonal similarity transformation that improves conventional numerical eigenvalue computations. Finally, the algorithm for the asymptotic spectrum is extended to the Toeplitz generalized eigenvalue problem which occurs, for example, in the stability of Pade type difference approximations to differential equations.

Beam, Richard M.↗
NASA NTRS1992

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier-Stokes-Boussinesq equations governing the transport of momentum, mass and heat in a nonisothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a nonorthogonal coordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.

Zhang, Yiqiang↗
NASA NTRS2018

The Effect of Forward-Facing Steps on Stationary Crossflow Instability Growth and Breakdown

The effect of a forward-facing step on stationary crossflow transition was studied using standard stereo particle image velocimetry (PIV) and time-resolved PIV. Step heights ranging from 53 to 71% of the boundary-layer thickness were studied in detail. The steps above a critical step height of approximately 60% of the boundary-layer thickness had a significant impact on the stationary crossflow growth downstream of the step. For the critical cases, the stationary crossflow amplitude grew suddenly downstream of the step, decayed for a short region, then grew again. The adverse pressure gradient upstream of the step resulted in a region of crossflow reversal. A secondary set of vortices, rotating in the opposite direction to the primary vortices, developed underneath the uplifted primary vortices. The wall-normal velocity disturbance (V' ) created by these secondary vortices impacted the step, and is believed to feed into the strong vortex that developed downstream of the step. A large but very short negative crossflow region formed for a short region downstream of the step due to a sharp inboard curvature of the streamlines near the wall. For the larger step height cases, a crossflow-reversal region formed just downstream of the strong negative crossflow region. This crossflow reversal region is believed to play an important role in the growth of the stationary crossflow vortices downstream of the step, and may be a good indication of the critical forward-facing step height.

Eppink, Jenna L.↗
NASA NTRS1983

An integral analysis of transonic normal shock wave/turbulent boundary layer interactions in internal flow

An approximate integral viscous-inviscid interaction method is presented for calculating the development of a turbulent boundary layer subjected to a normal shock wave induced adverse pressure gradient in an internal axisymmetric flow. The inflow conditions and the downstream pressure are provided for the computation. In the supersonic region of shock pressure rise, the Prandtl-Meyer function is used to couple the viscous and inviscid flows. An analytical model for the coupling process is postulated and appropriate equations are defined. Downstream of the sonic point, one-dimensional inviscid flow is assumed for coupling with the viscous flow. The turbulent boundary layer is calculated using Green's integral lag-entrainment method. Comparisons of the solutions with the experimental data are made for interactions which are unseparated, near separation and separated. For comparison purposes, solutions to the time-dependent, mass-averaged, Navier-Stokes equations incorporating a two-equation, Wilcox-Rubesin turbulence model are also shown. The computed results from the integral method show good agreement with experimental data for unseparated interactions and reasonable agreement with the trend of the viscous effects when the interaction becomes increasingly separated.

Om, D.↗
NASA NTRS1984

Shape design sensitivity analysis of built-up structures

Selection of the best shape of a fillet in a tension bar such that no yielding occurs has long attracted the attention of engineers. Dimensions and notations for the bar and fillet are shown. With symmetry, only the upper half of the bar is considered. The optimal design problem is to find a boundary shape to minimize the total area of the fillet such that no yielding occurs. The classical boundary value problem is reduced to a variational or energy related problem which not only has excellent properties of existence and uniqueness but also provides the mathematical foundation for finite element analysis. The variational formulation may be viewed as the principle of virtual work and the finite element method as a application of the Galerkin method to the variational equation for approximate solution of the boundary value problem.

Choi, K. K.↗
NASA NTRS2005

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗
NASA NTRS1981

A quasi-three-dimensional blade surface boundary layer analysis for rotating blade rows

A quasi-three-dimensional approximation has been developed for a blade boundary layer which involves the calculation of the effect of nonzero pressure gradients, turbulent flow, and blade twist, but includes only a simple coupling between streamlines. The resulting set of equations is solved using Keller's box scheme. The solution scheme is checked against available incompressible flow solutions and then applied to a NASA low aspect ratio transonic compressor stage for which extensive experimental and computational data are available. It is found that the three-dimensional boundary layer separates significantly sooner and has a much larger influence on rotor performance than would be expected from a two-dimensional analysis.

Thompkins, W. T., Jr.↗
NASA NTRS1990

Resolution of the 1D regularized Burgers equation using a spatial wavelet approximation

The Burgers equation with a small viscosity term, initial and periodic boundary conditions is resolved using a spatial approximation constructed from an orthonormal basis of wavelets. The algorithm is directly derived from the notions of multiresolution analysis and tree algorithms. Before the numerical algorithm is described these notions are first recalled. The method uses extensively the localization properties of the wavelets in the physical and Fourier spaces. Moreover, the authors take advantage of the fact that the involved linear operators have constant coefficients. Finally, the algorithm can be considered as a time marching version of the tree algorithm. The most important point is that an adaptive version of the algorithm exists: it allows one to reduce in a significant way the number of degrees of freedom required for a good computation of the solution. Numerical results and description of the different elements of the algorithm are provided in combination with different mathematical comments on the method and some comparison with more classical numerical algorithms.

Liandrat, J.↗
NASA NTRS1978

Summary of experimentally determined facts concerning the behavior of the boundary layer and performance of boundary layer measurements

A summary report of boundary layer studies is presented. Preliminary results of experimental measurements show that: (1) A very thin layer (approximately 0.4 mm) of the boundary layer seems to be accelerated; (2) the static pressure of the outer flow does not remain exactly constant through the boundary layer; and (3) an oncoming boundary layer which is already turbulent at the suction point can again become laminar behind this point without being completely sucked off.

Vanness, W.↗
NASA NTRS1994

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗
NASA NTRS2016

Experimental Investigation of Subsonic Turbulent Boundary Layer Flow Over a Wall-Mounted Axisymmetric Hill

An important goal for modern fluid mechanics experiments is to provide datasets which present a challenge for Computational Fluid Dynamics simulations to reproduce. Such "CFD validation experiments" should be well-characterized and well-documented, and should investigate flows which are difficult for CFD to calculate. It is also often convenient for the experiment to be challenging for CFD in some aspects while simple in others. This report is part of the continuing documentation of a series of experiments conducted to characterize the flow around an axisymmetric, modified-cosine-shaped, wall-mounted hill named "FAITH" (Fundamental Aero Investigates The Hill). Computation of this flow is easy in some ways - subsonic flow over a simple shape - while being complex in others - separated flow and boundary layer interactions. The primary set of experiments were performed on a 15.2 cm high, 45.7 cm base diameter machined aluminum model that was tested at mean speeds of 50 m/s (Reynolds Number based on height = 500,000). The ratio of model height to boundary later height was approximately 3. The flow was characterized using surface oil flow visualization, Cobra probe to determine point-wise steady and unsteady 3D velocities, Particle Image Velocimetry (PIV) to determine 3D velocities and turbulence statistics along specified planes, Pressure Sensitive Paint (PSP) to determine mean surface pressures, and Fringe Imaging Skin Friction (FISF) to determine surface skin friction magnitude and direction. A set of pathfinder experiments were also performed in a water channel on a smaller scale (5.1 cm high, 15.2 cm base diameter) sintered nylon model. The water channel test was conducted at a mean test section speed of 3 cm/s (Reynolds Number of 1500), but at the same ratio of model height to boundary layer thickness. Dye injection from both the model and an upstream rake was used to visualize the flow. This report summarizes the experimental set-up, techniques used, and data acquired. It also describes some details of the dataset that is being constructed for use by other researchers, especially the CFD community.

fluid machanics cfd validation wind tunnel↗
NASA NTRS1993

A three-dimensional pressure flux-split RNS application to sub/supersonic flow in inlets and ducts

The reduced Navier-Stokes (RNS) formulation is combined with a pressure based flux split procedure for the computational analysis of three dimensional flow. Only lowest order diffusion, that is required to satisfy the no-slip boundary condition is retained in this approximation. The governing equations are appropriately differenced such that the physical boundary conditions are closely coupled and combined to form a closed discrete system. No numerical or characteristic type boundary conditions are required. A global pressure relaxation procedure, i.e. multi-sweep PNS, is considered. A sparse matrix direct solver is applied to the crossplane solution. At the outflow boundary only the pressure or pressure gradient is prescribed. The resulting formulation has been applied to the solution of three-dimensional subsonic/supersonic internal flows in inlets and ducts.

Khosla, P. K.↗
NASA NTRS1974

Analysis of laminated, composite, circular cylindrical shells with general boundary conditions

This report develops: (1) a refined approximate theory for the static and dynamic analyses of finite, laminated, composite, circular cylindrical shells with general boundary conditions; (2) an exact three-dimensional analysis of simply supported, laminated, composite, circular cylindrical shells, and (3) a thin-shell theory for laminated, composite, circular cylindrical shells. In the refined approximate theory the displacements are assumed piecewise linear across the thickness and the effects of transverse shear deformations and transverse normal stress are included. A variational approach is followed to obtain the governing differential equations and boundary conditions. A general solution of the governing differential equations is also presented. The results obtained by using the refined approximate theory and the thin-shell theory are compared with the exact results for the case of free vibrations of simply supported, laminated, composite, circular cylindrical shells. The refined approximate theory is very accurate, even for thick shells with short nodal distances, whereas thin-shell theory is reasonably accurate only for thin shells at moderate nodal distances and wave number less than 2.

Srinivas, S.↗
DOE OSTI2023

A Fourth-Order Embedded Boundary Finite Volume Method for the Unsteady Stokes Equations with Complex Geometries

A fourth-order finite volume embedded boundary (EB) method is presented for the unsteady Stokes equations. The algorithm represents complex geometries on a Cartesian grid using EB, employing a technique to mitigate the ``small cut-cell"" problem without mesh modifications, cell merging, or state redistribution. Spatial discretizations are based on a weighted least-squares technique that has been extended to fourth-order operators and boundary conditions, including an approximate projection to enforce the divergence-free constraint. Solutions are advanced in time using a fourth-order additive implicit-explicit Runge-Kutta method, with the viscous and source terms treated implicitly and explicitly, respectively. Formal accuracy of the method is demonstrated with several grid convergence studies, and results are shown for an application with a complex bio-inspired material. In conclusion, the developed method achieves fourth-order accuracy and is stable despite the pervasive small cells arising from complex geometries.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1956

The Compressible Laminar Boundary Layer with Heat Transfer and Arbitrary Pressure Gradient

An approximate method for the calculation of the compressible laminar boundary layer with heat transfer and arbitrary pressure gradient, based on Thwaites' correlation concept, is presented. With the definition of dimensionless shear and heat-transfer parameters and an assumed correlation of these parameters in terms of a momentum parameter, a complete system of relations for calculating skin friction and heat transfer results. Knowledge of velocity or temperature profiles is not necessary in using this calculation method. When the method is applied to a convergent-divergent, axially symmetric rocket nozzle, it shows that high rates of heat transfer are obtained at the initial stagnation point and at the throat of the nozzle. Also indicated are negative displacement thicknesses in the convergent portion of the nozzle; these occur because of the high density within the lower portions of the cooled boundary layer. (author)

Cohen, Clarence B↗
NASA NTRS2003

A Reduced-Order Model for Efficient Simulation of Synthetic Jet Actuators

A new reduced-order model of multidimensional synthetic jet actuators that combines the accuracy and conservation properties of full numerical simulation methods with the efficiency of simplified zero-order models is proposed. The multidimensional actuator is simulated by solving the time-dependent compressible quasi-1-D Euler equations, while the diaphragm is modeled as a moving boundary. The governing equations are approximated with a fourth-order finite difference scheme on a moving mesh such that one of the mesh boundaries coincides with the diaphragm. The reduced-order model of the actuator has several advantages. In contrast to the 3-D models, this approach provides conservation of mass, momentum, and energy. Furthermore, the new method is computationally much more efficient than the multidimensional Navier-Stokes simulation of the actuator cavity flow, while providing practically the same accuracy in the exterior flowfield. The most distinctive feature of the present model is its ability to predict the resonance characteristics of synthetic jet actuators; this is not practical when using the 3-D models because of the computational cost involved. Numerical results demonstrating the accuracy of the new reduced-order model and its limitations are presented.

Yamaleev, Nail K.↗
DOE OSTI2023

Interface flux recovery framework for constructing partitioned heterogeneous time-integration methods

A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.

97 MATHEMATICS AND COMPUTING↗