Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unbounded domain”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Generalized Symbolic Execution for Model Checking and Testing

Modern software systems, which often are concurrent and manipulate complex data structures must be extremely reliable. We present a novel framework based on symbolic execution, for automated checking of such systems. We provide a two-fold generalization of traditional symbolic execution based approaches: one, we define a program instrumentation, which enables standard model checkers to perform symbolic execution; two, we give a novel symbolic execution algorithm that handles dynamically allocated structures (e.g., lists and trees), method preconditions (e.g., acyclicity of lists), data (e.g., integers and strings) and concurrency. The program instrumentation enables a model checker to automatically explore program heap configurations (using a systematic treatment of aliasing) and manipulate logical formulae on program data values (using a decision procedure). We illustrate two applications of our framework: checking correctness of multi-threaded programs that take inputs from unbounded domains with complex structure and generation of non-isomorphic test inputs that satisfy a testing criterion. Our implementation for Java uses the Java PathFinder model checker.

Khurshid, Sarfraz↗
DOE OSTI2024

Lattice Green’s Functions for High-Order Finite Difference Stencils

Lattice Green's Functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Shock/vortex interaction and vortex-breakdown modes

Computational simulation and study of shock/vortex interaction and vortex-breakdown modes are considered for bound (internal) and unbound (external) flow domains. The problem is formulated using the unsteady, compressible, full Navier-Stokes (NS) equations which are solved using an implicit, flux-difference splitting, finite-volume scheme. For the bound flow domain, a supersonic swirling flow is considered in a configured circular duct and the problem is solved for quasi-axisymmetric and three-dimensional flows. For the unbound domain, a supersonic swirling flow issued from a nozzle into a uniform supersonic flow of lower Mach number is considered for quasi-axisymmetric and three-dimensional flows. The results show several modes of breakdown; e.g., no-breakdown, transient single-bubble breakdown, transient multi-bubble breakdown, periodic multi-bubble multi-frequency breakdown and helical breakdown.

Kandil, Osama A.↗
NASA NTRS1992

Modal element method for scattering of sound by absorbing bodies

The modal element method for acoustic scattering from 2-D body is presented. The body may be acoustically soft (absorbing) or hard (reflecting). The infinite computational region is divided into two subdomains - the bounded finite element domain, which is characterized by complicated geometry and/or variable material properties, and the surrounding unbounded homogeneous domain. The acoustic pressure field is represented approximately in the finite element domain by a finite element solution, and is represented analytically by an eigenfunction expansion in the homogeneous domain. The two solutions are coupled by the continuity of pressure and velocity across the interface between the two subdomains. Also, for hard bodies, a compact modal ring grid system is introduced for which computing requirements are drastically reduced. Analysis for 2-D scattering from solid and coated (acoustically treated) bodies is presented, and several simple numerical examples are discussed. In addition, criteria are presented for determining the number of modes to accurately resolve the scattered pressure field from a solid cylinder as a function of the frequency of the incoming wave and the radius of the cylinder.

Baumeister, Kenneth J.↗
NASA NTRS1993

Modal element method for scattering and absorbing of sound by two-dimensional bodies

The modal element method for acoustic scattering from a 2-D body is presented. The body may be acoustically soft (absorbing) or hard (reflecting). The infinite computational region is divided into two subdomains - the bounded finite element domain, which is characterized by complicated geometry and/or variable material properties, and the surrounding unbounded homogeneous domain. The acoustic pressure field is represented approximately in the finite element domain by a finite element solution, and is represented analytically by an eigenfunction expansion in the homogeneous domain. The two solutions are coupled by the continuity of pressure and velocity across the interface between the two subdomains. Also, for hard bodies, a compact modal ring grid system is introduced for which computing requirements are drastically reduced. Analysis for 2-D scattering from solid and coated (acoustically treated) bodies is presented, and several simple numerical examples are discussed. In addition, criteria are presented for determining the number of modes to accurately resolve the scattered pressure field from a solid cylinder as a function of the frequency of the incoming wave and the radius of the cylinder.

Baumeister, K. J.↗
NASA NTRS2000

High-order Two-way Artificial Boundary Conditions for Nonlinear Wave Propagation with Backscattering

When solving linear scattering problems, one typically first solves for the impinging wave in the absence of obstacles. Then, by linear superposition, the original problem is reduced to one that involves only the scattered waves driven by the values of the impinging field at the surface of the obstacles. In addition, when the original domain is unbounded, special artificial boundary conditions (ABCs) that would guarantee the reflectionless propagation of waves have to be set at the outer boundary of the finite computational domain. The situation becomes conceptually different when the propagation equation is nonlinear. In this case the impinging and scattered waves can no longer be separated, and the problem has to be solved in its entirety. In particular, the boundary on which the incoming field values are prescribed, should transmit the given incoming waves in one direction and simultaneously be transparent to all the outgoing waves that travel in the opposite direction. We call this type of boundary conditions two-way ABCs. In the paper, we construct the two-way ABCs for the nonlinear Helmholtz equation that models the laser beam propagation in a medium with nonlinear index of refraction. In this case, the forward propagation is accompanied by backscattering, i.e., generation of waves in the direction opposite to that of the incoming signal. Our two-way ABCs generate no reflection of the backscattered waves and at the same time impose the correct values of the incoming wave. The ABCs are obtained for a fourth-order accurate discretization to the Helmholtz operator; the fourth-order grid convergence is corroborated experimentally by solving linear model problems. We also present solutions in the nonlinear case using the two-way ABC which, unlike the traditional Dirichlet boundary condition, allows for direct calculation of the magnitude of backscattering.

Fibich, Gadi↗
NASA NTRS1997

An Application of the Difference Potentials Method to Solving External Problems in CFD

Numerical solution of infinite-domain boundary-value problems requires some special techniques that would make the problem available for treatment on the computer. Indeed, the problem must be discretized in a way that the computer operates with only finite amount of information. Therefore, the original infinite-domain formulation must be altered and/or augmented so that on one hand the solution is not changed (or changed slightly) and on the other hand the finite discrete formulation becomes available. One widely used approach to constructing such discretizations consists of truncating the unbounded original domain and then setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The role of the ABC's is to close the truncated problem and at the same time to ensure that the solution found inside the finite computational domain would be maximally close to (in the ideal case, exactly the same as) the corresponding fragment of the original infinite-domain solution. Let us emphasize that the proper treatment of artificial boundaries may have a profound impact on the overall quality and performance of numerical algorithms. The latter statement is corroborated by the numerous computational experiments and especially concerns the area of CFD, in which external problems present a wide class of practically important formulations. In this paper, we review some work that has been done over the recent years on constructing highly accurate nonlocal ABC's for calculation of compressible external flows. The approach is based on implementation of the generalized potentials and pseudodifferential boundary projection operators analogous to those proposed first by Calderon. The difference potentials method (DPM) by Ryaben'kii is used for the effective computation of the generalized potentials and projections. The resulting ABC's clearly outperform the existing methods from the standpoints of accuracy and robustness, in many cases noticeably speed up the multigrid convergence, and at the same time are quite comparable to other methods from the standpoints of geometric universality and simplicity of implementation.

Ryaben 'Kii, Victor S.↗
NASA NTRS2003

Computational Methods for Jet Noise Simulation

The purpose of our project is to develop, analyze, and test novel numerical technologies central to the long term goal of direct simulations of subsonic jet noise. Our current focus is on two issues: accurate, near-field domain truncations and high-order, single-step discretizations of the governing equations. The Direct Numerical Simulation (DNS) of jet noise poses a number of extreme challenges to computational technique. In particular, the problem involves multiple temporal and spatial scales as well as flow instabilities and is posed on an unbounded spatial domain. Moreover, the basic phenomenon of interest, the radiation of acoustic waves to the far field, involves only a minuscule fraction of the total energy. The best current simulations of jet noise are at low Reynolds number. It is likely that an increase of one to two orders of magnitude will be necessary to reach a regime where the separation between the energy-containing and dissipation scales is sufficient to make the radiated noise essentially independent of the Reynolds number. Such an increase in resolution cannot be obtained in the near future solely through increases in computing power. Therefore, new numerical methodologies of maximal efficiency and accuracy are required.

Goodrich, John W.↗
NASA NTRS1989

Interaction dynamics of multiple mobile robots with simple navigation strategies

The global dynamic behavior of multiple interacting autonomous mobile robots with simple navigation strategies is studied. Here, the effective spatial domain of each robot is taken to be a closed ball about its mass center. It is assumed that each robot has a specified cone of visibility such that interaction with other robots takes place only when they enter its visibility cone. Based on a particle model for the robots, various simple homing and collision-avoidance navigation strategies are derived. Then, an analysis of the dynamical behavior of the interacting robots in unbounded spatial domains is made. The article concludes with the results of computer simulations studies of two or more interacting robots.

Wang, P. K. C.↗
NASA NTRS1989

Blow-up of unsteady two-dimensional Euler and Navier-Stokes solutions having stagnation-point form

The time-dependent form of the classic, two-dimensional stagnation-point solution of the Navier-Stokes equations is considered. If the viscosity is zero, a class of solutions of the initial-value problem can be found in closed form using Lagrangian coordinates. These solutions exhibit singular behavior in finite time, because of the infinite domain and unbounded initial vorticity. Thus, the blow-up found by Stuart in three dimensions using the stagnation-point form, also occurs in two. The singularity vanishes under a discrete, finite-dimensional 'point vortex' approximation, but is recovered as the number of vortices tends to infinity. We find that a small positive viscosity does not arrest the breakdown, but does strongly alter its form. Similar results are summarized for certain Boussinesq stratified flows.

Childress, S.↗
NASA NTRS1997

External Boundary Conditions for Three-Dimensional Problems of Computational Aerodynamics

We consider an unbounded steady-state flow of viscous fluid over a three-dimensional finite body or configuration of bodies. For the purpose of solving this flow problem numerically, we discretize the governing equations (Navier-Stokes) on a finite-difference grid. The grid obviously cannot stretch from the body up to infinity, because the number of the discrete variables in that case would not be finite. Therefore, prior to the discretization we truncate the original unbounded flow domain by introducing some artificial computational boundary at a finite distance of the body. Typically, the artificial boundary is introduced in a natural way as the external boundary of the domain covered by the grid. The flow problem formulated only on the finite computational domain rather than on the original infinite domain is clearly subdefinite unless some artificial boundary conditions (ABC's) are specified at the external computational boundary. Similarly, the discretized flow problem is subdefinite (i.e., lacks equations with respect to unknowns) unless a special closing procedure is implemented at this artificial boundary. The closing procedure in the discrete case is called the ABC's as well. In this paper, we present an innovative approach to constructing highly accurate ABC's for three-dimensional flow computations. The approach extends our previous technique developed for the two-dimensional case; it employs the finite-difference counterparts to Calderon's pseudodifferential boundary projections calculated in the framework of the difference potentials method (DPM) by Ryaben'kii. The resulting ABC's appear spatially nonlocal but particularly easy to implement along with the existing solvers. The new boundary conditions have been successfully combined with the NASA-developed production code TLNS3D and used for the analysis of wing-shaped configurations in subsonic (including incompressible limit) and transonic flow regimes. As demonstrated by the computational experiments and comparisons with the standard (local) methods, the DPM-based ABC's allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable increase of the convergence rate of multigrid iterations.

Tsynkov, Semyon V.↗
NASA NTRS1994

Modal element method for potential flow in non-uniform ducts: Combining closed form analysis with CFD

An analytical procedure is presented, called the modal element method, that combines numerical grid based algorithms with eigenfunction expansions developed by separation of variables. A modal element method is presented for solving potential flow in a channel with two-dimensional cylindrical like obstacles. The infinite computational region is divided into three subdomains; the bounded finite element domain, which is characterized by the cylindrical obstacle and the surrounding unbounded uniform channel entrance and exit domains. The velocity potential is represented approximately in the grid based domain by a finite element solution and is represented analytically by an eigenfunction expansion in the uniform semi-infinite entrance and exit domains. The calculated flow fields are in excellent agreement with exact analytical solutions. By eliminating the grid surrounding the obstacle, the modal element method reduces the numerical grid size, employs a more precise far field boundary condition, as well as giving theoretical insight to the interaction of the obstacle with the mean flow. Although the analysis focuses on a specific geometry, the formulation is general and can be applied to a variety of problems as seen by a comparison to companion theories in aeroacoustics and electromagnetics.

Baumeister, Kenneth J.↗
DOE OSTI2025

Structure-function studies of Vibrio cholerae quorum-sensing receptor CqsR signal recognition

Ethanolamine signaling through the transmembrane quorum-sensing receptor CqsR influencesVibrio choleraeniche recognition and host colonization. In this study, we present a comprehensive structure-function analysis of CqsR. Specifically, we have determined X-ray crystal structures of the CqsR periplasmic domain bound to the signaling agonist ethanolamine and its analogs, serinol and L-alaninol, as well as the ligand-free (apo) form of CqsR. The periplasmic ligand-binding domain of CqsR is a Cache domain, the most prevalent extracellular sensory module in prokaryotes. Our findings provide a rare structural comparison of ligand-bound and unbound states of a Cache domain receptor. Coupled with thermodynamic binding assays and genetic analyses, these structures elucidate the molecular basis of CqsR ligand specificity. This study not only advances the understanding of Cache domain function but also informs the identification of ligands for orphan Cache receptors and the rational design of signaling agonists and antagonists. Lastly, we discuss ligand-induced conformational changes in the CqsR Cache domains and explore the potential for the existence of additional regulatory ligands.

Microbiology↗
DOE OSTI2021

Theoretical analysis of key factors achieving reversed magnetic shear q-profiles sustained with lower hybrid waves on EAST

The phase space analysis technique is applied to provide new insights into the fully non-inductive lower hybrid current driven (LHCD) discharges on EAST (Garofalo et al 2017 Nucl. Fusion 57 076037). The analysis shows that there are bounded and unbounded topologies of lower hybrid (LH) waves in phase space. For typical parameters on EAST, the propagation domain for 4.6 GHz LH waves is bounded, while unbounded for 2.45 GHz LH waves and one of the conditions is recognized to be good for achieving an off-axis current profile driven by 4.6 GHz LH waves on EAST. The parametric analysis on the potential power deposition (PPD) region for those experiments demonstrates that the reversed magnetic shear dominates the wave behavior, and confines the LH power absorption to the far off-axis region as long as it occurs in the 4.6 GHz dominated LHCD discharge. GENRAY/CQL3D simulations also confirm this effect of reversed shear on the power absorption of 4.6 GHz waves. A bounded propagation domain and a positive feedback loop between magnetic shear reversal and off-axis LHCD profile could explain the LH-only-sustained strong reversed magnetic shear observed in experiments on multiple machines. In contrast, there is no such effect for 2.45 GHz waves on EAST since the PPD region is less sensitive to the change of reversed shear due to the unbounded propagation domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analysis of nonlinear features associated with radio-frequency sheaths

An analytic study is conducted of nonlinear features that arise when radio-frequency (RF) sheath boundary conditions are solved self-consistently with the electromagnetic fields in a magnetized plasma. An instability associated with RF sheath–plasma interactions is identified and shown to be closely linked with sheath–plasma resonance and the occurrence of multiple solutions. It is proposed that this instability serves as a physical mechanism that enables jumping from one solution branch to another, providing a physical basis for the prediction of an RF sheath hysteresis in the time domain. It is further shown that this instability can lead to non-intuitive features such as spontaneous symmetry breaking in an otherwise completely symmetric problem. In plane-stratified geometry, two analytically solvable examples are considered for the cases of an unbounded half-space and a bounded domain

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Mathematical and computational studies of equilibrium capillary free surfaces

The results of several independent studies are presented. The general question is considered of whether a wetting liquid always rises higher in a small capillary tube than in a larger one, when both are dipped vertically into an infinite reservoir. An analytical investigation is initiated to determine the qualitative behavior of the family of solutions of the equilibrium capillary free-surface equation that correspond to rotationally symmetric pendent liquid drops and the relationship of these solutions to the singular solution, which corresponds to an infinite spike of liquid extending downward to infinity. The block successive overrelaxation-Newton method and the generalized conjugate gradient method are investigated for solving the capillary equation on a uniform square mesh in a square domain, including the case for which the solution is unbounded at the corners. Capillary surfaces are calculated on the ellipse, on a circle with reentrant notches, and on other irregularly shaped domains using JASON, a general purpose program for solving nonlinear elliptic equations on a nonuniform quadrilaterial mesh. Analytical estimates for the nonexistence of solutions of the equilibrium capillary free-surface equation on the ellipse in zero gravity are evaluated.

Albright, N.↗
NASA NTRS1995

Numerical Stability In Hyperbolic Boundary-Value Problems

Technical memorandum discusses stability of numerical solutions involving semidiscrete approximations to hyperbolic partial differential equations in initial-and-boundary-value problems. Topic of practical significance because hyperbolic partial differential equations arise in mathematical modeling of waves and blasts. Solutions often needed over restricted regions of unbounded spaces. Outer boundaries artificial, introduced only to limit domains of numerical solutions. Conditions at such artificial boundaries cause numerical instabilities that degrade computed solutions.

Warming, Robert F.↗
DOE OSTI2020

Inverse point source location with the Helmholtz equation on a bounded domain

The problem of recovering acoustic sources, more specifically monopoles, from point-wise measurements of the corresponding acoustic pressure at a limited number of frequencies is addressed. To this purpose, a family of sparse optimization problems in measure space in combination with the Helmholtz equation on a bounded domain is considered in this study. A weighted norm with unbounded weight near the observation points is incorporated into the formulation. Optimality conditions and conditions for recovery in the small noise case are discussed, which motivates concrete choices of the weight. The numerical realization is based on an accelerated conditional gradient method in measure space and a finite element discretization.

97 MATHEMATICS AND COMPUTING↗