Explicit synchronous partitioned algorithms for interface problems based on Lagrange multipliers
Not Available
SEARCH · Search NASA
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.
Not Available
The problem of the transmission of sound in a duct with very thin shear layers at the walls is treated by an inner expansion method. The results show that the formulation of the problem of the transmission of sound in a duct with a shear layer at the wall converges, in the case of a vanishingly thin shear layer, to the formulation of the same problem when uniform flow is assumed and the wall boundary condition is that of continuity of particle displacement.
The optimal transmission switching (OTS) problem with AC physics represents a mixed integer non-linear non-convex optimization problem which can provide benefits to transmission level power system operations. In this paper we benchmark a set of open-source mixed integer non-linear programming (MINLP) solvers on the OTS problem with AC physics using the pglib set of power system test cases. Results characterizing the performance of the different solvers are reported and discussed.
Network transmission errors such as collisions, CRC errors, misalignment, etc. are statistical in nature. Although errors can vary randomly, a high level of errors does indicate specific network problems, e.g. equipment failure. In this project, we have studied the random nature of collisions theoretically as well as by gathering statistics, and established a numerical threshold above which a network problem is indicated with high probability.
Recent experience using high definition video on the International Space Station reveals camera pixel degradation due to particle radiation to be a much more significant problem with high definition cameras than with standard definition video. Although it may at first appear that increased pixel density on the imager is the logical explanation for this, the ISS implementations of high definition suggest a more complex causal and mediating factor mix. The degree of damage seems to vary from one type of camera to another, and this variation prompts a reconsideration of the possible factors in pixel loss, such as imager size, number of pixels, pixel aperture ratio, imager type (CCD or CMOS), method of error correction/concealment, and the method of compression used for recording or transmission. The problem of imager pixel loss due to particle radiation is not limited to out-of-atmosphere applications. Since particle radiation increases with altitude, it is not surprising to find anecdotal evidence that video cameras subject to many hours of airline travel show an increased incidence of pixel loss. This is even evident in some standard definition video applications, and pixel loss due to particle radiation only stands to become a more salient issue considering the continued diffusion of high definition video cameras in the marketplace.
Loosely coupled partitioned methods for multiphysics problems treat each subproblem as a separate entity and advance them independently in time. In so doing these methods enable code reuse, increase concurrency and provide a convenient framework for plug-and-play multiphysics simulations. However, mathematically loosely coupled schemes are equivalent to a single step of an iterative solution method, which can compromise their accuracy and stability. We present a new data-driven partitioned method for coupled parametric PDEs that can improve upon the accuracy of traditional loosely coupled methods without incurring a performance penalty. To that end, we replace conventional field transfers across the interface by a surrogate for the dynamics of the interface flux exchanged between the subdomains. To develop this surrogate we apply dynamic mode decomposition to a non-standard staggered-in-time state, comprising the interface flux and small solution patches near the interface. The new approach shifts the main computational burden to an offline training phase, whereas application of the surrogate in the online phase amounts to a single matrix–vector multiplication. In conclusion, we provide stability analysis of the surrogate-based partitioned scheme and include numerical results that demonstrate its potential.
The production of noise by turbulence in jets is an extremely complex problem. One aspect of that problem, the transmission of acoustic disturbances from the interior of the jet through the mean velocity profile and into the far field is studied. The jet (two-dimensional or circular cylindrical) is assumed infinitely long with mean velocity profile independent of streamwise location. The noise generator is a sequence of transient sources drifting with the surrounding fluid and confined to a short length of the jet.
This paper solves the problem of sound transmission through a system of two infinite concentric cylindrical sandwich shells. The shells are surrounded by external and internal fluid media and there is fluid (air) in the annular space between them. An oblique plane sound wave is incident upon the surface of the outer shell. A uniform flow is moving with a constant velocity in the external fluid medium. Classical thin shell theory is applied to the inner shell and first-order shear deformation theory is applied to the outer shell. A closed form for transmission loss is derived based on modal analysis. Investigations have been made for the impedance of both shells and the transmission loss through the shells from the exterior into the interior. Results are compared for double sandwich shells and single sandwich shells. This study shows that: (1) the impedance of the inner shell is much smaller than that of the outer shell so that the transmission loss is almost the same in both the annular space and the interior cavity of the shells; (2) the two concentric sandwich shells can produce an appreciable increase of transmission loss over single sandwich shells especially in the high frequency range; and (3) design guidelines may be derived with respect to the noise reduction requirement and the pressure in the annular space at a mid-frequency range.
With increasingly severe wildfire conditions driven by climate change, utilities must manage the risk of wildfire ignitions from electric power lines. During "public safety power shutoff'" events, utilities de-energize power lines to reduce wildfire ignition risk, which may result in load shedding. Distributed energy resources provide flexibility that can help support the system to reduce load shedding when lines are de-energized. We investigate a coordinated transmission-distribution optimization problem that balances wildfire risk mitigation and load shedding. We model distribution systems that include battery energy storage systems which may support loads when transmission lines are de-energized. This multi-period integrated transmission-distribution optimal switching problem jointly optimizes line switching decisions, the generators' setpoints, load shedding, and the batteries' states of charge, resulting in significant computational challenges. To improve scalability, we decompose the problem over both space and time and apply a distributed optimization algorithm. Using a large-scale synthetic California test case with realistic distribution models and real wildfire risk data, we show that distributed optimization can solve large-scale multi-period switching problems that are otherwise intractable for centralized solvers. We also discuss challenges and future directions for improving the distributed algorithm's convergence performance as the number of time periods increases.
The evaluation of the tensile strength of impact built-up ice on structural components has been prompted by such problems as electrical transmission line losses and catastrophic failures in Arctic regions, deicing problems with fixed-wing and rotary-wing aircraft, etc. It is demonstrated that the conventional tensile-testing technique furnishes adequate data on artificially refrigerated ice, and helps establish the influence of extrinsic factors on ice tensile strength.
The interaction of a vibrating plate with an adjacent acoustic medium is important in problems involving the radiation of sound from panels, in problems involving the transmission of sound through walls of buildings, aircraft, or launch vehicles; and in problems involving the estimation of damping and the stress amplitude of vibration for panel-fatigue predictions. There appear to have been no systematic studies of the effects on the plate of fluid coupling for an arbitrary fluid-mass/plate-mass loading ratio. An attempt is made to determine this effect for a wide range of fluid-plate mass ratios without resorting to the usual simplifications of light or heavy fluid loading. Emphasis is with the plate motion rather than the radiation of sound.
Two different aspects of active antenna array design were investigated. The transition between monolithic microwave integrated circuits and rectangular waveguides was studied along with crosstalk in multiconductor transmission lines. The boundary value problem associated with a discontinuity in a microstrip line is formulated. This entailed, as a first step, the derivation of the propagating as well as evanescent modes of a microstrip line. The solution is derived to a simple discontinuity problem: change in width of the center strip. As for the multiconductor transmission line problem. A computer algorithm was developed for computing the crosstalk noise from the signal to the sense lines. The computation is based on the assumption that these lines are terminated in passive loads.
Characterizing physical properties of faults, such as their transmissibility, is crucial for performing predictive numerical simulation of subsurface flows, such as those encountered in petroleum engineering and remediation of subsurface contamination. Here, this paper provides a complete investigation of the inverse problem for fault transmissibility in subsurface flow models, under appropriate assumptions on fault structure. In particular, the following aspects are considered: 1) fault modeling and well-posedness of the forward problem; 2) finite element (FEM) discretizations of the forward problem and their rigorous a priori convergence analysis; 3) Well-posedness of the Bayesian inverse problem, FEM discretization of the infinite dimensional Bayesian inverse formulation, and its rigorous a priori analysis. Moreover, computation of the maximum a posteriori (MAP) point via fast inexact Newton-conjugate gradient optimization and a Laplace approximation of the Bayesian posterior are also presented. Numerical results illustrate the use of the proposed fault model in forward and inverse problems for subsurface flows in two dimensional domains with multiple faults.
The main ideas of Shannon's (1948, 1960) theory of source encoding with a fidelity constraint, more commonly known as rate distortion theory, are summarized. The theory was specifically intended to provide a theoretical basis for efficient transmission of information such as images. What the theory has to contribute to the problem is demonstrated. Difficulties that impeded application of the theory to image transmission, and current efforts to solve these difficulties are discussed.
The problem of acoustic transmission through nonuniform ducts containing a high-speed subsonic flow is studied by means of the method of weighted residuals in the form of a modified Galerkin method and a Galerkin formulation of the finite element method. The method of weighted residuals is shown to employ the basis functions generated from eigenvalue calculations for the case of no flow, and is verified by comparison with exact eigenvalue calculations in the uniform duct case and numerical solutions of the one-dimensional form of the equations in the nonuniform duct case. The finite element scheme based on both the Galerkin method and the residual least squares method and employing eight-noded isoparametric elements is presented and used to investigate multimodal propagation by the coupling of the solution in the duct nonuniform section to modal expansions in uniform sections. Comparison of the results of the two methods reveals them to be in substantial agreement, and predicts the importance of multimodal interactions at high Mach numbers.
The Method of Weighted Residuals (MWR) and the Finite Element Method (FEM) are considered as computational schemes in the problem of acoustic transmission in nonuniform ducts. MWR is presented in an improved form which includes the interaction of acoustic modes (irrotational) and hydrodynamic modes (rotational). FEM is based on a weighted residuals formulation using eight noded isoparametric elements. Both are applicable to two-dimensional and axially symmetric problems. Calculations are made for several sample problems to demonstrate accuracy and relative efficiency. One test case has implications in the phenomenon of subsonic acoustic choking and it is found that a large transmission loss is not an automatic consequence of propagation against a high subsonic mean flow.
The extended method of equivalent inclusions is applied to study the specific wave problems: (1) the transmission of elastic waves in an infinite medium containing a layer of inhomogeneity, and (2) the scattering of elastic waves in an infinite medium containing a perfect spherical inhomogeneity. Eigenstrains are expanded as a geometric series and a method of integration based on the inhomogeneous Helmholtz operator is adopted. This study compares results, obtained by using limited number of terms in the eigenstrain expansion, with exact solutions for the layer problem and that for a perfect sphere.
The extended method of equivalent inclusion developed is applied to study the specific wave problems of the transmission of elastic waves in an infinite medium containing a layer of inhomogeneity, and of the scattering of elastic waves in an infinite medium containing a perfect spherical inhomogeneity. The eigenstrains are expanded as a geometric series and the method of integration for the inhomogeneous Helmholtz operator given by Fu and Mura is adopted. The results obtained by using a limited number of terms in the eigenstrain expansion are compared with exact solutions for the layer problem and for a perfect sphere. Two parameters are singled out for this comparison: the ratio of elastic moduli, and the ratio of the mass densities. General trends for three different situations are shown.