Engineering Papers⌕ Search

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.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Numerical tests of a modified full implicit continuous Eulerian (FICE) scheme with projected normal characteristic boundary conditions for MHD flows

A numerical method has been developed based on a modified full implicit continuous Eulerian (FICE) scheme and projected normal characteristic boundary conditions for simulating MHD flows which undergo a long process of evolution. An astrophysical flow is chosen for illustration of this procedure, and numerical tests are made to verify the computational stability and physically realistic solution. Three computational tests have been accomplished; they are tests of solving methods, characteristic boundary condition, and time steps. The tests show that the program from the modified FICE scheme with proper boundary conditions and time steps can be made numerically stable for a time long enough to obtain physically plausible solutions.

Wu, S. T.↗
NASA NTRS1988

A new method of imposing boundary conditions in pseudospectral approximations of hyperbolic equations

A new method to impose boundary conditions for pseudospectral approximations to hyperbolic equations is suggested. This method involves the collocation of the equation at the boundary nodes as well as satisfying boundary conditions. Stability and convergence results are proven for the Chebyshev approximation of linear scalar hyperbolic equations. The eigenvalues of this method applied to parabolic equations are shown to be real and negative.

Funaro, D.↗
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↗
DOE OSTI2019

Absorbing boundary condition (ABC) and perfectly matched layer (PML) in numerical beam propagation: a comparison

The absorbing boundary condition (ABC) and the perfectly matched layer (PML) are employed to suppress the numerical reflections in wave propagation methods. The PML has earlier been shown to fail at the grazing angle propagation of beam. The ABC is preferred over the PML where refractive index is not an analytic function. We have explored the performance of the ABC in reflection-less absorption of non-paraxial propagation of beams and have compared it against the PML. We have considered an example of the propagation of Gaussian beams through a homogeneous medium terminated with ABC or PML. We have used the collocation method for propagation. We have found that both the boundary conditions are comparable in performance. The ABC performed relatively better than the PML for beams at smaller angles. However, the performance of the PML is better than that of the ABC at larger tilt angles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1980

Nonreflecting far-field boundary conditions for unsteady transonic flow computation

The approximate nonreflecting far-field boundary condition, as proposed by Engquist and Majda, is implemented in the computer code LTRAN2. This code solves the implicit finite-difference representation of the small disturbance equations for unsteady transonic flows about airfoils. The nonreflecting boundary condition and the description of the algorithm for implementing these conditions in LTRAN2 are discussed. Various cases are computed and compared with results from the older, more conventional procedures. One concludes that the nonreflecting far-field boundary approximation allows the far-field boundary to be located closer to the airfoil; this permits a decrease in the computer time required to obtain the solution through the use of fewer mesh points.

Kwak, D.↗
NASA NTRS1981

Nonreflecting Far-Field Boundary Conditions for Unsteady Transonic Flow Computation

The approximate nonreflecting far-field boundary condition, as proposed by Engquisi and Majda, is implemented In the computer code LTRAN2. This code solves the Implicit finite-difference representation of the small-disturbance equations for unsteady transonic flows about airfoils. The nonreflecting boundary condition and the description of the algorithm for Implementing these conditions In LTRAN2 tire discussed. Various cases re computed and compared with results from the older, more conventional procedures. One concludes that the nonreflecting far-field boundary approximation allows the far-field boundary to be located closer to the airfoil; this permits a decrease in the computer lime required to obtain the solution through the use of fewer mesh points.

Kwak, D.↗
DOE OSTI2021

Periodic boundary conditions for arbitrary deformations in molecular dynamics simulations

A generalization of the Lees-Edwards periodic boundary conditions (gLE-PBC) for molecular dynamics (MD) simulations is developed to allow for arbitrary deformations to be applied to the domain. The gLE-PBC domain remains a rectangular cuboid regardless of the applied deformation in contrast with the Lagrangian-rhomboid periodic boundary conditions (LR-PBC) where the domain deforms according to the applied deformation. Furthermore, the kinematics of gLE-PBC are validated against pure shear. The gLE-PBC method for interacting systems is then validated against the LR-PBC method and analytical solutions for a solid under isotropic compression, one-dimensional shearing, three-dimensional extension and shearing and for a liquid under Couette flow. Bulk physical properties extracted from the gLE-PBC simulations agree well with values calculated from equilibrium MD simulations. Three dimensional shearing and a deformation with a full velocity gradient matrix are also simulated, showing the range of problems gLE-PBC can explore.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1994

A stable penalty method for the compressible Navier-Stokes equations. 1: Open boundary conditions

The purpose of this paper is to present asymptotically stable open boundary conditions for the numerical approximation of the compressible Navier-Stokes equations in three spatial dimensions. The treatment uses the conservation form of the Navier-Stokes equations and utilizes linearization and localization at the boundaries based on these variables. The proposed boundary conditions are applied through a penalty procedure, thus ensuring correct behavior of the scheme as the Reynolds number tends to infinity. The versatility of this method is demonstrated for the problem of a compressible flow past a circular cylinder.

Hesthaven, J. S.↗
NASA NTRS1973

Vibration of thermally stressed plates with various boundary conditions.

By discarding Lurie's (1952) assumption of mode identity, it is shown that linear theory correctly predicts the frequency of all modes of a thermally stressed cantilever plate as well as the frequency and modes of plates with other boundary conditions. The thermal stress distribution is obtained for whatever temperature distribution and boundary conditions that may be specified. Experimental results are compared to calculated results for several different plates. Boundary conditions for the plates range from a plate with edges completely clamped to a plate with edges completely free with various other combinations of mixed and uniform edge conditions. Comparison of calculated data to experimental data shows that accurate, quantitative results can be obtained from linear theory for 'as cut' real plates for a significant range of heating when the assumption of mode identity is discarded.

Bailey, C. D.↗
NASA NTRSDate not supplied

An Extension of the Truncated Ingard-Myers Impedance Boundary Condition for High Mach Number Grazing Flows

The Ingard-Myers impedance boundary condition is widely recognized to be prone to hydrodynamic Kelvin-Helmholtz-type instability, primarily due to its use of a vortex sheet to model the flow at the boundary. Recently, a stabilized approximation of the Ingard-Myers condition, called the Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC), has been introduced. This approximation offers a good representation of the Ingard-Myers condition for grazing flows characterized by low to mid subsonic Mach numbers. This paper explores an extension of the TIMIBC, referred to as TIMIBC-ext, designed for liners in grazing flows with high subsonic Mach numbers. The TIMIBC-ext introduces a tunable parameter and demonstrates that by judiciously selecting the parameter value, the accuracy of the TIMIBC can be improved for either upstream or downstream propagating waves in grazing flows with high subsonic Mach numbers. Consequently, the TIMIBC-ext represents an enhancement to the TIMIBC when the location of the source relative to the liner (upstream or downstream) is known. Time domain implementations of the TIMIBC-ext are also provided, employing a multipole expansion model for the impedance function. A numerical example is provided to assess the performance of the TIMIBC-ext, comparing computational results with measurements from a recent set of measurements in the NASA Langley Grazing Flow Impedance Tube.

time domain acoustic scattering↗
NASA NTRS1996

On the Effective Construction of Compactly Supported Wavelets Satisfying Homogenous Boundary Conditions on the Interval

We construct compactly supported wavelet bases satisfying homogeneous boundary conditions on the interval (0,1). The maximum features of multiresolution analysis on the line are retained, including polynomial approximation and tree algorithms. The case of H(sub 0)(sup 1)(0, 1)is detailed, and numerical values, required for the implementation, are provided for the Neumann and Dirichlet boundary conditions.

Chiavassa, G.↗
NASA NTRS2004

Unsteady Validation of a Mean Flow Boundary Condition for Computational Aeroacoustics

In this work, a previously developed mean flow boundary condition will be validated for unsteady flows. The test cases will be several reference benchmark flows consisting of vortical gusts convecting in a uniform mean flow, as well as the more realistic case of a vortical gust impinging on a loaded 2D cascade. The results will verify that the mean flow boundary condition both imposes the desired mean flow as well as having little or no effect on the instantaneous unsteady solution.

Hixon, R.↗
DOE OSTI2024

Conformal boundary conditions for a 4d scalar field

We construct unitary, stable, and interacting conformal boundary conditions for a free massless scalar in four dimensions by coupling it to edge modes living on a boundary. The boundary theories we consider are bosonic and fermionic QED _3 3 with N_f N f flavors and a Chern-Simons term at level k k , in the large- N_f N f limit with fixed k/N_f k / N f . We find that interacting boundary conditions only exist when k≠ 0 k ≠ 0 . To obtain this result we compute the \beta β functions of the classically marginal couplings at the first non-vanishing order in the large- N_f N f expansion, and to all orders in k/N_f k / N f and in the couplings. To check vacuum stability we also compute the large- N_f N f effective potential. We compare our results with the the known conformal bootstrap bounds.

Di Pietro, Lorenzo (ORCID:0000000174965276)↗
NASA NTRS1995

Solution of the three-dimensional Helmholtz equation with nonlocal boundary conditions

The Helmholtz equation is solved within a three-dimensional rectangular duct with a nonlocal radiation boundary condition at the duct exit plane. This condition accurately models the acoustic admittance at an arbitrarily-located computational boundary plane. A linear system of equations is constructed with second-order central differences for the Helmholtz operator and second-order backward differences for both local admittance conditions and the gradient term in the nonlocal radiation boundary condition. The resulting matrix equation is large, sparse, and non-Hermitian. The size and structure of the matrix makes direct solution techniques impractical; as a result, a nonstationary iterative technique is used for its solution. The theory behind the nonstationary technique is reviewed, and numerical results are presented for radiation from both a point source and a planar acoustic source. The solutions with the nonlocal boundary conditions are invariant to the location of the computational boundary, and the same nonlocal conditions are valid for all solutions. The nonlocal conditions thus provide a means of minimizing the size of three-dimensional computational domains.

Hodge, Steve L.↗
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↗
DOE OSTI2014

Plasmon excitation in metal slab by fast point charge: The role of additional boundary conditions in quantum hydrodynamic model

We study the wake effect in the induced potential and the stopping power due to plasmon excitation in a metal slab by a point charge moving inside the slab. Nonlocal effects in the response of the electron gas in the metal are described by a quantum hydrodynamic model, where the equation of electronic motion contains both a quantum pressure term and a gradient correction from the Bohm quantum potential, resulting in a fourth-order differential equation for the perturbed electron density. Thus, besides using the condition that the normal component of the electron velocity should vanish at the impenetrable boundary of the metal, a consistent inclusion of the gradient correction is shown to introduce two possibilities for an additional boundary condition for the perturbed electron density. We show that using two different sets of boundary conditions only gives rise to differences in the wake potential at large distances behind the charged particle. On the other hand, the gradient correction in the quantum hydrodynamic model is seen to cause a reduction in the depth of the potential well closest to the particle, and a reduction of its stopping power. Even for a particle moving in the center of the slab, we observe nonlocal effects in the induced potential and the stopping power due to reduction of the slab thickness, which arise from the gradient correction in the quantum hydrodynamic model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗