Search NASA⌕ Search

SEARCH · Search NASA

Results for “Integral boundary condition”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Integral boundary conditions in phase field models

Modeling the chemical, electric and thermal transport as well as phase transitions and the accompanying mesoscale microstructure evolution within a material in an electronic device setting involves the solution of partial differential equations often with integral boundary conditions. Employing the familiar Poisson equation describing the electric potential evolution in a material exhibiting insulator to metal transitions, we exploit a special property of such an integral boundary condition, and we properly formulate the variational problem and establish its well-posedness. Next, we compare our method with the commonly-used Lagrange multiplier method that can also handle such boundary conditions. Numerical experiments demonstrate that our new method achieves optimal convergence rate in contrast to the conventional Lagrange multiplier method. Furthermore, the linear system derived from our method is symmetric positive definite, and can be efficiently solved by Conjugate Gradient method with algebraic multigrid preconditioning.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Coupling conditions for integrating boundary layer and rotational inviscid flow

The matching of a boundary layer and a rotational inviscid flow is reexamined by extending the Johnson and Sockol (1979) coupling conditions to include the case where the boundary layer solution includes the second-order effects of the freestream vorticity and the total temperature gradient. It is pointed out that two of the three conditions are not independent. If the boundary layer solution satisfies the appropriate momentum and energy integral equations, it follows that the imposition of the normal mass flux condition insures that the conditions on a normal flux of streamwise momentum and total enthalpy will also be satisfied.

Sockol, P. M.↗
NASA NTRS1949

Two-Dimensional Potential Flows

Contents include the following: Characteristic differential equations - initial and boundary conditions. Integration of the second characteristic differential equations. Direct application of Meyer's characteristic hodograph table for construction of two-dimensional potential flows. Prandtl-Busemann method. Development of the pressure variation for small deflection angles. Numerical table: relation between deflection, pressure, velocity, mach number and mach angle for isentropic changes of state according to Prandtl-Meyer for air (k = 1.405). References.

Schaefer, Manfred↗
DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗
NASA NTRS2020

INVESTIGATING THE FEASIBILITY AND STABILITY FOR MODELING ACOUSTIC WAVE SCATTERING USING A TIME-DOMAIN BOUNDARY INTEGRAL EQUATION WITH IMPEDANCE BOUNDARY CONDITION

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet and environmentally friendly aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction and are characterized by a frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. This work considers using either an impedance or an admittance (inverse of impedance) boundary condition to allow for acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces. Three acoustic liner models are discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and the Broadband Impedance Model. In both the Helmholtz and Three-Parameter models, liner impedance is specified at a given frequency, whereas the Broadband model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for acoustic liners are derived for each model and coupled with a time-domain boundary integral equation. The scattering solution is obtained iteratively using a boundary element method with constant spatial and third-order temporal basis functions. Time-domain boundary integral equations are unfortunately prone to numerical instabilities due to resonant frequencies resulting from non-trivial solutions in the interior domain. When reformulated with the Burton-Miller method, the instabilities are eliminated. Using a Burton-Miller reformulation, the stability of the boundary element method assuming a liner boundary condition is assessed using eigenvalue analysis. The stability of each liner model is discussed, and it is shown that the Three-Parameter and Broadband models are sufficient for modeling an acoustic liner on the surface of scattering bodies. The Helmholtz model demonstrates strict limitations for stability, whereas the Three-Parameter and Broadband models are stable for most cases. Also included in this work is an assessment of the spatial accuracy of the time-domain boundary element method with respect to the surface element basis functions, as well as a performance study of the numerical algorithm.

computational mathematics↗
NASA NTRS1988

Investigation of viscous/inviscid interaction in transonic flow over airfoils with suction

The viscous/inviscid interaction over transonic airfoils with and without suction is studied. The streamline angle at the edge of the boundary layer is used to couple the viscous and inviscid flows. The potential flow equations are solved for the inviscid flow field. In the shock region, the Euler equations are solved using the method of integral relations. For this, the potential flow solution is used as the initial and boundary conditions. An integral method is used to solve the laminar boundary-layer equations. Since both methods are integral methods, a continuous interaction is allowed between the outer inviscid flow region and the inner viscous flow region. To avoid the Goldstein singularity near the separation point the laminar boundary-layer equations are derived in an inverse form to obtain solution for the flows with small separations. The displacement thickness distribution is specified instead of the usual pressure distribution to solve the boundry-layer equations. The Euler equations are solved for the inviscid flow using the finite volume technique and the coupling is achieved by a surface transpiration model. A method is developed to apply a minimum amount of suction that is required to have an attached flow on the airfoil. The turbulent boundary layer equations are derived using the bi-logarithmic wall law for mass transfer. The results are found to be in good agreement with available experimental data and with the results of other computational methods.

Vemuru, C. S.↗
NASA NTRS2005

Advanced Small Perturbation Potential Flow Theory for Unsteady Aerodynamic and Aeroelastic Analyses

An advanced small perturbation (ASP) potential flow theory has been developed to improve upon the classical transonic small perturbation (TSP) theories that have been used in various computer codes. These computer codes are typically used for unsteady aerodynamic and aeroelastic analyses in the nonlinear transonic flight regime. The codes exploit the simplicity of stationary Cartesian meshes with the movement or deformation of the configuration under consideration incorporated into the solution algorithm through a planar surface boundary condition. The new ASP theory was developed methodically by first determining the essential elements required to produce full-potential-like solutions with a small perturbation approach on the requisite Cartesian grid. This level of accuracy required a higher-order streamwise mass flux and a mass conserving surface boundary condition. The ASP theory was further developed by determining the essential elements required to produce results that agreed well with Euler solutions. This level of accuracy required mass conserving entropy and vorticity effects, and second-order terms in the trailing wake boundary condition. Finally, an integral boundary layer procedure, applicable to both attached and shock-induced separated flows, was incorporated for viscous effects. The resulting ASP potential flow theory, including entropy, vorticity, and viscous effects, is shown to be mathematically more appropriate and computationally more accurate than the classical TSP theories. The formulaic details of the ASP theory are described fully and the improvements are demonstrated through careful comparisons with accepted alternative results and experimental data. The new theory has been used as the basis for a new computer code called ASP3D (Advanced Small Perturbation - 3D), which also is briefly described with representative results.

Batina, John T.↗
NASA NTRS1977

Downstream boundary conditions for viscous flow problems

The problem of the specification of artificial outflow conditions in flow problems is studied. It is shown that for transport type equations incorrect outflow conditions will adversely affect the solution only in a small region near the outflow boundary, while for elliptic equations, e.g. those governing the streamfunction or pressure, a correct boundary specification is essential. In addition, integral outflow boundary conditions for fluid dynamical problems are considered. It is shown that such conditions are well posed, and their effect on the solutions of the Navier-Stokes equations is also considered.

Fix, G.↗
NASA NTRS1992

Electromagnetic scattering and radiation from microstrip patch antennas and spirals residing in a cavity

A new hybrid method is presented for the analysis of the scattering and radiation by conformal antennas and arrays comprised of circular or rectangular elements. In addition, calculations for cavity-backed spiral antennas are given. The method employs a finite element formulation within the cavity and the boundary integral (exact boundary condition) for terminating the mesh. By virtue of the finite element discretization, the method has no restrictions on the geometry and composition of the cavity or its termination. Furthermore, because of the convolutional nature of the boundary integral and the inherent sparseness of the finite element matrix, the storage requirement is kept very low at O(n). These unique features of the method have already been exploited in other scattering applications and have permitted the analysis of large-size structures with remarkable efficiency. In this report, we describe the method's formulation and implementation for circular and rectangular patch antennas in different superstrate and substrate configurations which may also include the presence of lumped loads and resistive sheets/cards. Also, various modelling approaches are investigated and implemented for characterizing a variety of feed structures to permit the computation of the input impedance and radiation pattern. Many computational examples for rectangular and circular patch configurations are presented which demonstrate the method's versatility, modeling capability and accuracy.

Volakis, J. L.↗
NASA NTRS1973

Introduction to gasdynamics of explosions

Questions concerning the genesis and sustenance of an explosion are investigated, giving attention to the mechanics of explosions, the gasdynamics of explosions, aspects of technological significance, and future prospects. The dynamics of exothermic processes is discussed together with the most prominent effects of explosions. Blast waves are considered, taking into account conservation principles, blast wave transformation, conservative equations in nondimensional form, the equation of state, Eulerian space profiles, Eulerian time profiles, Lagrangian time profiles, boundary conditions and integral relations, and self-similar flow fields.

Oppenheim, A. K.↗
NASA NTRS1976

Turbulent atmospheric flow over a backward-facing step

The phenomenon of atmospheric shear layer separation over a man-made structure such as a building (modeled as a backward-facing step) has been analyzed theoretically by (1) solving the two-dimensional equations of motion in the two variables, stream function and vorticity, and by (2) employing an approximate integral technique. Boundary conditions for the undisturbed flow are that of the turbulent atmospheric shear flow over a rough terrain. In the first approach a two-equation model of turbulence was used. In the second approach an approximate technique was utilized in an attempt to describe the details of the flow in the recirculation zone behind the step. The results predict velocity profiles in sufficient detail that the presence of the corner eddy in the region of negative surface pressure gradient is evident. The magnitude of the reversed flow velocity in the recirculation eddy has been found to agree with that found from experiments. Also, a surface eddy viscosity distribution has been an outgrowth of the method which realistically follows the magnitude of the surface pressure gradient distribution as found experimentally.

Kaul, U. K.↗
NASA NTRS1982

Vibration in Planetary Gear Systems with Unequal Planet Stiffnesses

An algorithm suitable for a minicomputer was developed for finding the natural frequencies and mode shapes of a planetary gear system which has unequal stiffnesses between the Sun/planet and planet/ring gear meshes. Mode shapes are represented in the form of graphical computer output that illustrates the lateral and rotational motion of the three coaxial gears and the planet gears. This procedure permits the analysis of gear trains utilizing nonuniform mesh conditions and user specified masses, stiffnesses, and boundary conditions. Numerical integration of the equations of motion for planetary gear systems indicates that this algorithm offers an efficient means of predicting operating speeds which may result in high dynamic tooth loads.

Frater, J. L.↗
NASA NTRS1985

The Control of Float Zone Interfaces by the Use of Selected Boundary Conditions

Since the float zone process involves two solid-melt interfaces, possible gas interfaces, heat and mass transfers, various driving forces and complex heating sources, an analysis of the entire process would be very complex. For an initial investigation, a more feasible approach is to examine each component of the process separately. The three principal components are: (1) the shapes of the melt and solid-melt interfaces, (2) the heat and mass transfers, and (3) the heating and cooling sources. This study combined facets of all three components. The purpose of this effort was to study and compute the surface boundary conditions required to give flat (float zone) FZ solid-melt interfaces. This study was undertaken in two phases. The first phase was to investigate the solid zones surface boundary conditions required for flat solid-melt interfaces when given the melt zone surface boundary conditions. The second phase complemented the first and was to investigate the melt zone surface boundary conditions required for flat solid-melt interfaces if given the solid zones surface boundary conditions. Dual integral transform methods were used in both phases; in addition, the use of various numerical methods for differential equations and linear systems of equations were required.

Foster, L. M.↗
NASA NTRS1986

Transfer function analysis of thermospheric perturbations

Applying perturbation theory, a spectral model in terms of vectors spherical harmonics (Legendre polynomials) is used to describe the short term thermospheric perturbations originating in the auroral regions. The source may be Joule heating, particle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is considered, allowing for the collisional momentum exchange between species including Ar, O2, N2, O, He and H. The coupled equations of energy, mass and momentum conservation are solved simultaneously for the major species N2 and O. Applying homogeneous boundary conditions, the integration is carred out from the Earth's surface up to 700 km. In the analysis, the spherical harmonics are treated as eigenfunctions, assuming that the Earth's rotation (and prevailing circulation) do not significantly affect perturbations with periods which are typically much less than one day. Under these simplifying assumptions, and given a particular source distribution in the vertical, a two dimensional transfer function is constructed to describe the three dimensional response of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this transfer function reveals five components. To compile the transfer function, the numerical computations are very time consuming (about 100 hours on a VAX for one particular vertical source distribution). However, given the transfer function, the atmospheric response in space and time (using Fourier integral representation) can be constructed with a few seconds of a central processing unit. This model is applied in a case study of wind and temperature measurements on the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in the polar region and then are reflected to propagate toward lower latitudes.

Mayr, H. G.↗
NASA NTRS1989

Theoretical investigation of EM wave generation and radiation in the ULF, ELF, and VLF bands by the electrodynamic orbiting tether

The problem of electromagnetic wave generation by an electrodynamic tethered satellite system is important both for the ordinary operation of such systems and for their possible application as orbiting transmitters. The tether's ionospheric circuit closure problem is closely linked with the propagation of charge-carrying electromagnetic wave packets away from the tethered system. Work is reported which represents a step towards a solution to the problem that takes into account the effects of boundaries and of vertical variations in plasma density, collision frequencies, and ion species. The theory of Alfen wave packet generation by an electrodynamic tethered system in an infinite plasma medium is reviewed, and brief summary of previous work on the problem is given. The consequences of the presence of the boundaries and the vertical nonuniformity are then examined. One of the most significant new features to emerge when ion-neutral collisions are taken into account is the coupling of the Alfven waves to the fast magnetosonic wave. This latter wave is important, as it may be confined by vertical variations in the Alfven speed to a sort of leaky ionospheric wave guide, the resonances of which could be of great importance to the signal received on the Earth's surface. The infinite medium solution for this case where the (uniform) geomagnetic field makes an arbitrary angle with the vertical is taken as the incident wave-packet. Even without a full solution, a number of conclusions can be drawn, the most important of which may be that the electromagnetic field associated with the operation of a steady-current tethered system will probably be too weak to detect on the Earth's surface, even for large tethered currents. This is due to the total reflection of the incident wave at the atmospheric boundary and the inability of a steady-current tethered system to excite the ionospheric wave-guide. An outline of the approach to the numerical problem is given. The use of numerical integrations and boundary conditions consistent with a conducting Earth is proposed to obtain the solution for the horizontal electromagnetic field components at the boundary of the ionosphere with the atmospheric cavity.

Estes, Robert D.↗
NASA NTRS1990

Two-dimensional Euler and Navier-Stokes Time accurate simulations of fan rotor flows

Two numerical methods are presented which describe the unsteady flow field in the blade-to-blade plane of an axial fan rotor. These methods solve the compressible, time-dependent, Euler and the compressible, turbulent, time-dependent, Navier-Stokes conservation equations for mass, momentum, and energy. The Navier-Stokes equations are written in Favre-averaged form and are closed with an approximate two-equation turbulence model with low Reynolds number and compressibility effects included. The unsteady aerodynamic component is obtained by superposing inflow or outflow unsteadiness to the steady conditions through time-dependent boundary conditions. The integration in space is performed by using a finite volume scheme, and the integration in time is performed by using k-stage Runge-Kutta schemes, k = 2,5. The numerical integration algorithm allows the reduction of the computational cost of an unsteady simulation involving high frequency disturbances in both CPU time and memory requirements. Less than 200 sec of CPU time are required to advance the Euler equations in a computational grid made up of about 2000 grid during 10,000 time steps on a CRAY Y-MP computer, with a required memory of less than 0.3 megawords.

Boretti, A. A.↗