Satisfaction of asymptotic boundary conditions in the numerical solution of boundary-layer equations.
Satisfaction of asymptotic boundary conditions in numerical solution of boundary layer equations
ASYMPTOTIC METHOD↗
SEARCH · Engineering Papers
Results for “Boundary Conditions”
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.
Time dependent boundary conditions for hyperbolic systems
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.
Thompson, Kevin W.↗
Implementation of a Compressor Face Boundary Condition Based on Small Disturbances
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.
Slater, John W.↗
Surface boundary conditions for the numerical solution of the Euler equations
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.
Dadone, A.↗
A Discrete Analysis of Non-reflecting Boundary Conditions for Discontinuous Galerkin Method
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.
Hu, Fang Q.↗
Supersonic far-field boundary conditions for transonic small-disturbance theory
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.
Gibbons, Michael D.↗
Absorbing boundary conditions for second-order hyperbolic equations
A uniform approach to construct absorbing artificial boundary conditions for second-order linear hyperbolic equations is proposed. The nonlocal boundary condition is given by a pseudodifferential operator that annihilates travelling waves. It is obtained through the dispersion relation of the differential equation by requiring that the initial-boundary value problem admits the wave solutions travelling in one direction only. Local approximation of this global boundary condition yields an nth-order differential operator. It is shown that the best approximations must be in the canonical forms which can be factorized into first-order operators. These boundary conditions are perfectly absorbing for wave packets propagating at certain group velocities. A hierarchy of absorbing boundary conditions is derived for transonic small perturbation equations of unsteady flows. These examples illustrate that the absorbing boundary conditions are easy to derive, and the effectiveness is demonstrated by the numerical experiments.
Jiang, Hong↗
Absorbing boundary conditions for second-order hyperbolic equations
A uniform approach to construct absorbing artificial boundary conditions for second-order linear hyperbolic equations is proposed. The nonlocal boundary condition is given by a pseudodifferential operator that annihilates travelling waves. It is obtained through the dispersion relation of the differential equation by requiring that the initial-boundary value problem admits the wave solutions travelling in one direction only. Local approximation of this global boundary condition yields an nth-order differential operator. It is shown that the best approximations must be in the canonical forms which can be factorized into first-order operators. These boundary conditions are perfectly absorbing for wave packets propagating at certain group velocities. A hierarchy of absorbing boundary conditions is derived for transonic small perturbation equations of unsteady flows. These examples illustrate that the absorbing boundary conditions are easy to derive, and the effectiveness is demonstrated by the numerical experiments.
Jiang, Hong↗
A direct approach to finding unknown boundary conditions in steady heat conduction
The capability of the boundary element method (BEM) in determining thermal boundary conditions on surfaces of a conducting solid where such quantities are unknown was demonstrated. The method uses a non-iterative direct approach in solving what is usually called the inverse heat conduction problem (IHCP). Given any over-specified thermal boundary conditions such as a combination of temperature and heat flux on a surface where such data is readily available, the algorithm computes the temperature field within the object and any unknown thermal boundary conditions on surfaces where thermal boundary values are unavailable. A two-dimensional, steady-state BEM program was developed and was tested on several simple geometries where the analytic solution was known. Results obtained with the BEM were in excellent agreement with the analytic values. The algorithm is highly flexible in treating complex geometries, mixed thermal boundary conditions, and temperature-dependent material properties and is presently being extended to three-dimensional and unsteady heat conduction problems. The accuracy and reliability of this technique was very good but tended to deteriorate when the known surface conditions were only slightly over-specified and far from the inaccessible surface.
Martin, Thomas J.↗
Boundary Condition for Modeling Semiconductor Nanostructures
A recently proposed boundary condition for atomistic computational modeling of semiconductor nanostructures (particularly, quantum dots) is an improved alternative to two prior such boundary conditions. As explained, this boundary condition helps to reduce the amount of computation while maintaining accuracy.
Lee, Seungwon↗
Characteristic boundary conditions for the Euler equations
The boundary conditions are demonstrated for the quasi-one-dimensional Euler equations with the extension to two and three dimensions being straightforward. In this application an implicit finite-difference scheme is employed with the boundary conditions being applied implicitly. The boundary application uses both characteristic extrapolations and evaluations which distinguishes it from other theories. Flow fields with shocks are calculated with inflow-outflow conditions of supersonic-subsonic and subsonic-subsonic flow.
T H Pulliam↗
Numerical Study of Outlet Boundary Conditions for Unsteady Turbulent Internal Flows Using the NCC
This paper presents the results of studies on the outlet boundary conditions for turbulent internal flow simulations. Several outlet boundary conditions have been investigated by applying the National Combustion Code (NCC) to the configuration of a LM6000 single injector flame tube. First of all, very large eddy simulations (VLES) have been performed using the partially resolved numerical simulation (PRNS) approach, in which both the nonlinear and linear dynamic subscale models were employed. Secondly, unsteady Reynolds averaged Navier- Stokes (URANS) simulations have also been performed for the same configuration to investigate the effects of different outlet boundary conditions in the context of URANS. Thirdly, the possible role of the initial condition is inspected by using three different initial flow fields for both the PRNS/VLES simulation and the URANS simulation. The same grid is used for all the simulations and the number of mesh element is about 0.5 million. The main purpose of this study is to examine the long-time behavior of the solution as determined by the imposed outlet boundary conditions. For a particular simulation to be considered as successful under the given initial and boundary conditions, the solution must be sustainable in a physically meaningful manner over a sufficiently long period of time. The commonly used outlet boundary condition for steady Reynolds averaged Navier-Stokes (RANS) simulation is a fixed pressure at the outlet with all the other dependent variables being extrapolated from the interior. The results of the present study suggest that this is also workable for the URANS simulation of the LM6000 injector flame tube. However, it does not work for the PRNS/VLES simulation due to the unphysical reflections of the pressure disturbances at the outlet boundary. This undesirable situation can be practically alleviated by applying a simple unsteady convection equation for the pressure disturbances at the outlet boundary. The numerical results presented in this paper suggest that this unsteady convection of pressure disturbances at the outlet works very well for all the unsteady simulations (both PRNS/VLES and URANS) of the LM6000 single injector flame tube.
Liu, Nan-Suey↗
Effect of Inflow Boundary Conditions on Hovering Tilt-Rotor Flows
A study on the effect of inflow turbulence boundary conditions on the local flow on and around a body in flight has been carried out. The study has been carried out using OVERFLOW2 flow solver using the default Spalart-Allmaras turbulence model in OVERFLOW. Many OVERFLOW turbulent flow simulations have been reported using the SA-fv3 model1 over the years. The present study demonstrates that the turbulence levels imposed as a boundary condition (b.c.) in the far-field as implemented in OVERFLOW are not correct. In fact, very low level of turbulence at the far-field boundary as implemented in OVERFLOW results in the SA-fv3 model predicting transition-like profiles on a given body. By choosing sufficiently high levels of the Reynolds number of turbulence, Rt, as a boundary condition, this anomalous behavior of SA-fv3 model is eliminated. Since numerous papers using OVERFLOW have been presented in the literature including that by the author2 using low level of inflow turbulence, it will be beneficial to the CFD community at large and in particular to the OVERFLOW community to understand this effect of the inflow b.c. in the SA-fv3 model. Various results reported over the years in domains such as high-lift applications, drag prediction applications and rotorcraft flow applications can be revisited using the new turbulence inflow boundary condition prescription as suggested in this study. It has been demonstrated in the literature that the effect of inflow turbulence level on the downstream development of flow is too important to be ignored. Demonstrative results are shown and compared with experiment for the case of hovering XV-15 rotor flow at a tip Reynolds number of 4.9 x 10(exp 6) and Mach number of 0.69 in Figs. 1 and 2 below. Fig. 1(a,b) shows skin friction predictions at radial station of r/R = 0.28 corresponding to 5 different inflow boundary conditions obtained with the SA-fv3 model. As the inflow turbulence is increased monotonically, from a value of Rt = 0.1 to 100.0 as shown in Fig. 1, a fully developed turbulent skin friction profile is realized progressively as shown in Fig. 1(b). Similarly skin friction results shown in Fig. 2(a,b) corresponding to radial stations of r/R = 0.72 reinforce the postulate made here that as the inflow turbulence level is increased gradually, fully turbulent on-the-blade simulation is realized uniformly. In the full paper extensive results on skin friction, pressure and velocity profiles as well as the Figure of Merit will be presented and compared with experiment that will illustrate the substantial effect of the inflow boundary conditions on the development of hovering rotor flow downstream.
Kaul, Upender K.↗
Testing horizontal momentum boundary conditions in MPAS-Ocean
A rising concern within E3SM’s MPAS-Ocean is the implementation of boundary conditions. Currently, MPAS-Ocean utilizes a no-slip boundary condition, meaning that the velocity of fluid at the boundaries is zero; however, the model should should satisfy both no-slip and free-slip boundary conditions given the correct parameters. This study focuses on testing these boundary conditions in depth, determining the accuracy of the current implementation. To accomplish this, we focus on two test cases: the barotropic gyre and barotropic channel. By manipulating parameters and comparing the numerical to corresponding analytical solutions (when available), we show that MPAS-Ocean’s ability to satisfy both boundary conditions is promising.
54 ENVIRONMENTAL SCIENCES↗
Effects of Horizontal Resolution, Domain Size, Boundary Conditions, and Surface Heterogeneity on Coarse LES of a Convective Boundary Layer
Atmospheric properties in a convective boundary layer vary over a wide range of spatial scales and are commonly studied using large-eddy simulations (LES) in various configurations. We examine how the boundary layer depth and distribution of variability across scales are affected by LES grid spacing, domain size, inhomogeneity of surface properties, and external forcing. Two different setups of the Weather Research and Forecasting (WRF) model are analyzed. A semi-idealized configuration uses a periodic domain, flat surface, prescribed homogeneous surface heat fluxes, and horizontally uniform profiles of large-scale advective tendencies. A nested LES setup employs a larger domain and realistic initial and boundary conditions, including an interactive land surface model with representative topography and vegetation and soil types. Subdomains of identical size are analyzed for all simulations. Characteristic structure sizes are quantified using the variability scales L 50 and L 95 , defined such that features smaller than that contain 50% and 95% of the total variance, respectively. Progressive increase in L 50 from vertical velocity to temperature and moisture structures is systematically reproduced in all simulation configurations. This dependence of L 50 on the considered variable complicates the development of scale-aware parameterizations for models with grid spacing in the “terra incognita”. In simulations using a larger domain with heterogeneous surface properties, the development of internal mesoscale patterns significantly affects variance distributions inside analyzed subdomains. Sizes of boundary layer structures also strongly depend on the LES grid spacing and, in case of heterogeneous surface and topography, on location of the subdomain inside a larger computational domain.
54 ENVIRONMENTAL SCIENCES↗
Long-time behaviour of absorbing boundary conditions
A new class of computational far-field boundary conditions for hyperbolic partial differential equations was recently introduced by the authors. These boundary conditions combine properties of absorbing conditions for transient solutions and properties of far-field conditions for steady states. This paper analyses the properties of the wave equation coupled with these new boundary conditions: well-posedness, dissipativity and convergence in time.
Engquist, B.↗
On the accurate long-time solution of the wave equation in exterior domains: Asymptotic expansions and corrected boundary conditions
We consider the solution of scattering problems for the wave equation using approximate boundary conditions at artificial boundaries. These conditions are explicitly viewed as approximations to an exact boundary condition satisfied by the solution on the unbounded domain. We study the short and long term behavior of the error. It is provided that, in two space dimensions, no local in time, constant coefficient boundary operator can lead to accurate results uniformly in time for the class of problems we consider. A variable coefficient operator is developed which attains better accuracy (uniformly in time) than is possible with constant coefficient approximations. The theory is illustrated by numerical examples. We also analyze the proposed boundary conditions using energy methods, leading to asymptotically correct error bounds.
Hagstrom, Thomas↗
Satisfaction of asymptotic boundary conditions in the numerical solution of boundary-layer equations.
Boundary layer differential equations numerical solution with asymptotic boundary conditions satisfaction, using least squares convergence criterion
Nachtsheim, P. R.↗