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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
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 OSTI2019

Analysis of the Spectral Stability of the Generalized Runge–Kutta Methods Applied to Initial-Boundary-Value Problems for Equations of the Parabolic Type. II. Implicit Methods

We consider specific realizations of different implicit generalized Runge–Kutta methods as applied to the numerical integration with respect to time of initial-boundary-value problems for the second-order parabolic equations and investigate their spectral stability. It is shown that all implicit generalized Runge–Kutta methods are unconditionally spectrally stable but some of them have the conditional property of monotonicity of the numerical solution with respect to time. The functions of spectral stability of the implicit generalized Runge–Kutta methods are rational. We compare the analytic solution of the nonstationary one-dimensional problem of heat conduction with the numerical solutions of this problem obtained by different implicit generalized Runge–Kutta methods. It is shown that, in this case, the application of the one-stage Radau methods with subsequent discretization of the problem with respect to the space variable leads to the classical forward finite difference scheme (Laasonen scheme), whereas the use of the one-stage Gauss–Legendre method leads to a six-point symmetric scheme (Crank–Nicolson scheme). It is shown that diagonally implicit generalized Nørsett and Burrage methods are realized in almost the same way as the one-stage Radau and Gauss–Legendre methods but their accuracy in the time step is 10–1000 times higher. On the basis of comparison of the numerical and analytic solutions, we conclude that, in order to get practically suitable numerical solutions without any restrictions on the time step, it is reasonable to use one- and three-stage generalized Radau methods or two- and four-stage Lobatto IIIC methods. All other explicit and implicit generalized Runge–Kutta methods require certain restrictions imposed on the time step.

42 ENGINEERING↗
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

Initial-Boundary Value Problem for the Nonstationary Radiative Transfer Equation with Diffuse Reflection and Refraction Conditions

We consider the initial-boundary value problem for the nonstationary radiative transfer equation in a system of semitransparent bodies with the conditions of diffuse reflection and refraction of radiation. We establish the unique solvability of the problem with boundary and initial data in the complete scale of Lebesgue spaces. We obtain estimates for solutions.

42 ENGINEERING↗
DOE OSTI2018

Smoothness in the Dini Space of a Single Layer Potential for a Parabolic System in the Plane

We prove estimates in the Dini space for the second order spatial derivative of a vector parabolic single layer potential. We obtain estimates for higher order derivatives of solutions to initial-boundary-value problems for one-dimensional parabolic systems with Dini-continuous coefficients in a curvilinear domain with nonsmooth lateral boundary.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Solution of A. Einstein’s Problem on the Density of Matter in the Universe

On the basis of the special theory of relativity and a solution of the Schwarzschild equations of general relativity, a generalized equation of Friedmann–Einstein–de Sitter type has been obtained for the expanding Universe. As a result of our work on the boundary value problem for this equation, we have found the eigenvalue of the resulting eigenvalue problem, linking the density of the Universe in the comoving coordinates, the Hubble constant, and the gravitational constant. We have calculated a value for the density of matter in the Universe which is in agreement with the observational data. We have shown that the missing dark energy is the relativistic kinematic and potential energy of the expanding Universe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
NASA NTRS2009

Investigation of Boundary Conditions for Flexible Multibody Spacecraft Dynamics

In support of both the Space Shuttle and International Space Station programs, a set of generic multibody dynamics algorithms integrated within the Trick simulation environment have addressed the variety of on-orbit manipulator simulation requirements for engineering analysis, procedures development and crew familiarization/training at the NASA Johnson Space Center (JSC). Enhancements to these dynamics algorithms are now being driven by a new set of Constellation program requirements for flexible multibody spacecraft simulation. One particular issue that has been discussed within the NASA community is the assumption of cantilever-type flexible body boundary conditions. This assumption has been commonly utilized within manipulator multibody dynamics formulations as it simplifies the computation of relative motion for articulated flexible topologies. Moreover, its use for modeling of space-based manipulators such as the Shuttle Remote Manipulator System (SRMS) and Space Station Remote Manipulator System (SSRMS) has been extensively validated against flight data. For more general flexible spacecraft applications, however, the assumption of cantilever-type boundary conditions may not be sufficient. This paper describes the boundary condition assumptions that were used in the original formulation, demonstrates that this formulation can be augmented to accommodate systems in which the assumption of cantilever boundary conditions no longer applies, and verifies the approach through comparison with an independent model previously validated against experimental hardware test data from a spacecraft flexible dynamics emulator.

Boundary-value problems↗
NASA NTRS1990

Computing Deformations Of Rubbery Materials

Better use made of experimental data in finite-element computations. New formulation of constitutive equations of rubbery, nonlinearly elastic material enables finite-element analysis of boundary-value stress-and-strain problems involving arbitrary shapes and loads. In development of formulation, principal stretches used as arguments of strain-energy-density function.

Peng, Steven T. J.↗
NASA NTRS1990

Solving time-dependent two-dimensional eddy current problems

Transient eddy current calculations are presented for an EM wave-scattering and field-penetrating case in which a two-dimensional transverse magnetic field is incident on a good (i.e., not perfect) and infinitely long conductor. The problem thus posed is of initial boundary-value interface type, where the boundary of the conductor constitutes the interface. A potential function is used for time-domain modeling of the situation, and finite difference-time domain techniques are used to march the potential function explicitly in time. Attention is given to the case of LF radiation conditions.

Lee, Min Eig↗
NASA NTRS1989

Singular trajectories for time-optimal half-loop maneuvers of a high alpha fighter aircraft

Consideration is given to the problem of deriving a time-optimal open-loop control for the half-loop maneuver of a high-alpha aircraft, with initial conditions Mach 0.6 and 15,000 feet. Pontriagin's maximum principle is used to derive candidate optimal solutions. Using the two-point boundary-value algorithm, the flight path angle is maximized for various increasing specified final times until a final time of 13.6 sec yields a 180-deg flight-path angle. As the final time increased from 0.0 to 13.6 sec, the optimization process revealed 13 distinct switching structures of the control law, of which 11 contained singular arcs, and two had double singular arcs.

Hoffman, Eric↗
NASA NTRS1990

Spurious frequencies as a result of numerical boundary treatments

The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.

Abarbanel, Saul↗
DOE OSTI2019

Combination of the Laguerre Transform with the Boundary-element Method for the Solution of Integral Equations with Retarded Kernel

We apply the Laguerre transform with respect to time to a time-dependent boundary-value integral equation encountered in the solution of three-dimensional Dirichlet initial-boundary-value problems for the homogeneous wave equation with homogeneous initial conditions by using the retarded potential of single layer. The obtained system of boundary integral equations is reduced to a sequence of Fredholm integral equations of the first kind that differ solely by the recursively dependent right-hand sides. To find their numerical solution, we use the boundary-element method. We establish an asymptotic estimate of the error of numerical solution and present the results of numerical simulations aimed at finding the solutions of retarded-potential integral equations for model examples.

97 MATHEMATICS AND COMPUTING↗
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.↗
NASA NTRS1992

Computational aspects of helicopter trim analysis and damping levels from Floquet theory

Helicopter trim settings of periodic initial state and control inputs are investigated for convergence of Newton iteration in computing the settings sequentially and in parallel. The trim analysis uses a shooting method and a weak version of two temporal finite element methods with displacement formulation and with mixed formulation of displacements and momenta. These three methods broadly represent two main approaches of trim analysis: adaptation of initial-value and finite element boundary-value codes to periodic boundary conditions, particularly for unstable and marginally stable systems. In each method, both the sequential and in-parallel schemes are used and the resulting nonlinear algebraic equations are solved by damped Newton iteration with an optimally selected damping parameter. The impact of damped Newton iteration, including earlier-observed divergence problems in trim analysis, is demonstrated by the maximum condition number of the Jacobian matrices of the iterative scheme and by virtual elimination of divergence. The advantages of the in-parallel scheme over the conventional sequential scheme are also demonstrated.

Gaonkar, Gopal H.↗