Search NASA⌕ Search

SEARCH · Search NASA

Results for “Partial integral differential equations”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Two-dimensional integrating matrices on rectangular grids

The use of integrating matrices in solving differential equations associated with rotating beam configurations is examined. In vibration problems, by expressing the equations of motion of the beam in matrix notation, utilizing the integrating matrix as an operator, and applying the boundary conditions, the spatial dependence is removed from the governing partial differential equations and the resulting ordinary differential equations can be cast into standard eigenvalue form. Integrating matrices are derived based on two dimensional rectangular grids with arbitrary grid spacings allowed in one direction. The derivation of higher dimensional integrating matrices is the initial step in the generalization of the integrating matrix methodology to vibration and stability problems involving plates and shells.

Lakin, W. D.↗
NASA NTRS1986

Optimal control of aeroassisted plane change maneuver using feedback expansions

A guidance law for an aeroassisted plane change maneuver is developed by an asymptotic expansion technique using a small parameter which essentially represents the ratio of the inertial forces to the atmospheric forces. This guidance law minimizes the energy loss while meeting terminal constraints on the altitude, flight path angle, and heading angle. By neglecting the inertial forces, the resulting optimization problem is integrable and can be determined in closed form. This zeroth-order solution is the first term in an asymptotic series solution of the Hamilton-Jacobi-Bellman equation. The remaining terms are determined from the solution of a first-order, linear partial differential equation whose solution requires only quadrature integration. Our initial results in using this guidance scheme are encouraging.

Mishne, D.↗
NASA NTRS1982

Recent advances in electromagnetic synthesis and analysis of dual-shaped reflector antennas

Dual-shaped reflectors have been used for many years. Thus, these reflectors have been used as high gain antennas on Voyagers 1 and 2. The objectives of the geometrical optics (GO) dual shaped synthesis are considered. Concerning the synthesis of dual shaped reflectors, it has been shown for circular symmetric reflectors that an exact GO solution can be found to the problem of transforming, by two reflections, any feed pattern into any aperture distribution. This problem involves solving two simultaneous nonlinear ordinary differential equations. The same approach for offset geometry leads to two simultaneous nonlinear partial differential equations. It is shown that these equations could also be integrated numerically, except that in general these equations are not total and therefore, in general, they do not have a 'smooth' solution. It is further shown that the offset partial differentials often very nearly form a total differential in many cases of practical importance.

Galindo-Israel, V.↗
NASA NTRS1977

Derivation of an integral equation for transonic flows

It is noted that the nonlinear partial differential equation for the perturbation velocity potential and boundary conditions describing steady inviscid compressible transonic flow past a thin two-dimensional airfoil can be transformed into a singular integrodifferential equation and that differentiation of the latter yields an integral equation. Two forms of this integral equation currently exist: one for the singularity that is enclosed in an infinitely long strip of vanishing thickness and the other for the singularity that is enclosed in a vanishing circle. In the present article, a more general integral equation is derived by enclosing the singularity in a vanishing rectangular cavity of arbitrary aspect ratio. The two existing forms of this equation are deduced as special cases distinguished by the respective values for the aspect ratio (infinity for the first form and unity for the second).

Ogana, W.↗
NASA NTRS1947

Voltera's Solution of the Wave Equation as Applied to Three-Dimensional Supersonic Airfoil Problems

A surface integral is developed which yields solutions of the linearized partial differential equation for supersonic flow. These solutions satisfy boundary conditions arising in wing theory. Particular applications of this general method are made, using acceleration potentials, to flat surfaces and to uniformly loaded lifting surfaces. Rectangular and trapezoidal plan forms are considered along with triangular forms adaptable to swept-forward and swept-back wings. The case of the triangular plan form in sideslip is also included. Emphasis is placed on the systematic application of the method to the lifting surfaces considered and on the possibility of further application.

Heslet, Max A↗
NASA NTRS1978

Remarks on boundary-fitted coordinate system generation

The essential part of numerical solutions of partial differential equations is the representation of gradients and integrals by, respectively, differences between points and summations over points. In order for such numerical representations to be accurate, it is necessary that these points be more closely spaced in regions of large gradients. The need for accurate representation is particularly acute near body surfaces, since the boundary conditions are generally the most influential part of a partial differential equation solution. This is especially true of viscous solutions at high Reynolds number, where very large gradients occur in the boundary layer.

Thompson, J. F.↗
NASA NTRS1979

A simplified computer program for the prediction of the linear stability behavior of liquid propellant combustors

A program for predicting the linear stability of liquid propellant rocket engines is presented. The underlying model assumptions and analytical steps necessary for understanding the program and its input and output are also given. The rocket engine is modeled as a right circular cylinder with an injector with a concentrated combustion zone, a nozzle, finite mean flow, and an acoustic admittance, or the sensitive time lag theory. The resulting partial differential equations are combined into two governing integral equations by the use of the Green's function method. These equations are solved using a successive approximation technique for the small amplitude (linear) case. The computational method used as well as the various user options available are discussed. Finally, a flow diagram, sample input and output for a typical application and a complete program listing for program MODULE are presented.

Mitchell, C. E.↗
NASA NTRS1979

Theoretical determination of nozzle admittances using a finite element method

A finite element method (FEM) is used to predict the admittances of axisymmetric nozzles. The flow in the nozzle is assumed to be isentropic and the disturbances are assumed to be small so that linear analyses apply. An approximate two dimensional compressible flow model is used to describe the steady flow in the nozzle. The propagation of acoustic disturbances is governed by the complete linear acoustic wave equation. This partial differential wave equation is transformed to an integral equation using Galerkin's method and Green's theorem is applied so that the acoustic boundary conditions can be introduced through the boundary residuals. A two dimensional finite element method using linear triangular elements is used to solve the integral acoustic equation. A one dimensional FEM is used to solve the reduced nozzle acoustic equation developed by Crocco and the solution is used to verify the sufficiency of the boundary residual formation. It is shown that agreement between predicted values of the admittance and experimental data is quite good.

Sigman, R. K.↗
NASA NTRS1989

Semidiscrete Galerkin modelling of compressible viscous flow past a circular cone at incidence

A numerical study of the laminar and compressible boundary layer, about a circular cone in a supersonic free stream, is presented. It is thought that if accurate and efficient numerical schemes can be produced to solve the boundary layer equations, they can be joined to numerical codes that solve the inviscid outer flow. The combination of these numerical codes is competitive with the accurate, but computationally expensive, Navier-Stokes schemes. The primary goal is to develop a finite element method for the calculation of 3-D compressible laminar boundary layer about a yawed cone. The proposed method can, in principle, be extended to apply to the 3-D boundary layer of pointed bodies of arbitrary cross section. The 3-D boundary layer equations governing supersonic free stream flow about a cone are examined. The 3-D partial differential equations are reduced to 2-D integral equations by applying the Howarth, Mangler, Crocco transformations, a linear relation between viscosity, and a Blasius-type of similarity variable. This is equivalent to a Dorodnitsyn-type formulation. The reduced equations are independent of density and curvature effects, and resemble the weak form of the 2-D incompressible boundary layer equations in Cartesian coordinates. In addition the coordinate normal to the wall has been stretched, which reduces the gradients across the layer and provides high resolution near the surface. Utilizing the parabolic nature of the boundary layer equations, a finite element method is applied to the Dorodnitsyn formulation. The formulation is presented in a Petrov-Galerkin finite element form and discretized across the layer using linear interpolation functions. The finite element discretization yields a system of ordinary differential equations in the circumferential direction. The circumferential derivatives are solved by an implicit and noniterative finite difference marching scheme. Solutions are presented for a 15 deg half angle cone at angles of attack of 5 and 10 deg. The numerical solutions assume a laminar boundary layer with free stream Mach number of 7. Results include circumferential distribution of skin friction and surface heat transfer, and cross flow velocity distributions across the layer.

Meade, Andrew James, Jr.↗
DOE OSTI2025

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Assessment of a transitional boundary layer theory at low hypersonic Mach numbers

An investigation was carried out to assess the accuracy of a transitional boundary layer theory in the low hypersonic Mach number regime. The theory is based upon the simultaneous numerical solution of the boundary layer partial differential equations for the mean motion and an integral form of the turbulence kinetic energy equation which controls the magnitude and development of the Reynolds stress. Comparisions with experimental data show the theory is capable of accurately predicting heat transfer and velocity profiles through the transitional regime and correctly predicts the effects of Mach number and wall cooling on transition Reynolds number. The procedure shows promise of predicting the initiation of transition for given free stream disturbance levels. The effects on transition predictions of the pressure dilitation term and of direct absorption of acoustic energy by the boundary layer were evaluated.

Shamroth, S. J.↗
NASA NTRS1985

A Survey of Composite Grid Generation for General Three-dimensional Sections

The generation and use of composite grids for general three-dimensional physical boundary configurations is discussed, and the availability of several codes or procedures is noted. With the composite framework, the physical region is segmented into sub-regions, each bounded by six curved sides, and a grid is generated in each sub-region. These grids may be joined at the interfaces between the sub-regions with various degrees of continuity. This structure allows codes to be constructed to operate on rectangular blocks in computational space, so that existing solution procedures can be readily incorporated in the construction of codes for general configurations. Numerical grid generation is an integral part of the numerical solution of partial differential equations and is one of the pacing items in the development of codes for general configurations. The numerically generated grid frees the computational simulation from restriction to certain boundary shapes and allows general codes to be written in which the boundary shape is specified simply by input. The numerically generated grid allows all computation to be done on a fixed square grid in the computational space, which is always rectangular by construction.

Thompson, J. F.↗
NASA NTRS2012

Assessing Tsunami Vulnerabilities of Geographies with Shallow Water Equations

Tsunami preparedness is crucial for saving human lives in case of disasters that involve massive water movement. In this work, we develop a framework for visual assessment of tsunami preparedness of geographies. Shallow water equations (also called Saint Venant equations) are a set of hyperbolic partial differential equations that are derived by depth-integrating the Navier-Stokes equations and provide a great abstraction of water masses that have lower depths compared to their free surface area. Our specific contribution in this study is to use Microsoft's XNA Game Studio to import underwater and shore line geographies, create different tsunami scenarios, and visualize the propagation of the waves and their impact on the shore line geography. Most importantly, we utilized the computational power of graphical processing units (GPUs) as HLSL based shader files and delegated all of the heavy computations to the GPU. Finally, we also conducted a validation study, in which we have tested our model against a controlled shallow water experiment. We believe that such a framework with an easy to use interface that is based on readily available software libraries, which are widely available and easily distributable, would encourage not only researchers, but also educators to showcase ideas.

Aras, Rifat↗
DOE OSTI2026

Parallel-in-Time Solution of Allen-Cahn Equations by Integrating Operator Learning into the Parareal Method

While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Examples Of Synthesis Of Dual-Shaped Reflectors

Report presents examples to demonstrate validity and utility of method of synthesis of offset dual-shape reflectors. Method of synthesis described by the authors in previous journal article. Current report reviews derivation of partial differential equations and iterative method of numerical solution. Discusses significance of starting point of numerical integration on each reflector surface; this point could be at center, on outer rim, or at interior point. Emphasizes that one of notable attributes of partial differential equations is speed with which they can be solved.

Galindo, Victor↗
NASA NTRS1982

Integrating matrix formulations for vibrations of rotating beams including the effects of concentrated masses

By expressing partial differential equations of motion in matrix notation, utilizing the integrating matrix as a spatial operator, and applying the boundary conditions, the resulting ordinary differential equations can be cast into standard eigenvalue form upon assumption of the usual time dependence. As originally developed, the technique was limited to beams having continuous mass and stiffness properties along their lengths. Integrating matrix methods are extended to treat the differential equations governing the flap, lag, or axial vibrations of rotating beams having concentrated masses. Inclusion of concentrated masses is shown to lead to the same kind of standard eigenvalue problem as before, but with slightly modified matrices.

Lakin, W. D.↗