Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary value problem”

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 343 records · Page 19

On New Structures in the Theory of Fully Nonlinear Equations

We describe the current state of the theory of equations with m-Hessian stationary and evolution operators. It is quite important that new algebraic and geometric notions appear in this theory. In the present work, a list of those notions is provided. Among them, the notion of m-positivity of matrices is quite important; we provide a proof of an analog of Sylvester’s criterion for such matrices. From this criterion, we easily obtain necessary and sufficient conditions for existence of classical solutions of the first initial boundary-value problem for m-Hessian evolution equations. The asymptotic behavior of m-Hessian evolutions in a semibounded cylinder is considered as well.

97 MATHEMATICS AND COMPUTING↗

Meso-Scale Approximations of Fields Around Clusters of Defects

We present a review of recent results on asymptotic approximations of solutions to boundary value problems in domains with large clusters of small defects. There are no assumptions which require periodicity within the cluster. The asymptotic approximations we discuss are uniform and include the boundary layers occurring in neighborhoods of singularly perturbed boundaries of the domains concerned. The term “meso-scale” is used to describe these approximations since they go beyond the conventional constraints of the homogenization theory.

97 MATHEMATICS AND COMPUTING↗

Enforcing global constraints for the dispersion closure problem: τ 2 -SIMPLE algorithm

Permeability and effective dispersion tensors are critical parameters to characterize flow and transport in porous media at the continuum scale. Homogenization theory defines a framework in which such effective properties are first computed from solving a closure problem in a repeating unit cell of the periodic microstructure and then used in a macroscopic formulation for efficient computation. The closure problem is formulated as a local boundary value problem subjected to global constraints, which guarantee the uniqueness of the solution and can be difficult to satisfy for complex geometries and at high flow conditions. These constraints also ensure that pore-scale pressure, velocity, and concentration fields can be accurately reconstructed from the closure variable. Building on a previous work, here we present a framework that allows to satisfy global constraints associated to both the permeability and the dispersion closure problems by introducing two artificial time scales. The algorithm, called τ 2 -SIMPLE, computes both permeability and effective dispersion given an arbitrarily complex geometry and flow condition. Furthermore, this algorithm is demonstrated to be accurate for both 2D and 3D geometries across varying flow conditions, and thus it can be used to quickly characterize effective properties from porous media images in many applications.

97 MATHEMATICS AND COMPUTING↗

Application of the inverse fast multipole method as a preconditioner in a 3D Helmholtz boundary element method

Here we investigate an efficient preconditioning of iterative methods (such as GMRES) for solving dense linear systems Ax = b that follow from a boundary element method (BEM) for the 3D Helmholtz equation, focusing on the low-frequency regime. While matrix–vector products in GMRES can be accelerated through the low-frequency fast multipole method (LFFMM), the BEM often remains computationally expensive due to the large number of GMRES iterations. We propose the application of the inverse fast multipole method (IFMM) as a preconditioner to accelerate the convergence of GMRES. The IFMM is in essence an approximate direct solver that uses a multilevel hierarchical decomposition and low-rank approximations. The proposed IFMM-based preconditioning has a tunable parameter ε that balances the cost to construct a preconditioner M, which is an approximation of A -1 , and the cost to perform the iterative process by means of M. Namely, using a small (respectively, large) value of ε takes a long (respectively, short) time to construct M, while the number of iterations can be small (respectively, large). A comprehensive set of numerical examples involving various boundary value problems with complicated geometries and mixed boundary conditions is presented to validate the efficiency of the proposed method. We show that the IFMM preconditioner (with a nearly optimal ε of 10 -2 ) clearly outperforms some common preconditioners for the BEM, achieving 1.2–10.8 times speed-up of the computations, in particular when the scale of the underlying scatterer is about five wavelengths or more. In addition, the IFMM preconditioner is capable of solving complicated problems (in a reasonable amount of time) that BD preconditioner can not.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Investigation of technical problems related to retrieval of uncooperative orbiting objects

Preliminary work on high altitude deployment of payloads was carried out. The problem of optimal transfer of the retrieval package from space shuttle to target was studied. The problem was characterized as a nonlinear boundary value problem with a terminal constraint on approach direction. Methods of identifying spin axis and direction are considered.

Kaplan, M. H.↗

Electromagnetic wave propagation in a magneto-plasma filled coaxial structure. I - Theoretical. II - Experimental

This study is concerned with the problem of electromagnetic wave propagation in a magneto-plasma filled coaxial structure. The problem is formulated using the classical boundary value problem approach. A numerical investigation shows the existence of propagating slow modes, backward modes, a quasi-TEM mode, and waveguide-type modes in a magneto-plasma filled coaxial structure. Dispersion curves for these different modes are presented. Measurements have been made of electromagnetic propagation in a coaxial electrode structure filled with longitudinally magnetized plasma. The annular plasma region had a 9.55 cm outer diameter, a 3.82 cm inner diameter and was approximately 60 cm long. A magnetic field of 300 gauss was employed. Electromagnetic wave frequencies were in the range .5 to 2.4 GHz. The plasma was generated by a continuous glow discharge. The resulting dispersion curves closely follow the predicted curves for the quasi-TEM mode.

Askins, H. W., Jr.↗

Numerical quadrature methods for integrals of singular periodic functions and their application to singular and weakly singular integral equations

High accuracy numerical quadrature methods for integrals of singular periodic functions are proposed. These methods are based on the appropriate Euler-Maclaurin expansions of trapezoidal rule approximations and their extrapolations. They are used to obtain accurate quadrature methods for the solution of singular and weakly singular Fredholm integral equations. Such periodic equations are used in the solution of planar elliptic boundary value problems, elasticity, potential theory, conformal mapping, boundary element methods, free surface flows, etc. The use of the quadrature methods is demonstrated with numerical examples.

Sidi, A.↗

Optimal Angle of Attack Control for Aerocapture Maneuvers

Recent investigations of aerocapture maneuvers have considered direct force control as a promising alternative to bank angle control. By modulating the angle of attack and sideslip angle independently, the longitudinal and lateral channels can be decoupled and thus solved separately. However, the optimal angle of attack profile for minimizing post-atmospheric propellant expenditure was previously unknown. This research applies principles of optimal control theory to numerically solve for the optimal angle of attack profile for control of the longitudinal channel. Direct methods are first employed to discretize the angle of attack profile and show that an optimal solution follows a nearly bang-bang (full lift up then full lift down) structure under a variety of conditions along the flight path angle corridor. To corroborate these results, indirect methods were also applied to scope the Optimal Control Problem (OCP) as a Two Point Boundary Value Problem (2P-BVP) which could then be solved numerically. The indirect optimization results aligned closely with the direct results, providing further evidence that a bang-bang structure provides a good approximation for the optimal angle of attack profile for aerocapture maneuvers using direct force control.

aerocapture↗

Quantum models with energy-dependent potentials solvable in terms of exceptional orthogonal polynomials

We construct energy-dependent potentials for which the Schrödinger equations admit solutions in terms of exceptional orthogonal polynomials. Our method of construction is based on certain point transformations, applied to the equations of exceptional Hermite, Jacobi and Laguerre polynomials. We present several examples of boundary-value problems with energy-dependent potentials that admit a discrete spectrum and the corresponding normalizable solutions in closed form.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Characteristically Closed Domains for First Order Strictly Hyperbolic Systems in the Plane

We consider a first order strictly hyperbolic system of n equations with constant coefficients in a bounded domain. It is assumed that the domain is strictly convex relative to characteristics, so that the projection along each characteristic is an involution having two fixed singular points. The natural statement of boundary value problems for such systems requires that singular points go to singular points under such transformations. We present a necessary and sufficient condition for the existence of such domains, called characteristically closed.

97 MATHEMATICS AND COMPUTING↗

Structural zooming research and development of an interactive computer graphical interface for stress analysis of cracks

Engineering problems sometimes involve the numerical solution of boundary value problems over domains containing geometric feature with widely varying scales. Often, a detailed solution is required at one or more of these features. Small details in large structures may have profound effects upon global performance. Conversely, large-scale conditions may effect local performance. Many man-hours and CPU-hours are currently spent in modeling such problems. With the structural zooming technique, it is now possible to design an integrated program which allows the analyst to interactively focus upon a small region of interest, to modify the local geometry, and then to obtain highly accurate responses in that region which reflect both the properties of the overall structure and the local detail. A boundary integral equation analysis program, called BOAST, was recently developed for the stress analysis of cracks. This program can accurately analyze two-dimensional linear elastic fracture mechanics problems with far less computational effort than existing finite element codes. An interactive computer graphical interface to BOAST was written. The graphical interface would have several requirements: it would be menu-driven, with mouse input; all aspects of input would be entered graphically; the results of a BOAST analysis would be displayed pictorially but also the user would be able to probe interactively to get numerical values of displacement and stress at desired locations within the analysis domain; the entire procedure would be integrated into a single, easy to use package; and it would be written using calls to the graphic package called HOOPS. The program is nearing completion. All of the preprocessing features are working satisfactorily and were debugged. The postprocessing features are under development, and rudimentary postprocessing should be available by the end of the summer. The program was developed and run on a VAX workstation, and must be ported to the SUN workstation. This activity is currently underway.

Gerstle, Walter↗

Numerical solution of the stochastic parabolic equation with the dependent operator coefficient

In the present paper, a single step implicit difference scheme for the numerical solution of the stochastic parabolic equation with the dependent operator coefficient is presented. Theorem on convergence estimates for the solution of this difference scheme is established. In applications, this abstract result permits us to obtain the convergence estimates for the solution of difference schemes for the numerical solution of initial boundary value problems for parabolic equations. The theoretical statements for the solution of this difference scheme are supported by the results of numerical experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Image-Driven Hybrid Structural Analysis Based on Continuum Point Cloud Method with Boundary Capturing Technique

Conventional approaches for the structural health monitoring of infrastructures often rely on physical sensors or targets attached to structural members, which require considerable preparation, maintenance, and operational effort, including continuous on-site adjustments. This paper presents an image-driven hybrid structural analysis technique that combines digital image processing (DIP) and regression analysis with a continuum point cloud method (CPCM) built on a particle-based strong formulation. Polynomial regressions capture the boundary shape change due to the structural loading and precisely identify the edge and corner coordinates of the deformed structure. The captured edge profiles are transformed into essential boundary conditions. This allows the construction of a strongly formulated boundary value problem (BVP), classified as the Dirichlet problem. Capturing boundary conditions from the digital image is novel, although a similar approach was applied to the point cloud data. It was shown that the CPCM is more efficient in this hybrid simulation framework than the weak-form-based numerical schemes. Unlike the finite element method (FEM), it can avoid aligning boundary nodes with regression points. A three-point bending test of a rubber beam was simulated to validate the developed technique. The simulation results were benchmarked against numerical results by ANSYS and various relevant numerical schemes. The technique can effectively solve the Dirichlet-type BVP, yielding accurate deformation, stress, and strain values across the entire problem domain when employing a linear strain model and increasing the number of CPCM nodes. In addition, comparative analysis with conventional displacement tracking techniques verifies the developed technique’s robustness. The proposed technique effectively circumvents the inherent limitations of traditional monitoring methods resulting from the reliance on physical gauges or target markers so that a robust and non-contact solution for remote structural health monitoring in real-scale infrastructures can be provided, even in unfavorable experimental environments.

Chemistry↗

Wave-driven torques to drive current and rotation

In the classic Landau damping initial value problem, where a planar electrostatic wave transfers energy and momentum to resonant electrons, a recoil reaction occurs in the nonresonant particles to ensure momentum conservation. To explain how net current can be driven in spite of this conservation, the literature often appeals to mechanisms that transfer this nonresonant recoil momentum to ions, which carry negligible current. However, this explanation does not allow the transport of net charge across magnetic field lines, precluding E X B rotation drive. In this study, we show that in steady state, this picture of current drive is incomplete. Using a simple Fresnel model of the plasma, we show that for lower hybrid waves, the electromagnetic energy flux (Poynting vector) and momentum flux (Maxwell stress tensor) associated with the evanescent vacuum wave become the Minkowski energy flux and momentum flux in the plasma and are ultimately transferred to resonant particles. Thus, the torque delivered to the resonant particles is ultimately supplied by the electromagnetic torque from the antenna, allowing the nonresonant recoil response to vanish and rotation to be driven. We present a warm fluid model that explains how this momentum conservation works out locally, via a Reynolds stress that does not appear in the one-dimensional initial value problem. This model is the simplest that can capture both the nonresonant recoil reaction in the initial-value problem, and the absence of a nonresonant recoil in the steady-state boundary value problem, thus forbidding rotation drive in the former while allowing it in the latter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Approximation of Two-Dimensional Viscoelastic Flows of General Form

We consider the initial-boundary value problem for approximations of the system of integro-differential equations generalizing the equations of motion for viscoelastic fluids. We prove the existence and convergence theorems and give some examples of non-Newtonian fluids described by the model under consideration.

97 MATHEMATICS AND COMPUTING↗