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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multistage Stabilized Continuation for Indirect Optimal Control of Three-Dimensional Hypersonic Trajectories

This work presents the development of a multistaged stabilized continuation for the three-dimensional unpowered hypersonic trajectory planning problem using indirect optimal control methods. The stabilized continuation method is noniterative and guaranteed to terminate within a finite number of floating point operations, thereby making it well suited for onboard autonomous implementations. Here, we present a multistage formulation of the stabilized continuation scheme that involves starting with a “loose” integration tolerance during the first stage and ramping up toward a “strict” integration tolerance through subsequent stages. An important benefit of this approach is that even when the solution to the underlying optimal control problem is numerically unstable, such as with the hypersonic vehicle footprint generation problem, the stabilized continuation algorithm is shown to be successful in finding a solution while providing some additional insights into the underlying cause of the numerical instability.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2016

Fully nonlinear δf gyrokinetics for scrape-off layer parallel transport

Edge plasmas present a few challenges for gyrokinetic simulations that are absent in tokamak cores. Among them are large fluctuation amplitudes and plasma-wall interactions in the open field line region. In this work, the widely used core turbulence code GENE, which employs a δf -splitting technique, is extended to simulate open systems with large electrostatic fluctuations. With inclusion and proper discretization of the parallel nonlinear term, it becomes equivalent to a full- f code and the δf -splitting causes no fundamental difficulty in handling large fluctuations. The loss of particles to the wall is accounted for by using a logical sheath boundary, which is implemented in the context of a finite-volume method. The extended GENE code is benchmarked for the well-established one-dimensional parallel transport problem in the scrape-off layer during edge-localized modes. The parallel heat flux deposited onto the divertor target is compared with previous simulation results and shows good agreement.

43 PARTICLE ACCELERATORS↗
NASA NTRS1983

Altitude/path-angle transitions in fuel-optimal problems for transport aircraft

Altitude/path angle transitions to fuel-optimal, energy climb and descent paths are examined for subsonic transport aircraft. These transitions are termed boundary layers in singular perturbation theory, the framework used herein to simplify solutions to the state-Euler system. The energy climb and descent paths provide equilibrium or stationary points for the boundary layer system, and thus for moderate transitions, presumed valid for the class of aircraft considered, the boundary layer system may be linearized. This simplification allows derivation of an explicit solution to the two-point boundary-value, boundary layer problem and, as a consequence, yields a nearly optimal control law, in feedback form, for the transitions. Numerical simulation results using the feedback control law are presented for a Boeing 737 airframe (NASA-ATOPS Vehicle) employing twin JT8D-7-7A engines.

Gracey, C.↗
DOE OSTI2022

Verification of a specialized hydrodynamic simulation code for modeling deflagration and detonation of high explosives

A specialized hydrodynamic simulation code has been developed and verified for the simulation of one-dimensional unsteady problems involving the detonation and deflagration of high explosives. To model all the relevant physical processes in these problems, a code is required to simulate compressible hydrodynamics, unsteady thermal conduction, and chemical reactions with complex rate laws. Several verification exercises are presented which test the implementation of these capabilities. The code also requires models for physics processes such as equations of state and conductivity for pure materials and mixtures as well as rate laws for chemical reactions. Additional verification tests are required to ensure that these models are implemented correctly. Though this code is limited in the types of problems it can simulate, its computationally efficient formulation allows it to be used in calibration studies for reactive burn models for high explosives. Furthermore, this study demonstrates how a series of verification tests can be used to ensure that the various physics processes needed to simulate complex phenomenon can be tested to ensure that they are correctly implemented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

High-beta analytic equilibria in circular, elliptical, and D-shaped large aspect ratio axisymmetric configurations with poloidal and toroidal flows

The Grad-Shafranov-Bernoulli system of equations is a single fluid magnetohydrodynamical description of axisymmetric equilibria with mass flows. Using a variational perturbative approach [E. Hameiri, Phys. Plasmas 20, 024504 (2013)], analytic approximations for high-beta equilibria in circular, elliptical, and D-shaped cross sections in the high aspect ratio approximation are found, which include finite toroidal and poloidal flows. Assuming a polynomial dependence of the free functions on the poloidal flux, the equilibrium problem is reduced to an inhomogeneous Helmholtz partial differential equation (PDE) subject to homogeneous Dirichlet conditions. An application of the Green's function method leads to a closed form for the circular solution and to a series solution in terms of Mathieu functions for the elliptical case, which is valid for arbitrary elongations. In order to extend the elliptical solution to a D-shaped domain, a boundary perturbation in terms of the triangularity is used. A comparison with the code FLOW [L. Guazzotto et al., Phys. Plasmas 11(2), 604-614 (2004)] is presented for relevant scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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

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↗