Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary value problems”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Objective and algorithm considerations when optimizing the number and placement of turbines in a wind power plant

Abstract. Optimizing turbine layout is a challenging problem that has been extensively researched in the literature. However, optimizing the number of turbines within a given boundary has not been studied as extensively and is a difficult problem because it introduces discrete design variables and a discontinuous design space. An essential step in performing wind power plant layout optimization is to define the objective function, or value, that is used to express what is valuable to a wind power plant developer, such as annual energy production, cost of energy, or profit. In this paper, we demonstrate the importance of selecting the appropriate objective function when optimizing a wind power plant in a land-constrained site. We optimized several different wind power plants with different wind resources and boundary sizes. Results show that the optimal number of turbines varies drastically depending on the objective function. For a simple, one-dimensional, land-based scenario, we found that a wind power plant optimized for minimal cost of energy produced just 72 % of the profit compared to the wind power plant optimized for maximum profit, which corresponded to a loss of about USD 2 million each year. This paper also compares the performance of several different optimization algorithms, including a novel repeated-sweep algorithm that we developed. We found that the performance of each algorithm depended on the number of design variables in the problem as well as the objective function.

17 WIND ENERGY↗
NASA NTRS1976

A new method for solving problems of sound radiation from a duct for the case of low Mach number

A method for rapidly solving duct end-plane impedance problems is presented which is based on deriving a form of the Helmholtz integral formula expressing the normal derivative of the acoustic pressure at a boundary point. A set of integral expressions is obtained by satisfying the existing boundary conditions and then solved using a general collocation method. The results obtained for the Levine-Schwinger (1948) problem agree well with the exact values. In the case of a duct flow having a temperature mismatch with the surroundings, the results show an increase in the magnitude of the reflection coefficient for the case of a cold core, and a decrease in its value for the case of a hot core, as compared with the case of exhaust into a uniform medium. A phase change of between pi/2 and pi is nearly always achieved, which indicates the tendency towards maintaining a constant pressure at the exit plane. The method can easily be extended to arbitrary duct shapes as long as they are axisymmetric.

Wahbah, M. M.↗
NASA NTRS1978

Numerical spatial marching techniques for estimating duct attenuation and source pressure profiles

A numerical method was developed that could predict the pressure distribution of a ducted source from far field pressure inputs. Using an initial value formulation, the two-dimensional homogeneous Helmholtz wave equation (no steady flow) was solved using explicit marching techniques. The Von Neumann method was used to develop relationships which describe how sound frequency and grid spacing effect numerical stability. At the present time, stability considerations limit the approach to high frequency sound. Sample calculations for both hard and soft wall ducts compare favorably to known boundary value solutions. In addition, assuming that reflections in the duct are small, this initial value approach was successfully used to determine the attenuation of a straight soft wall duct. Compared to conventional finite difference or finite element boundary value approaches, the numerical marching technique is orders of magnitude shorter in computation time and required computer storage and can be easily employed in problems involving high frequency sound.

Baumeister, K. J.↗
NASA NTRS1978

Numerical spatial marching techniques for estimating duct attenuation and source pressure profiles

A numerical method is developed that could predict the pressure distribution of a ducted source from far-field pressure inputs. Using an initial value formulation, the two-dimensional homogeneous Helmholtz wave equation (no steady flow) is solved using explicit marching techniques. The Von Neumann method is used to develop relationships which describe how sound frequency and grid spacing effect numerical stability. At the present time, stability considerations limit the approach to high frequency sound. Sample calculations for both hard and soft wall ducts compare favorably to known boundary value solutions. In addition, assuming that reflections in the duct are small, this initial value approach is successfully used to determine the attenuation of a straight soft wall duct. Compared to conventional finite difference or finite element boundary value approaches, the numerical marching technique is orders of magnitude shorter in computation time and required computer storage and can be easily employed in problems involving high frequency sound.

Baumeister, K. J.↗
NASA NTRS1991

Generating Three-Dimensional Grids About Anything

Three-Dimensional Grids About Anything by Poisson's Equation (3DGRAPE) computer program designed to make computational grids in or about almost any shape. Generated by solution of Poisson's differential equations in three dimensions. Program automatically finds its own values for inhomogeneous terms giving near-orthogonality and controlled grid-cell height at boundaries. Grids generated applied to both viscous and inviscid aerodynamic problems, and to problems in other areas of fluid dynamics. Written in 100 percent FORTRAN 77.

Sorenson, Reese L.↗
NASA NTRS1991

A study of three dimensional turbulent boundary layer separation and vortex flow control using the reduced Navier Stokes equations

A reduced Navier Stokes (RNS) initial value space marching solution technique was applied to vortex generator and separated flow problems and demonstrated good predictions of the engine face flow field. This RNS solution technique using FLARE approximations can adequately describe the topological and topographical structure flow separation associated with vortex liftoff, and this conclusion led to the concept of a subclass of separations which can be called vorticity separations: separations dominated by the transport of vorticity. Adequate near wall resolution of vorticity separations appears necessary for good predictions of these flows.

Anderson, Bernhard H.↗
NASA NTRS1998

Theory of the Shape of the Heliospheric Termination Shock: Axial Symmetry

The heliospheric termination shock must exhibit asymmetry in its shape, due in part to internal latitude variations in the solar wind, and in part to the special directions defined by the external interstellar flow and/or the galactic magnetic field. This asymmetry shows up naturally in numerical simulations of the interaction between the heliosphere and local interstellar medium. However, to date only one analytical treatment of the theory has appeared, which discussed the modification of the shock shape due to solar wind latitude variations in the presence of spherically symmetric outer boundary conditions. In the present report, we discuss an extension and generalization of this theory to an axially symmetric gasdynamic system, in which departure from spherical symmetry may be due either to internal solar wind variations or to the directional properties of the external interstellar medium. It is shown that for steady flow the post-shock region is characterized by an infinite set of quantities that are conserved along streamlines; among these invariants are the stagnation pressure and a quantity closely related to vorticity. Moreover, for a given latitude profile of the (supersonic) solar wind, the geometry of the termination shock uniquely determines the valued of these invariants at the points where the streamlines emerge from the shock. A self-consistent solution of the problem thus requires a match between the shock geometry and distant boundary conditions. The theory is applied to give estimates of termination-shock asymmetry for various internal and external conditions.

Barnes, Aaron↗
DOE OSTI2016

Study of a mixed dispersal population dynamics model

In this study, we consider a mixed dispersal model with periodic and Dirichlet boundary conditions and its corresponding linear eigenvalue problem. This model describes the time evolution of a population which disperses both locally and non-locally. We investigate how long time dynamics depend on the parameter values. Furthermore, we study the minimization of the principal eigenvalue under the constraints that the resource function is bounded from above and below, and with a fixed total integral. Biologically, this minimization problem is motivated by the question of determining the optimal spatial arrangement of favorable and unfavorable regions for the species to die out more slowly or survive more easily. Our numerical simulations indicate that the optimal favorable region tends to be a simply-connected domain. Numerous results are shown to demonstrate various scenarios of optimal favorable regions for periodic and Dirichlet boundary conditions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

A method for determining acoustic-liner admittance in ducts with sheared flow in two-cross-sectional directions

A method is developed for determining the acoustic admittance of a test liner installed in the wall of a grazing flow impedance tube. The mean flow is permitted flow gradients in both cross-sectional directions of the tube. The unknown admittance value is obtained by solving an eigenvalue problem. This eigenvalue problem results from the application of the finite-element method to the partial differential equation and boundary conditions governing the acoustic field. The credibility of the method is established by comparing results with exact solutions obtained for a constant mean-flow profile and with previous results for cases involving shear in only one cross-sectional direction. Excellent comparisons are obtained in both cases. The analysis is used in conjunction with a limited amount of experimental data and shows that the flow must be accurately modeled in order to determine the acoustic-liner properties.

Watson, W. R.↗
NASA NTRS1987

How good is the impedance boundary condition?

The impedance boundary condition (IBC) is often used in scattering problems involving material-coated conducting bodies. It is shown that for some commonly encountered coating configurations, the value of the impedance varies significantly as functions of the incident angle and polarization. Hence, the use of IBC in a rigorously formulated problem may affect the accuracy of the final solution.

Lee, Shung-Wu↗
NASA NTRS1999

Receptivity and Bypass Dynamics

Problems concerning laminar-turbulent transition are addressed by solving a series of initial value problems. Solutions to the temporal, initial-value problem .with an inhomogeneous forcing term imposed upon the flow are sought. It is shown that: (1) A transient disturbance lying located outside of the boundary layer can lead to the growth of an unstable Tollmein-Schlicting wave; (2) A resonance with the continuous spectrum may provide a mechanism for bypass transition; and (3) The continuum modes of a disturbance feed directly into the Tollmein-Schlicting wave downstream through non-parallel effects.

Lasseigne, D. G.↗
NASA NTRS2012

Correction of Excessive Precipitation over Steep Mountains in a General Circulation Model (GCM)

Excessive precipitation over steep and high mountains (EPSM) is a well-known problem in GCMs and regional climate models even at a resolution as high as 19km. The affected regions include the Andes, the Himalayas, Sierra Madre, New Guinea and others. This problem also shows up in some data assimilation products. Among the possible causes investigated in this study, we found that the most important one, by far, is a missing upward transport of heat out of the boundary layer due to the vertical circulations forced by the daytime subgrid-scale upslope winds, which in turn is forced by heated boundary layer on the slopes. These upslope winds are associated with large subgrid-scale topographic variance, which is found over steep mountains. Without such subgrid-scale heat ventilation, the resolvable-scale upslope flow in the boundary layer generated by surface sensible heat flux along the mountain slopes is excessive. Such an excessive resolvable-scale upslope flow in the boundary layer combined with the high moisture content in the boundary layer results in excessive moisture transport toward mountaintops, which in turn gives rise to excessive precipitation over the affected regions. We have parameterized the effects of subgrid-scale heated-slope-induced vertical circulation (SHVC) by removing heat from the boundary layer and depositing it in the layers higher up when topographic variance exceeds a critical value. Test results using NASA/Goddard's GEOS-5 GCM have shown that the EPSM problem is largely solved.

Chao, Winston C.↗
NASA NTRS1991

A study of three dimensional turbulent boundary layer separation and vortex flow control using the reduced Navier Stokes equations

A reduced Navier-Stokes (RNS) initial value space marching solution technique was applied to vortex generator and separated flow problems and demonstrated good predictions of the engine face flow field. This RNS solution technique using FLARE approximations can adequately describe the topological and topographical structure flow separation associated with vortex liftoff, and this conclusion led to the concept of a subclass of separations which can be called vorticity separations: separations dominated by the transport of vorticity. Adequate near wall resolution of vorticity separations appears necessary for good predictions of these flows.

Anderson, Bernhard H.↗
NASA NTRS1966

Numerical Studies of the Non-linear Vlasov Equation

The subject of this dissertation is the numerical integration of the initial-value problem for the non-linear Vlasov equation. The Vlasov equation is used to describe the dynamics of a "collisionless", one-dimensional, classical electron gas confined between two perfectly reflecting boundaries. Only the long-range Coulomb interactions of the electrons are considered; effects associated with the discrete structure are neglected. The numerical results obtained for non-linear Landau damping compare well with similar results obtained by Knorr. A general statement of the results on stable initial conditions is: As the degree of non-linearity of the initial conditions is increased, the deviation from linear Landau damping appears sooner and is more severe. In some cases damping was observed to cease. Curves showing the time dependence of the damping decrement are derived and compared with predictions of non-linear theories. New results obtained in this study include the observation that for strongly non-linear cases, the damping of the electric field causes an initially Maxwellian fo (v, O) to develop a peak in the neighborhood of the phase velocity; strong growth of the second harmonic is seen after fo (v, t) develops such a peak. Also new in this study is the interpretation of the development of a certain class of strongly unstable initial conditions as approaching an inhomogeneous equilibrium.

NUMERICAL INTEGRATION↗
NASA NTRS2018

Analysis of High-Speed Aerodynamics of a Swept Wing with Seamless Flaps

Computational fluid dynamics (CFD) analysis was conducted to evaluate the high-speed aerodynamics of a Gulfstream GIII airplane (Gulfstream Aerospace Corporation, Savannah, Georgia) swept wing modified with an experimental seamless, compliant flap called the Adaptive Compliant Trailing Edge (ACTE) flap for airworthiness. The high-speed aerodynamics of the modified ACTE wing were analyzed at a Mach number of 0.85 and an altitude of 40,000 ft. A polyhedral finite-volume unstructured Reynolds-averaged Navier-Stokes CFD code was used with a two-equation turbulence model. Wing aerodynamics sensitivities to lift coefficients, airplane engine operation, vortex generators (VGs), and ACTE flap deflections were examined. High-speed flow over the clean GIII wing without VGs was found to be complex with outboard shock and boundary-layer interactions causing flow separation over the upper wing surface upstream of the aileron, which could result in buffet and controllability problems at higher lift coefficient values. The VGs were found to improve the high-speed aerodynamics of the clean wing, as expected; however, ACTE flap deflections of as little as 1 to 2 deg were found to negate all of the aerodynamic benefits provided by the VGs at high speeds.

wings↗
NASA NTRS1995

Supersonic solutions of the solar-wind equations

We re-examine the inviscid solar-wind equations with heat conduction from below, and establish a fundamentally new approach for finding stellar solutions. Although the problem is fourth order, only two independent integration constants can be assigned, since two boundary conditions that are required to specify well-behaved supersonic solutions determine the values of the other two constants. The solar-wind-model of Noble and Scarf are essentially accurate for practical purposes, but in a fundamental sense they are not self-consistent. At the supersonic level, the ratio of thermal energy to gravitational potential, (kTr/GMm)(sub s), must lie in a narrow range, between 0.4375 and about 0.40, to permit well-behaved solutions for ionized hydrogen. For a supersonic solution this ratio uniquely determines the constants A and epsilon, which in turn govern not only the escape flux and the energy flux per particle but also T'(lambda approaches 0), the temperature gradient at infinity. It is very difficult to find solutions by trial and error, without employing T'(0) as a guide to the iteration of A and epsilon, because not only must (kTr/GMm)(sub s) be chosen within the acceptable range but then either A or epsilon must be guessed with phenomenal accuracy. Otherwise the numerical integration of temperature and expansion velocity do not yield acceptable values at infinity.

Chamberlain, Joseph W.↗