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.
The classic slotted-wall boundary-condition coefficient for rod-wall wind tunnels is derived by approximating the potential flow solution through a cascade of two staggered rows of rods. A comparison with the corrected Chen and Mears solution for flow through an unstaggered cascade is made.
An outstanding issue in the computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.
Eleven sensitivity experiments that were made with general circulation models to see how land-surface boundary conditions can influence the rainfall, temperature, and motion fields of the atmosphere are discussed. In one group of experiments, different soil moistures or albedos are prescribed as time-invariant boundary conditions. In a second group, different soil moistures or different albedos are initially prescribed, and the soil moisture (but not the albedo) is allowed to change with time according to the governing equations for soil moisture. In a third group, the results of constant versus time-dependent soil moistures are compared.
Eleven sensitivity experiments that were made with general circulation models to see how land-surface boundary conditions can influence the rainfall, temperature, and motion fields of the atmosphere are discussed. In one group of experiments, different soil moistures or albedos are prescribed as time-invariant boundary conditions. In a second group, different soil moistures or different albedos are initially prescribed, and the soil moisture (but not the albedo) is allowed to change with time according to the governing equations for soil moisture. In a third group, the results of constant versus time-dependent soil moistures are compared. Previously announced in STAR as N83-27536
The idea of replacing a divergence constraint by a divergence boundary condition is investigated. The connections between the formulations are considered in detail. It is shown that the most common methods of using divergence boundary conditions do not always work properly. Necessary and sufficient conditions for the equivalence of the formulations are given.
In this work, we derive the interface exchange boundary conditions for the classical linear dynamics of magnetization in ferromagnetic layers with the interface Dzyaloshinskii-Moriya interaction (IDMI). We show that IDMI leads to pinning of dynamic magnetization at the interface. An unusual peculiarity of the IDMI-based pinning is that its scales as the spin-wave wave number. We incorporate these boundary conditions into an existing numerical model for the dynamics of the Damon-Eshbach spin wave in ferromagnetic films. IDMI affects the dispersion and the frequency non-reciprocity of the travelling Damon-Eshbach spin wave. For a broad range of film thicknesses L and wave numbers, the results of the numerical simulations of the spin wave dispersion are in a good agreement with a simple analytical expression, which shows that the contribution of IDMI to the dispersion scales as 1/L, similarly to the effect of other types of interfacial anisotropy. Suggestions to experimentalists how to detect the presence of IDMI in a spin wave experiment are given.
We discuss in this paper the computational efficiency of approximate discrete ordinates (SN) albedo boundary conditions for two-energy group eigenvalue problems in X,Y-geometry. The non-standard SN albedo substitutes approximately the reflector system around the active domain, as we neglect the transverse leakage terms within the non-multiplying reflector region. Should the problem have no transverse leakage terms, i.e., one-dimensional slab geometry, then the offered albedo boundary conditions are exact. By computational efficiency we mean analyzing the accuracy of the numerical results versus the CPU execution time of each run for a given model problem. Numerical results to a typical test problem are shown to illustrate this efficiency analysis. (authors)
An isotope-enabled ocean-atmosphere general circulation model (GISS ModelE -R) is used to estimate the spatial gradients of the oxygen isotopic composition of seawater (delta O-18(sub sw), where delta is the deviation from a known reference material in per mil) during the early Paleogene (45.65 Ma). Understanding the response of delta O-18(sub sw) to changes in climatic and tectonic boundary conditions is important because records of carbonate delta O-18 document changes in hydrology, as well as changes in temperature and global ice -volume. We present results from an early Paleogene configuration of ModelE -R which indicate that spatial gradients of surface ocean delta O-18(sub sw) during this period could have been significantly different to those in the modern ocean. The differences inferred from ModelE -R are sufficient to change early Paleogene sea surface temperature estimates derived from primary carbonate delta O-18 signatures by more than +/-2 C in large areas of the ocean. In the North Atlantic, Indian, and Southern Oceans, the differences in d18Osw inferred from our simulation with ModelE -R are in direct contrast with those from another d18O ]tracing model study which used different, but equally plausible, early Paleogene boundary conditions. The large differences in delta O-18(sub sw) between preindustrial and early Paleogene simulations, and between models, emphasizes the sensitivity of d18Osw to climatic and tectonic boundary conditions. For this reason, absolute estimates of Eocene/ Paleocene temperature derived from carbonate delta O-18 alone are likely to have larger uncertainties than are usually assumed.
Users manual for a U.S.A. FORTRAN 4 computer program which determines boundary conditions for a thin wing lifting surface program is described. This program, the geometry program, and several other programs are used together in the analysis of lifting, thin wings in steady, subsonic flow according to a kernel function lifting surface theory. The program calculates specific types of boundary conditions automatically such as those necessary to determine pitch and roll damping derivatives. The program also accepts descriptions of the camber or downwash and twist in the form of tables and/or coefficients of equations. The program performs interpolations so that tables and/or coefficients can apply at stations selected by the user and not at stations dictated by the control point locations.
Absorption of the electromagnetic energy in a semi-infinite electron plasma is calculated for an arbitrary degree of the electron gas degeneracy. Absorption is determined by solving the boundary-value problem on the oscillations of electron plasma in a half-space with mirror boundary conditions for electrons. The Vlasov−Boltzmann kinetic equation with the Bhatnagar–Gross–Krook collision integral for the electron distribution function and Maxwell’s equation for the electric field are employed. The electron distribution function and the electric field inside plasma are searched for in the form of expansions in the eigenfunctions of the initial set of equations. The expansion coefficients are found for the case of mirror boundary conditions. The contribution of the plasma surface to absorption is analyzed. Cases with different degrees of electron gas degeneracy are considered. It is shown that absorption of the electromagnetic energy near the surface depends substantially on the ratio between the electric field frequency and the volumetric electron collision frequency.
Inconsistent discrete expressions in the boundary treatment of Navier–Stokes solvers and in the definition of force objective functionals can lead to discrete-adjoint boundary treatments that are not a valid representation of the boundary conditions to the corresponding adjoint partial differential equations. The underlying problem is studied for an elementary 1D advection–diffusion problem first using a node-centred finite-volume discretisation. The defect of the boundary operators in the inconsistently defined discrete-adjoint problem leads to oscillations and becomes evident with the additional insight of the continuous-adjoint approach. A homogenisation of the discretisations for the primal boundary treatment and the force objective functional yields second-order functional accuracy and eliminates the defect in the discrete-adjoint boundary treatment. Subsequently, the issue is studied for aerodynamic Reynolds-averaged Navier–Stokes problems in conjunction with a standard finite-volume discretisation on median-dual grids and a strong implementation of noslip walls, found in many unstructured general-purpose flow solvers. Going out from a base-line discretisation of force objective functionals which is independent of the boundary treatment in the flow solver, two improved flux-consistent schemes are presented; based on either body wall-defined or farfield-defined control-volumes they resolve the dual inconsistency. The behaviour of the schemes is investigated on a sequence of grids in 2D and 3D.
For quasi-geostraphic stationary long waves forced by topography, the nonlinear lower boundary condition is derived in terms of the geopotential height and compared with the linearized version. The common practice of replacing terms describing the flow over and around a mountain by upstream zonal flow over the mountain and evaluating the resulting condition at sea level is found to be a good approximation for the cases considered and does not need to be modified as sometimes suggested. Specifically, it is found that this approximation does not affect, for most cases, the lower boundary condition expressed in terms of the geopotential height provided that the stationary wave is not near resonance. At resonance, the eddy advection terms may become important for large-amplitude waves when dissipation and surface diabatic heating are taken into account.
A number of researchers have suggested using an inter-blade row boundary condition to extend isolated blade row flow solvers to multiple blade row configurations. This suggestion is worth consideration for it appears to result in codes that are computationally more efficient than those based on other schemes that were suggested to accomplish the same task. The work is concerned with the development of a mathematical constraint which this boundary condition must satisfy to insure the proper transfer of momentum and vorticity across the plane. Using experimental data, the work quantifies the error in the time-averaged vorticity field which results from simply requiring continuity across the boundary plane of the momentum based on the time-averaged velocity fields associated with a multiple blade row configuration.
This paper will present results in computational nonlinear optics. An algorithm will be described that provides absorbing boundary conditions for optical pulses in dispersive, nonlinear materials. A new numerical absorber at the boundaries has been developed that is responsive to the spectral content of the pulse. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of "light bullet" like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. Comparisons will be shown of calculations that use the standard boundary conditions and the new ones.
For applications of ICRF power in fusion devices, control of RF sheath interactions is of great importance. A sheath boundary condition (SBC) was previously developed to provide an effective surface impedance for the interaction of the RF sheath with the waves. The SBC enables the surface power flux and rectified potential energy available for sputtering to be calculated. For legacy codes which cannot easily implement the SBC, or to speed convergence in codes which do implement it, we consider here an approximate method to simulate SBCs by post-processing results obtained using other, e.g. conducting wall, boundary conditions. The basic approximation is that the modifications resulting from the generalized SBC are driven by a fixed incoming wave which could be either a fast wave or a slow wave. Finally, the method is illustrated in slab geometry and compared with exact numerical solutions; it is shown to work very well.
A numerical scheme suitable for the computation of both the near field acoustic sources and the far field sound produced by turbulent free shear flows utilizing the Navier-Stokes equations is presented. To produce stable numerical schemes in the presence of shear, damping terms must be added to the boundary conditions. The numerical technique and boundary conditions are found to give stable results for computations of spatially evolving mixing layers.
Computational fluid dynamic (CFD) simulations of models tested in wind tunnels require a high level of fidelity and accuracy particularly for the purposes of CFD validation efforts. Considerable effort is required to ensure the proper characterization of both the physical geometry of the wind tunnel and recreating the correct flow conditions inside the wind tunnel. The typical trial-and-error effort used for determining the boundary condition values for a particular tunnel configuration are time and computer resource intensive. This paper describes a method for calculating and updating the back pressure boundary condition in wind tunnel simulations by using a proportional-integral-derivative controller. The controller methodology and equations are discussed, and simulations using the controller to set a tunnel Mach number in the NASA Langley 14- by 22-Foot Subsonic Tunnel are demonstrated.
A method is presented for formulating the boundary conditions in implicit finite-difference form needed for obtaining solutions to the compressible Navier-Stokes equations by the Beam and Warming implicit factored method. The usefulness of the method was demonstrated (a) by establishing the boundary conditions applicable to the analysis of the flow inside an axisymmetric piston-cylinder configuration and (b) by calculating velocities and mass fractions inside the cylinder for different geometries and different operating conditions. Stability, selection of time step and grid sizes, and computer time requirements are discussed in reference to the piston-cylinder problem analyzed.