Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary value problems”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems

In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Projection-Based Model Reduction for Coupled Conduction—Enclosure Radiation Systems

We report a projection-based reduced order model (pROM) methodology has been developed for transient heat transfer problems involving coupled conduction and enclosure radiation. The approach was demonstrated on two test problems of varying complexity. The reduced order models demonstrated substantial speedups (up to 185×) relative to the full order model with good accuracy (less than 3% L ∞ error). An attractive feature of pROMs is that there is a natural error indicator for the ROM solution: the final residual norm at each time-step of the converged ROM solution. Using example test cases, we discuss how to interpret this error indicator to assess the accuracy of the ROM solution. The approach shows promise for many-query applications, such as uncertainty quantification and optimization. The reduced computational cost of the ROM relative to the full-order model (FOM) can enable the analysis of larger and more complex systems as well as the exploration of larger parameter spaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Stability of Solutions of Initial Boundary-Balue Problems of Aerohydroelasticity

In designing structures and devices interacting with the flow of gas or liquid, it is necessary to solve the problems associated with the investigation of the stability required for their functioning and operational reliability. The definition of stability of an elastic body taken in the article corresponds to Lyapunov’s concept of stability of a dynamical system. On the base of the proposed nonlinear mathematical model the dynamic stability of the elastic aileron of the wing taking into account the incident subsonic flow of gas or liquid (in an ideal model of an incompressible environment) is investigated. Also a nonlinear mathematical model of the device relating to the vibration technique, which is intended for intensification of technological processes, for example, the process of mixing, is considered. The action of these devices is based on the oscillations of elastic elements during the gas or liquid flow. The dynamic stability of the elastic element located on one wall of the flow channel with the subsonic flow of gas or liquid (in an ideal model of a compressible environment) is investigated. Both models are described by a coupled nonlinear system of differential equations for the unknown functions — the potential of the gas velocity and deformation of the elastic element. On the basis of the construction of functionals, we obtain the sufficient conditions of the stability, the imposed restrictions on the free-stream velocity of the gas, the flexural stiffness of the elastic element, and other parameters of the mechanical system. The examples of construction of the stability regions for particular parameters of the mechanical system are presented.

42 ENGINEERING↗
DOE OSTI2018

Optimal Portfolio Management in a Modified Constant Elasticity of Variance Model

We investigate the optimal management of a portfolio consisting of a risk asset and a risk-free bond. The asset dynamics is specified by the M-CEV model. Assuming that the investor has a power utility function, we derive an explicit analytical portfolio-management formula that contains confluent hypergeometric functions. The asymptotic form of the proposed portfolio management strategy is obtained, as well as approximate formulas consisting only of elementary functions. We additionally present an application of our results in the context of algorithmic statistical-arbitrage strategies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Unique Determination of 3-Connected Plane Domains by Relative Conformal Moduli of Pairs of Boundary Components

We study the boundary values of conformal mappings of plane finitely connected domains. We present a complete exposition of a theorem giving an exhaustive description of the boundary values of such mappings in terms of conformal moduli (extremal lengths) of pairs of boundary components of a domain of connectedness at most 3. Bibliography: 12 titles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

The influence of interval arithmetic on the shape of uncertainly defined domains modelled by closed curves

The paper presents an unambiguous modelling of imprecisely defined shapes of closed curves using classical and directed interval arithmetic. The authors focus on the development of an effective strategy of modelling interval smooth closed curves (which enforce C{sup 2} continuity in points at which adjacent interval segments join) using interval cubic Bézier segments. For this purpose, algebraic relationships between Bézier and de Boor control points, formerly known for precisely defined curves, are generalized. We obtain interval control points that define interval closed curves. Additionally, the reliability of such way of modelling of closed curves is examined. We directly apply classical and directed interval arithmetic to mentioned relationships and try to solve obtained interval systems of algebraic equations. However, we obtain ambiguous solutions. Therefore, we propose our new strategy of modification of directed interval arithmetic to obtain reliable and unambiguous shapes of interval closed curves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Reconstructing the membrane detection of a 1D electrostatic-driven MEMS device by the shooting method: convergence analysis and ghost solutions identification

In this paper, in the domain of 1D-membrane micro-electro-mechanical systems in which the electrostatic field is expressed in terms of geometric curvature of the membrane, we present a numerical approach based on shooting techniques to reconstruct the membrane profile in the device in steady-state case. In particular, starting from known results in literature about existence achieved by Schauder–Tychonoff’s fixed point approach and uniqueness, and focusing on two physical–mathematical parameters appropriately indicative of the applied voltage and electromechanical properties of the membrane, respectively, we will discuss what operation parameters (applied voltage, amplitude of electrostatic field) and for which electromechanical membrane characteristic of the device is permitted or not a convergence of the method with respect to analytical results. Finally, we will discuss in detail the detected ghost solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Variance reduction through robust design of boundary conditions for stochastic hyperbolic systems of equations

We consider a hyperbolic system with uncertainty in the boundary and initial data. Our aim is to show that different boundary conditions give different convergence rates of the variance of the solution. This means that we can with the same knowledge of data get a more or less accurate description of the uncertainty in the solution. A variety of boundary conditions are compared and both analytical and numerical estimates of the variance of the solution are presented. As an application, we study the effect of this technique on Maxwell's equations as well as on a subsonic outflow boundary for the Euler equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A regular version of Smilansky model

We discuss a modification of Smilansky model in which a singular potential “channel” is replaced by a regular, below unbounded potential which shrinks as it becomes deeper. We demonstrate that, similarly to the original model, such a system exhibits a spectral transition with respect to the coupling constant, and determine the critical value above which a new spectral branch opens. The result is generalized to situations with multiple potential “channels.”.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Bell-Plesset effects in Rayleigh-Taylor instability of finite-thickness spherical and cylindrical shells

Bell-Plesset (BP) effects account for the influence of global convergence or divergence of the fluid flow on the evolution of the interfacial perturbations embedded in the flow. The development of the Rayleigh-Taylor instability in radiation-driven spherical capsules and magnetically-driven cylindrical liners necessarily includes a significant contribution from BP effects due to the time dependence of the radius, velocity, and acceleration of the unstable surfaces or interfaces. An analytical model is presented that, for an ideal incompressible fluid and small perturbation amplitudes, exactly evaluates the BP effects in finite-thickness shells through acceleration and deceleration phases. The time-dependent dispersion equations determining the “instantaneous growth rate” are derived. It is demonstrated that by integrating this approximate growth rate over time, one can accurately evaluate the number of perturbation e-foldings during the inward acceleration phase of the implosion. As a result, in the limit of small shell thickness, exact thin-shell perturbationequations and approximate thin-shell dispersion equations are obtained, generalizing the earlier results [E. G. Harris, Phys. Fluids 5, 1057 (1962); E. Ott, Phys. Rev. Lett. 29, 1429 (1972); A. B. Bud'ko et al., Phys. Fluids B 2, 1159 (1990)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Two-fluid and finite Larmor radius effects on helicity evolution in a plasma pinch

The evolution of magnetic energy, helicity, and hybrid helicity during nonlinear relaxation of a driven-damped plasma pinch is compared in visco-resistive magnetohydrodynamics and two-fluid models with and without the ion gyroviscous stress tensor. Magnetic energy and helicity are supplied via a boundary electric field which initially balances the resistive dissipation, and the plasma undergoes multiple relaxation events during the nonlinear evolution. The magnetic helicity is well conserved relative to the magnetic energy over each event, which is short compared with the global resistive diffusion time. The magnetic energy decreases by roughly 1.5% of its initial value over a relaxation event, while the magnetic helicity changes by at most 0.2% of the initial value. The hybrid helicity is dominated by magnetic helicity in low-β pinch conditions and is also well conserved. Differences of less than 1% between magnetic helicity and hybrid helicity are observed with two-fluid modeling and result from cross helicity evolution. The cross helicity is found to change appreciably due to the first-order finite Larmor radius effects which have not been included in contemporary relaxation theories. The plasma current evolves towards the flat parallel current state predicted by Taylor relaxation theory but does not achieve it. Plasma flow develops significant structure for two-fluid models, and the flow perpendicular to the magnetic field is much more substantial than the flow along it.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Structure of cyano-anion ionic liquids: X-ray scattering and simulations

Ionic liquids with cyano anions have long been used because of their unique combination of low-melting temperatures, reduced viscosities, and increased conductivities. Recently we showed that cyano anions in ionic liquids are particularly interesting for their potential use as electron donors to excited state photo-acceptors [B. Wu et al., J. Phys. Chem. B 119, 14790–14799 (2015)]. Here we report on bulk structural and quantum mechanical results for a series of ionic liquids based on the 1-ethyl-3-methylimidazolium cation, paired with the following five cyano anions: SeCN - , SCN - , N(CN) $-\atop{2}$ , C(CN) $-\atop{3}$ , and B(CN) $-\atop{4}$ . By combining molecular dynamics simulations, high-energy X-ray scattering measurements, and periodic boundary condition DFT calculations, we are able to obtain a comprehensive description of the liquid landscape as well as the nature of the HOMO-LUMO states for these ionic liquids in the condensed phase. Features in the structure functions for these ionic liquids are somewhat different than the commonly observed adjacency, charge-charge, and polarity peaks, especially for the bulkiest B(CN) $-\atop{4}$ anion. While the other four cyano-anion ionic liquids present an anionic HOMO, the one for Im $+\atop{2,1}$ /B(CN) $-\atop{4}$ is cationic.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Plasma electron hole kinematics. II. Hole tracking Particle-In-Cell simulation

The kinematics of a 1-D electron hole is studied using a novel Particle-In-Cell simulation code. A hole tracking technique enables us to follow the trajectory of a fast-moving solitary hole and study quantitatively hole acceleration and coupling to ions. We observe a transient at the initial stage of hole formation when the hole accelerates to several times the cold-ion sound speed. Artificially imposing slow ion speed changes on a fully formed hole causes its velocity to change even when the ion stream speed in the hole frame greatly exceeds the ion thermal speed, so there are no reflected ions. The behavior that we observe in numerical simulations agrees very well with our analytic theory of hole momentum conservation and the effects of “jetting.”.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Real-space pseudopotential method for computing the vibrational Stark effect

The vibrational Stark shift is an important effect in determining the electrostatic environment for molecular or condensed matter systems. However, accurate ab initio calculations of the vibrational Stark effect are a technically demanding challenge. In this paper, we make use of density functional theory constructed on a real-space grid to expedite the computation of this effect. Our format is especially advantageous for the investigation of small molecules in finite fields as cluster boundary conditions eliminate spurious supercell interactions and allow for charged systems, while convergence is controlled by a single parameter, the grid spacing. The Stark tuning rate is highly sensitive to the interaction between anharmonicity in a vibrational mode and the applied field. To ensure this subtle interaction is fully captured, we apply three parallel approaches: a direct finite field, a perturbative method, and a molecular dynamics method. Finally, we illustrate this method by applying it to several small molecules containing C–O and C–N bonds and show that a consistent result can be obtained.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

A Semi-Analytical Solution Approach for Solving Constant-Coefficient First-Order Partial Differential Equations

Simulation and control of many dynamic systems involve solving partial differential equations (PDE). This letter proposes a semi-analytical solution (SAS) approach for fast and high-quality solution of first-order PDEs. The region of interest of the studied PDE is divided into a grid, and an SAS is derived for each grid cell in the form of the multivariate polynomials, of which the coefficients are identified using initial value and boundary value conditions. The solutions are solved in a “time-stepping” manner, i.e. within one time step, the coefficients of the SAS are identified and the initial value of the next time step is evaluated. This approach achieves a significantly larger grid cell than the widely used finite difference method, and thus enhances the computational efficiency significantly. Furthermore, the simulation result on the natural gas pipeline model demonstrates the advantages of SAS in accuracy and computational efficiency.

97 MATHEMATICS AND COMPUTING↗