SEARCH · Engineering Papers
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.
Boundary condition iteration methods in chemical process plant optimal design and control, considering simplicity and ease of programming, computation time and storage requirements
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.
The influence of the exit boundary conditions (vanishing first derivative of the velocity components and constant pressure) on the large eddy simulation of the fully developed turbulent channel flow has been investigated for equidistant and stretched grids at the channel exit. Results show that the chosen exit boundary conditions introduce some small disturbance which is mostly damped by the grid stretching. The difference between the fully developed turbulent channel flow obtained with LES with periodicity condition and the inlet and exit and the LES with fully developed flow at the inlet and the exit boundary condition is less than 10% for equidistant grids and less than 5% for the case grid stretching. The chosen boundary condition is of interest because it may be used in complex flows with backflow at exit.
A time domain boundary integral equation with Burton-Miller reformulation is presented for acoustic scattering by surfaces with liners in a uniform mean flow. The Ingard-Myers impedance boundary condition is implemented using a broadband multipole impedance model and converted into time domain differential equations to augment the boundary integral equation. The coupled integral-differential equations are solved numerically by a March-On-in-Time (MOT) scheme. While the Ingard-Myers condition is known to support Kelvin-Helmholtz instability due to its use of a vortex sheet interface between the flow and the liner surface, it is found that by neglecting a second derivative term in the current time domain impedance boundary condition formulation, the instability an be effectively suppressed in computation. The proposed formulation and implementation are validated using a NASA Langley Research Center Grazing Flow Impedance Tube (GFIT) experimental dataset with satisfactory results. Moreover, a minimization procedure for finding the poles and coefficients of the broadband multiple impedance model is formulated in this paper by which, unlike the commonly used vector-fitting method, passivity of the model is ensured. Numerical tests show the proposed minimization approach is effective for modeling liners that are commonly used in aeroacoustic applications.
Cell-centered finite-volume (CCFV) schemes have certain attractive properties for the solution of the equations governing compressible fluid flow. Among others, they provide a natural vehicle for specifying flux conditions at the boundaries of the physical domain. Unfortunately, they lead to slow convergence for numerical programs utilizing them. In this report a method for investigating and improving the convergence of CCFV schemes is presented, which focuses on the effect of the numerical boundary conditions. The key to the method is the computation of the spectral radius of the iteration matrix of the entire demoralized system of equations, not just of the interior point scheme or the boundary conditions.
A nonreflecting boundary condition is presented for the numerical solution of the time-dependent compressible Navier-Stokes equations when these equations are used to obtain a steady state. This boundary condition is shown to be effective in reducing reflections at a subsonic outflow boundary. Numerical calculations using a model problem were made to compare this boundary condition with other outflow boundary conditions. The nonreflecting boundary condition contains a parameter whose optimal value is estimated using the analysis of a simplified set of equations.
A general boundary condition formalism is developed for all types of boundary conditions to which hyperbolic systems are subject; the formalism makes possible a 'cookbook' approach to boundary conditions, by means of which novel boundary 'recipes' may be derived and previously devised ones may be consulted as required. Numerous useful conditions are derived for such CFD problems as subsonic and supersonic inflows and outflows, nonreflecting boundaries, force-free boundaries, constant pressure boundaries, and constant mass flux. Attention is given to the computation and integration of time derivatives.
A boundary condition for the linear acoustic equation has been developed that allows the acoustic pressure waves to propagate out of the computational domain boundary, just as they would propagate in an infinitely long duct. The problem is divided into two domains: numerical and analytical. The boundary condition provides a matching of the two domains. Examples show this method works well in an acoustic panel program for a model problem (simple source in a rectangular duct with several propagating modes present). The boundary condition is described so that it can be used with various duct various duct geometries and numerical methods.
Implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solid objects subjected to nonlinear boundary conditions is shown to generate an accurate and economical approach for the computation of shape sensitivities for this class of problems. This approach involves the employment of analytical derivatives of boundary-element kernel functions with respect to shape design variables. A formulation is presented that can consistently account for both temperature-dependent convection and radiation boundary conditions. Several iterative strategies are presented for the solution of the resulting sets of nonlinear equations and the computational performances examined in detail. Multizone analysis and zone condensation strategies are demonstrated to provide substantive computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear example problems are presented that have closed-form solutions.
Surface impedance boundary conditions are used to reduce the solution volume during the analysis of scattering from lossy dielectric objects. In a finite difference solution, they also can be used to avoid using small cells, made necessary by shorter wavelengths in conducting media, throughout the solution volume. A one dimensional implementation is presented for a surface impedance boundary condition for good conductors in the Finite Difference Time Domain (FDTD) technique. In order to illustrate the FDTD surface impedance boundary condition, a planar air-lossy dielectric interface is considered.
Accurate numerical non-reflecting boundary conditions are important in all the proposed benchmark problems of the Second Workshop. Recently, a new absorbing boundary condition has been developed using Perfectly Matched Layer (PML) equations for the Euler equations. In this approach, a region with a width of a few grid points is introduced adjacent to the non-reflecting boundaries. In the added region, Perfectly Matched Layer equations are constructed and applied so that the out-going waves are absorbed inside the layer with little reflection to the interior domain. It will be demonstrated in the present paper that the proposed absorbing boundary condition is quite general and versatile, applicable to radiation boundaries as well as inflow and outflow boundaries. It is also easy to implement. The emphasis of the paper will be on the application of the PML absorbing boundary condition to problems in Categories 1, 2, and 3. In Category 1, solutions of problems 1 and 2 are presented. Both problems are solved using a multi-domain polar grid system. Perfectly Matched Layer equations for a circular boundary are constructed and their effectiveness assessed. In Category 2, solutions of problem 2 are presented. Here, in addition to the radiation boundary conditions at the far field in the axisymmetric coordinate system, the inflow boundary condition at the duct inlet is also dealt with using the proposed Perfectly Match Layer equations. At the inlet, a PML domain is introduced in which the incident duct mode is simulated while the waves reflected from the open end of the duct are absorbed at the same time. In Category 3, solutions of all three problems are presented. Again, the PML absorbing boundary condition is used at the inflow boundary so that the incoming vorticity wave is simulated while the outgoing acoustic waves are absorbed with very little numerical reflection. All the problems are solved using central difference schemes for spatial discretizations and the optimized Low-Dissipation and Low-Dispersion Runge-Kutta scheme for the time integration. Issues of numerical accuracy and efficiency are also addressed.
A new class of computational far field boundary conditions for hyperbolic partial differential equations is developed. These boundary conditions combine properties of absorbing boundary conditions for transient solutions and properties of far field boundary conditions for steady-state problems. The conditions can be used to limit the computational domain when both traveling waves and evanescent waves are present. Boundary conditions for scalar wave equations are derived and analyzed. Extensions to systems of equations are discussed and results from numerical experiments are presented.
Satisfaction of asymptotic boundary conditions in numerical solution of boundary layer equations
Nonreflecting boundary conditions are defined for multidimensional fluid dynamics problems where waves enter and leave the interior of a domain modeled by hyperbolic equations. Separate equations are defined for each type of incoming and outgoing wave. Temporally varying problems are considered in terms of a nonreflecting boundary condition which permit the amplitude of incoming waves to remain constant over time. Conservative expressions are presented that include dissipative terms. Applications of the computational techniques are illustrated with sample results for a traveling shock wave, a shock tube, a spherical explosion and expansion problems on one- and two-dimensions.
A compressor-face boundary condition that models the unsteady interactions of acoustic and convective velocity disturbances with a compressor has been implemented into a three-dimensional computational fluid dynamics code. Locally one-dimensional characteristics along with a small-disturbance model are used to compute the acoustic response as a function of the local stagger angle and the strength and direction of the disturbance. Simulations of the inviscid flow in a straight duct, a duct coupled to a compressor, and a supersonic inlet demonstrate the behavior of the boundary condition in relation to existing boundary conditions. Comparisons with experimental data show a large improvement in accuracy over existing boundary conditions in the ability to predict the reflected disturbance from the interaction of an acoustic disturbance with a compressor.
We consider the implementation of boundary conditions at solid walls in inviscid Euler solutions by upwind, finite-volume methods. We review some current methods for the implementation of surface boundary conditions and examine their behavior for the problem of an oblique shock reflecting off a planar surface. We show the importance of characteristic boundary conditions for this problem and introduce a method of applying the classical flux-difference splitting of Roe as a characteristic boundary condition. Consideration of the equivalent problem of the intersection of two (equal and opposite) oblique shocks was very illuminating on the role of surface boundary conditions for an inviscid flow and led to the introduction of two new boundary-condition procedures, denoted as the symmetry technique and the curvature-corrected symmetry technique. Examples of the effects of the various surface boundary conditions considered are presented for the supersonic blunt body problem and the subcritical compressible flow over a circular cylinder. Dramatic advantages of the curvature-corrected symmetry technique over the other methods are shown, with regard to numerical entropy generation, total pressure loss, drag and grid convergence.
We present a discrete analysis of non-reflecting boundary conditions for the discontinuous Galerkin method. The boundary conditions considered in this paper include the recently proposed Perfectly Matched Layer absorbing boundary condition for the linearized Euler equation and two non-reflecting boundary conditions based on the characteristic decomposition of the flux on the boundary. The analyses for the three boundary conditions are carried out in a unifled way. In each case, eigensolutions of the discrete system are obtained and applied to compute the numerical reflection coefficients of a specified out-going wave. The dependencies of the reflections at the boundary on the out-going wave angle and frequency as well as the mesh sizes arc? studied. Comparisons with direct numerical simulation results are also presented.
Characteristic far-field boundary conditions for supersonic freestream flow have been developed and implemented within a transonic small-disturbance code. The boundary conditions have been implemented within the CAP-TSD code which has been developed recently for aeroelastic analysis of complete aircraft configurations. These boundary conditions improve the accuracy of the solutions for supersonic freestream applications. They also allow the extent of the grid to be much smaller, thus providing savings in the computational time required to obtain solutions. Comparisons are shown between surface pressures computed using large and small grid extents for the NACA 0012 airfoil and the F-5 wing at various Mach numbers and angles of attack. Both steady and unsteady results are presented and comparisons are made with Euler results and with experimental data to assess the accuracy of the new far-field boundary conditions. Comparisons of these results show that the supersonic boundary conditions allow a much smaller grid to be used without losing accuracy.