SEARCH · Engineering Papers

Results for “Boundary value problem”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Thrust vectoring effect on time-optimal 90 degrees angle of attack pitch up maneuvers of a high alpha fighter aircraft

Using thrust vectoring, the problem of pitching up a high alpha fighter aircraft to 90 degrees angle of attack in minimum time is considered. Pontryagin's maximum principle is used together with a two-point boundary value numerical algorithm to derive open-loop controls for various parameterized limits on thrust vectoring angles. Without thrust vectoring, a high alpha fighter with initial conditions 0.6 Mach number and 15,000 feet requires 7.0 seconds to pitch-up to 90 degrees. However, with 20-degree angle thrust vectoring, the pitch-up takes only 1.5 seconds with a final Mach number of 0.47. Also, the gain in altitude is less than 100 feet. More important, additional thrust vectoring yields litle benefit beyond the 20-degree angles. This paper describes control solutions in detail for thrust vectoring angles between 0 and 90 degrees.

Stalford, Harold↗

NASA NTRS1971

The problem of uniqueness in thermosphere dynamics

The unique solution of tidal wave propagation within the thermosphere depends on the boundary conditions of the model. It is shown that the radiation condition leads to such a unique solution. Any other boundary values give rise to deviations from the physical solution with unrealistic physical parameters in the environment of the boundaries. The thickness of these boundary layers with unrealistic solutions is a few scale heights below the upper boundary of the model and a few tens of scale heights above the lower boundary of the model.

Volland, H.↗

NASA NTRS1993

Improved Finite Element Modeling of the Turbofan Engine Inlet Radiation Problem

Improvements have been made in the finite element model of the acoustic radiated field from a turbofan engine inlet in the presence of a mean flow. The problem of acoustic radiation from a turbofan engine inlet is difficult to model numerically because of the large domain and high frequencies involved. A numerical model with conventional finite elements in the near field and wave envelope elements in the far field has been constructed. By employing an irrotational mean flow assumption, both the mean flow and the acoustic perturbation problem have been posed in an axisymmetric formulation in terms of the velocity potential; thereby minimizing computer storage and time requirements. The finite element mesh has been altered in search of an improved solution. The mean flow problem has been reformulated with new boundary conditions to make it theoretically rigorous. The sound source at the fan face has been modeled as a combination of positive and negative propagating duct eigenfunctions. Therefore, a finite element duct eigenvalue problem has been solved on the fan face and the resulting modal matrix has been used to implement a source boundary condition on the fan face in the acoustic radiation problem. In the post processing of the solution, the acoustic pressure has been evaluated at Gauss points inside the elements and the nodal pressure values have been interpolated from them. This has significantly improved the results. The effect of the geometric position of the transition circle between conventional finite elements and wave envelope elements has been studied and it has been found that the transition can be made nearer to the inlet than previously assumed.

Roy, Indranil Danda↗

NASA NTRS1987

Nonlinear force-free magnetic fields: Calculation and applicatin to astrophysics

The problem concerned in this work is that of calculating magnetic field configurations in which the Lorentz force (vector)j x (vector)B is everywhere zero, subject to specified boundary conditions. The magnetic field is represented in terms of Clebsch variables in the form (vector)B = del alpha x del beta. These variables are constant on any field line. The most appropriate choice of boundary conditions is to specify the values of alpha and beta on the bounding surface. It is proposed that the field lines move in the direction of local Lorentz force and relax towards a force-free field configuration. This concept leads to an iteration procedure for modifying the variables alpha and beta that tends asymptotically towards the force-free state. This method is first applied to a simple problem in two rectangular dimensions; the calculation shows that the convergence of magnetic field energy to a minimum state (force-free) is close to exponential. This method is then applied to study some astrophysical force-free magnetic fields, such as the structures and evolution of magnetic fields of rotating sunspots and accretion disks. The implication of the results, as related to the mechanisms of solar flares, extragalactic radio sources and radio jets, are discussed.

Yang, Wei-Hong↗

NASA NTRS2011

A Source-Term Based Boundary Layer Bleed/Effusion Model for Passive Shock Control

A modeling framework for boundary layer effusion has been developed based on the use of source (or sink) terms instead of the usual practice of specifying bleed directly as a boundary condition. This framework allows the surface boundary condition (i.e. isothermal wall, adiabatic wall, slip wall, etc.) to remain unaltered in the presence of bleed. This approach also lends itself to easily permit the addition of empirical models for second order effects that are not easily accounted for by simply defining effective transpiration values. Two effusion models formulated for supersonic flows have been implemented into this framework; the Doerffer/Bohning law and the Slater formulation. These models were applied to unit problems that contain key aspects of the flow physics applicable to bleed systems designed for hypersonic air-breathing propulsion systems. The ability of each model to predict bulk bleed properties was assessed, as well as the response of the boundary layer as it passes through and downstream of a porous bleed system. The model assessment was performed with and without the presence of shock waves. Three-dimensional CFD simulations that included the geometric details of the porous plate bleed systems were also carried out to supplement the experimental data, and provide additional insights into the bleed flow physics. Overall, both bleed formulations fared well for the tests performed in this study. However, the sample of test problems considered in this effort was not large enough to permit a comprehensive validation of the models.

Baurle, Robert A.↗

NASA NTRS1987

Linear stability analysis of three-dimensional compressible boundary layers

A compressible stability analysis computer code is developed. The code uses a matrix finite-difference method for local eigenvale solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used inital-value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed that may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black-box stability code. Sample stability calculations are presented for the boundary layer profiles of an LFC swept wing.

Malik, Mujeeb R.↗

NASA NTRS1973

High-speed leading edge problem.

The sharp leading edge problem has been studied for both monatomic and diatomic gases using the Boltzmann equation with the Bhatnagar-Gross-Krook type models as the governing equation and the discrete ordinate method with a closed-boundary value approach as a tool. Plate length relative to the freestream mean free path is taken to be 52. The gas-surface interaction law is assumed to be diffuse reflection. The local distribution functions of molecular velocities and internal energies (for the diatomic gas) for the entire flowfield have been calculated for a freestream Mach number of 6.1. Comparisons are made between the calculated results and experimental data.

Huang, A. B.↗

NASA NTRS1974

A finite-element method for lifting surfaces in steady incompressible subsonic flow

The problem of potential steady subsonic flow for lifting surfaces is considered. This problem requires the solution of an integral equation relating the value of the potential discontinuity on the lifting surface and its wake to the values of the normal derivative of the potential which are known from the boundary conditions. The lifting surface is divided into small (quadrilateral hyperboloidal) surface elements, which are described in terms of the Cartesian components of the four corner points. The values of the potential discontinuity and the normal derivative of the potential are assumed to be constant within each element and equal to their values at the centroids of the elements. This yields a set of linear algebraic equations. Numerical results are in good agreement with existing ones.

Morino, L.↗

NASA NTRS2011

Inflow/Outflow Boundary Conditions with Application to FUN3D

Several boundary conditions that allow subsonic and supersonic flow into and out of the computational domain are discussed. These boundary conditions are demonstrated in the FUN3D computational fluid dynamics (CFD) code which solves the three-dimensional Navier-Stokes equations on unstructured computational meshes. The boundary conditions are enforced through determination of the flux contribution at the boundary to the solution residual. The boundary conditions are implemented in an implicit form where the Jacobian contribution of the boundary condition is included and is exact. All of the flows are governed by the calorically perfect gas thermodynamic equations. Three problems are used to assess these boundary conditions. Solution residual convergence to machine zero precision occurred for all cases. The converged solution boundary state is compared with the requested boundary state for several levels of mesh densities. The boundary values converged to the requested boundary condition with approximately second-order accuracy for all of the cases.

Carlson, Jan-Renee↗

NASA NTRS1994

On the nonlinear stability of viscous modes within the Rayleigh problem on an infinite flat plate

The stability has been investigated of the unsteady flow past an infinite flat plate when it is moved impulsively from rest, in its own plane. For small times the instantaneous stability of the flow depends on the linearized equations of motion which reduce in this problem to the Orr-Sommerfeld equation. It is known that the flow for certain values of Reynolds number, frequency and wave number is unstable to Tollmien-Schlichting waves, as in the case of the Blasius boundary layer flow past a flat plate. With increase in time, the unstable waves only undergo growth for a finite time interval, and this growth rate is itself a function of time. The influence of finite amplitude effects is studied by solving the full Navier-Stokes equations. It is found that the stability characteristics are markedly changed both by the consideration of the time evolution of the flow, and by the introduction of finite amplitude effects.

Webb, J. C.↗

NASA NTRS1974

An eddy-viscosity treatment of the unsteady turbulent boundary layer on a flat plate in an expansion tube

An analysis is presented for the relaxation of a turbulent boundary layer on a semiinfinite flat plate after passage of a shock wave and a trailing driver gas-driven gas interface. The problem has special application to expansion tube flows. The flow-governing equations have been transformed into the Lamcrocco variables. The numerical results indicate that a fully turbulent boundary layer relaxes faster to the final steady-state values of heat transfer and skin-friction than a fully laminar boundary layer.

Gupta, R. N.↗

NASA NTRS1986

The effect of a small initial curvature on the free vibration of clamped, rectangular plates

An analytical method of obtaining the natural frequencies and mode shapes of clamped, rectangular plates having a small initial curvature is presented. Specifically, the singular perturbation technique is used to reduce the fourth-order plate vibration problem to the simpler membrane problem with modified boundary conditions that account for the bending effects. The eigenfrequencies for plates with inverse aspect ratios varying between 0.1 and 1.0 and for the dimensionless normal prestress between 0.1 and 1.0 have been presented for values of epsilon, the normalized bending rigidity, ranging between 0.0010 and 0.2500. It is established that a small initial curvature has no effect on the frequency of vibration of the plate. However, its effect is manifested in the eigenmodes.

Adeniji-Fashola, A. A.↗

NASA NTRS1989

Shear buckling of specially orthotropic plates with centrally located cutouts

The shear loading of a composite rectangular plate with a centrally located circular cutout was analyzed in order to predict the buckling load of the plate. The first step in this analysis is to calculare the in-plane stress distribution of the plate prior to buckling. This problem can be solved using finite element methods, but a quicker and more cost efficient method is desired. The method chosen to determine the pre-buckling stresses is that of boundary collocation using complex variables. The real valued force functions are written in terms of two complex valued functions, each of which is a function of a different complex variable. The solution of the generalized biharmonic function is a superposition of functions of these two complex variables. For this analysis, the force functions are represented with a Laurent series. The constants in these series are found using boundary collocation. Two force equations for complex variables are satisfied over the plate boundaries to obtain these constants. The stresses in the plate are then known for any particular locations on the plate given the applied shear on the edges of the plate.

Klang, Eric C.↗

NASA NTRS1993

Stability of the Halley cometosheath with resistivity and plasma motion

The MHD stability of the cometary inner sheath determined by the balance between the inward Lorentz body force and the outward ion-neutral drag force is investigated by numerically solving the wave equations which include resistivity, plasma motion and plasma pressure with the help of two-point boundary value method. The eigenvalues and the eigenfunctions are obtained numerically by treating the cometary inner sheath as a layer of finite thickness, bounded by the contact surface, that is, the diamagnetic cavity boundary. To gain insight into the problem, certain limiting cases of the wave equations are also discussed. The diamagnetic cavity boundary and the adjacent layer of about 100-km thickness of Comet Halley is found to be unstable. The effects of finite plasma pressure, dissociative recombination, mass loading due to photoionization, resistivity, and plasma motion are found to be stabilizing but are unable to quench the instability completely. Motion of the Halley ionopause has been confirmed by observations.

Srivastava, Krishna M.↗

NASA NTRS1990

Vortex instabilities in 3D boundary layers: The relationship between Goertler and crossflow vortices

The inviscid and viscous stability problems are addressed for a boundary layer which can support both Goertler and crossflow vortices. The change in structure of Goertler vortices is found when the parameter representing the degree of three-dimensionality of the basic boundary layer flow under consideration is increased. It is shown that crossflow vortices emerge naturally as this parameter is increased and ultimately become the only possible vortex instability of the flow. It is shown conclusively that at sufficiently large values of the crossflow there are no unstable Goertler vortices present in a boundary layer which, in the zero crossflow case, is centrifugally unstable. The results suggest that in many practical applications Goertler vortices cannot be a cause of transition because they are destroyed by the 3-D nature of the basic state. In swept wing flows the Goertler mechanism is probably not present for typical angles of sweep of about 20 degrees. Some discussion of the receptivity problem for vortex instabilities in weakly 3-D boundary layers is given; it is shown that inviscid modes have a coupling coefficient marginally smaller than those of the fastest growing viscous modes discussed recently by Denier, Hall, and Seddougui (1990). However the fact that the growth rates of the inviscid modes are the largest in most situations means that they are probably the most likely source of transition.

Bassom, Andrew↗

NASA NTRS1986

Transformation of two and three-dimensional regions by elliptic systems

Efforts in transferring computational work from the LRC computer to the IRIS Graphics Workstation at MSU are reported and the computation of a conservative solution of a simple hyperbolic equation on an overlapping grid is discussed. Several conclusions concerning computations on overlapping grids are apparent. Problems only occur when there is a major difference in grid spacing on the individual component grids. In the case of hyperbolic equations, it is necessary that both interpolation and extrapolation be applied at the grid boundaries. When interpolated values are used at outflow boundary points, excessive oscillations in the numerical solution may be the result. The same conclusions would be valid for more complicated systems of hyperbolic equations such as the Euler equations for inviscid flow. Some of the solution values would be extrapolated at the overlap boundary, the exact number depending on the number of characteristics pointing out of the overlap region. It is also possible that similar boundary conditions may be needed for some parabolic equations such as high Reynolds number viscous flow equations. Efforts were also expended on the development of three-dimensional conservative interpolation procedures. Finally, the investigation of grid smoothing procedures were initiated during this reporting period. It was decided that the first grid smoothing algorithms will be based on the concepts of variational grid generation.

Mastin, C. W.↗

NASA NTRS1990

Time development of a perturbed-spherical nucleus in a pure supercooled liquid. I - Power-law growth of morphological instabilities. II - Nonlinear development

The linear growth stage of the morphological instabilities developing at the liquid-solid interface during crystal growth from an axisymmetric spherical nucleus is analyzed. The corresponding growth rate parameters are calculated numerically, and it is shown that morphological instabilities for free growth evolve according to a power law in agreement with WKB results and contrary to an exponential law found in a quasi-stationary approximation. The time evolution of an arbitrary perturbed liquid-solid interface from a linear into the nonlinear stage is studied. The initial perturbations include the eigenfunctions for the linear problem, localized perturbations of the grain-boundary type, and a stochastic noise. It is shown that the perturbations grow and spread in a wavelike manner. The formation of a predendritic growth stage is characterized by establishment of constant values of tip radius and velocity.

Pines, V.↗

NASA NTRS1996

Thermocapillary Bubble Migration: Thermal Boundary Layers for Large Marangoni Numbers

The migration of an isolated gas bubble in an immiscible liquid possessing a temperature gradient is analyzed in the absence of gravity. The driving force for the bubble motion is the shear stress at the interface which is a consequence of the temperature dependence of the surface tension. The analysis is performed under conditions for which the Marangoni number is large, i.e. energy is transferred predominantly by convection. Velocity fields in the limit of both small and large Reynolds numbers are used. The thermal problem is treated by standard boundary layer theory. The outer temperature field is obtained in the vicinity of the bubble. A similarity solution is obtained for the inner temperature field. For both small and large Reynolds numbers, the asymptotic values of the scaled migration velocity of the bubble in the limit of large Marangoni numbers are calculated. The results show that the migration velocity has the same scaling for both low and large Reynolds numbers, but with a different coefficient. Higher order thermal boundary layers are analyzed for the large Reynolds number flow field and the higher order corrections to the migration velocity are obtained. Results are also presented for the momentum boundary layer and the thermal wake behind the bubble, for large Reynolds number conditions.

Balasubramaniam, R.↗