Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential algebraic 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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Development of the general interpolants method for the CYBER 200 series of supercomputers

The General Interpolants Method (GIM) is a 3-D, time-dependent, hybrid procedure for generating numerical analogs of the conservation laws. This study is directed toward the development and application of the GIM computer code for fluid dynamic research applications as implemented for the Cyber 200 series of supercomputers. An elliptic and quasi-parabolic version of the GIM code are discussed. Turbulence models, algebraic and differential equations, were added to the basic viscous code. An equilibrium reacting chemistry model and an implicit finite difference scheme are also included.

Stalnaker, J. F.↗
NASA NTRS2002

Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation

The ice accrued on airfoils can have enormously complicated shapes with multiple protruded horns and feathers. In this paper, several blocking topologies are proposed and evaluated on their ability to produce high-quality structured multi-block grid systems. A transition layer grid is introduced to ensure that jaggedness on the ice-surface geometry do not to propagate into the domain. This is important for grid-generation methods based on hyperbolic PDEs (Partial Differential Equations) and algebraic transfinite interpolation. A 'thick' wrap-around grid is introduced to ensure that grid lines clustered next to solid walls do not propagate as streaks of tightly packed grid lines into the interior of the domain along block boundaries. For ice shapes that are not too complicated, a method is presented for generating high-quality single-block grids. To demonstrate the usefulness of the methods developed, grids and CFD solutions were generated for two iced airfoils: the NLF0414 airfoil with and without the 623-ice shape and the B575/767 airfoil with and without the 145m-ice shape. To validate the computations, the computed lift coefficients as a function of angle of attack were compared with available experimental data. The ice shapes and the blocking topologies were prepared by NASA Glenn's SmaggIce software. The grid systems were generated by using a four-boundary method based on Hermite interpolation with controls on clustering, orthogonality next to walls, and C continuity across block boundaries. The flow was modeled by the ensemble-averaged compressible Navier-Stokes equations, closed by the shear-stress transport turbulence model in which the integration is to the wall. All solutions were generated by using the NPARC WIND code.

Chi, X.↗
NASA NTRS1974

Algebraic Riccati equations in zero-sum differential games

The procedure for finding the closed-loop Nash equilibrium solution of two-player zero-sum linear time-invariant differential games with quadratic performance criteria and classical information pattern may be reduced in most cases to the solution of an algebraic Riccati equation. Based on the results obtained by Willems, necessary and sufficient conditions for existence of solutions to these equations are derived, and explicit conditions for a scalar example are given.

Johnson, T. L.↗
DOE OSTI2022

Simulation of crumpled sheets via alternating quasistatic and dynamic representations

In this work, we present a method for simulating the large-scale deformation and crumpling of thin, elastoplastic sheets. Motivated by the physical behavior of thin sheets during crumpling, two different formulations of the governing equations of motion are used: (1) a quasistatic formulation that effectively describes smooth deformations, and (2) a fully dynamic formulation that captures large changes in the sheet's velocity. The former is a differential-algebraic system of equations integrated implicitly in time, while the latter is a set of ordinary differential equations (ODEs) integrated explicitly. We adopt a hybrid integration scheme to adaptively alternate between the quasistatic and dynamic representations as appropriate. Further, we demonstrate the capacity of this method to effectively simulate a variety of crumpling phenomena. Finally, we show that statistical properties, notably the accumulation of creases under repeated loading, as well as the area distribution of facets, are consistent with experimental observations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Comparison of level-1 and level-2 buckling and postbuckling solutions

The accuracy of the solutions obtained with the Level-1 computational modules of DISDECO (Delft Interactive Shell Design Code) are investigated by comparing them with solutions obtained with ANILISA, a Level-2 computational module. Whereas the Level-1 solutions are based on the use of trigonometric functions to represent the out-of-plane displacement W and involve Galerkin type approximations to reduce the governing partial differential equations to simple algebraic expressions, the Level-2 solutions employ a truncated Fourier expansion in the circumferential direction and solve the resulting ordinary differential equations forming a nonlinear eigenvalue problem by a Stodola like technique. This involves the numerical integration of the set of ordinary differential equations, which in turn makes it possible to satisfy the specified boundary conditions exactly. It is shown that in many applications the Level-1 predictions concerning the critical buckling load of the perfect structure and its imperfection sensitivity are quite accurate and reliable. For those cases where this is not true, it is shown how the pre and postprocessing facilities of DISDECO can be used efficiently to arrive at a reliable prediction. This involves always the use of the Level-2 computational module ANILISA with its rigorous prebuckling solution and its feature which makes it possible to satisfy the specified boundary conditions exactly.

Arbocz, J.↗
NASA NTRS1978

An extension of A-stability to alternating direction implicit methods

An alternating direction implicit (ADI) scheme was constructed by the method of approximate factorization. An A-stable linear multistep method (LMM) was used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation. Sufficient conditions for the A-stability of the LMM were determined by applying the theory of positive real functions to reduce the stability analysis of the partial differential equations to a simple algebraic test. A linear test equation for partial differential equations is defined and then used to analyze the stability of approximate factorization schemes. An ADI method for the three-dimensional heat equation is also presented.

Warming, R. F.↗
NASA NTRS1983

A unified mathematical theory of electrophoretic processes

A mathematical theory is presented which shows that each of the four classical electrophoretic modes (zone electrophoresis, moving boundary electrophoresis, isotachophoresis, and isoelectric focusing) is based on the same general principles and can collectively be described in terms of a single set of equations. This model can predict the evolution of the four electrophoretic modes as a function of time. The model system is one-dimensional, neglecting the effects of electroosmosis, temperature gradients, and any bulk flows of liquid. The model is based on equations which express the components' dissociation equilibria, the mass transport due to electromigration and diffusion, electroneutrality, and the conservation of mass and charge. The model consists of a system of coupled partial differential and nonlinear algebraic equations which can be solved numerically by use of a computer. The versatility of this model was verified using an example of a three-component system containing cacodylate, tris hydroxylmethylaminomethane, and histidine. Results show that this model not only correctly predicts the characteristic features of each electrophoretic mode, but also gives details of the concentration, pH, and conductivity profiles not easily amenable to direct experimental measurement.

Bier, M.↗
DOE OSTI2023

Parallel Algebraic Multigrid for Fusion and Higher-Order PDEs

Multigrid methods play a key role in large-scale scientific simulation because they are among the fastest and most scalable approaches for solving the underlying sparse linear systems of equations that arise from a wide array of Partial Differential Equation (PDE) discretizations. Algebraic multigrid (AMG) is a special type of multigrid method that depends only on the description of the linear system, giving it better portability and broader applicability than geometric multigrid, as it requires no explicit knowledge of the problem geometry. Even though these methods are widely used today, there are still applications where further development is needed. In this report, we focus on PDEs with higher-order terms (e.g., fourth order), concentrating on a PDE that arises in tokamak edge plasma simulations (a tokamak is a machine that confines a plasma using magnetic fields and is believed to be the leading plasma confinement concept for future fusion power plants). General multigrid relaxes a linear system on coarser grids and reverses this process with interpolation, but standard AMG methods struggle with the aforementioned higher-order PDEs. We investigate cyclic coarsening and interpolation heuristics, as well as new iterative approximation methods of refining the solution at each grid to improve the existing multigrid approach. To this end, we ensure that these techniques are transferable to a parallelized setting with LLNL’s supercomputers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A gradient-based deep neural network model for simulating multiphase flow in porous media

We report simulation of multiphase flow in porous media is crucial for the effective management of subsurface energy and environment-related activities. The numerical simulators used for modeling such processes rely on spatial and temporal discretization of the governing mass and energy balance partial-differential equations (PDEs) into algebraic systems via finite-difference/volume/element methods. These simulators usually require dedicated software development and maintenance, and suffer low efficiency from a runtime and memory standpoint for problems with multi-scale heterogeneity, coupled-physics processes or fluids with complex phase behavior. Therefore, developing cost-effective, data-driven models can become a practical choice, and in this work, we choose deep learning approaches as they can handle high dimensional data and accurately predict state variables with strong nonlinearity. In this paper, we describe a gradient-based deep neural network (GDNN) constrained by the physics related to multiphase flow in porous media. We tackle the nonlinearity of flow in porous media induced by rock heterogeneity, fluid properties, and fluid-rock interactions by decomposing the nonlinear PDEs into a dictionary of elementary differential operators. We use a combination of operators to handle rock spatial heterogeneity and fluid flow by advection. Since the augmented differential operators are inherently related to the physics of fluid flow, we treat them as first principles prior knowledge to regularize the GDNN training. We use the example of pressure management at geologic CO 2 storage sites, where CO 2 is injected in saline aquifers and brine is produced, and apply GDNN to construct a predictive model that is trained with physics-based simulation data and emulates the physics process. We demonstrate that GDNN can effectively predict the nonlinear patterns of subsurface responses, including the temporal and spatial evolution of the pressure and saturation plumes. We also successfully extend the GDNN to convolutional neural network (CNN), namely gradient-based CNN (GCNN), and validate its capability to improve the prediction accuracy. GDNN has great potential to tackle challenging problems that are governed by highly nonlinear physics and enable the development of data-driven models with higher fidelity.

42 ENGINEERING↗
NASA NTRS1976

The method of Ritz applied to the equation of Hamilton

Without any reference to the theory of differential equations, the initial value problem of the nonlinear, nonconservative double pendulum system is solved by the application of the method of Ritz to the equation of Hamilton. Also shown is an example of the reduction of the traditional eigenvalue problem of linear, homogeneous, differential equations of motion to the solution of a set of nonhomogeneous algebraic equations. No theory of differential equations is used. Solution of the time-space path of the linear oscillator is demonstrated and compared to the exact solution.

Bailey, C. D.↗
DOE OSTI2025

Parametric matrix models

We present a general class of machine learning algorithms called parametric matrix models. In contrast with most existing machine learning models that imitate the biology of neurons, parametric matrix models use matrix equations that emulate physical systems. Similar to how physics problems are usually solved, parametric matrix models learn the governing equations that lead to the desired outputs. Parametric matrix models can be efficiently trained from empirical data, and the equations may use algebraic, differential, or integral relations. While originally designed for scientific computing, we prove that parametric matrix models are universal function approximators that can be applied to general machine learning problems. After introducing the underlying theory, we apply parametric matrix models to a series of different challenges that show their performance for a wide range of problems. For all the challenges tested here, parametric matrix models produce accurate results within an efficient and interpretable computational framework that allows for input feature extrapolation.

Computational science↗
DOE OSTI2021

pnnl/neural_ODE_ICLR2020

We show how to model discrete ordinary differential equations (ODE) with algebraic nonlinearities as deep neural networks with varying degrees of prior knowledge. We derive the stability guarantees of the network layers based on the implicit constraints imposed on the weight's eigenvalues. Moreover, we show how to use barrier methods to generically handle additional inequality constraints. We demonstrate the prediction accuracy of learned neural ODEs evaluated on open-loop simulations compared to ground truth dynamics with bi-linear terms.

Tuor, Aaron↗
DOE OSTI2023

Four lectures on Euler integrals

These lecture notes provide a self-contained introduction to Euler integrals, which are frequently encountered in applications. In particle physics, they arise as Feynman integrals or string amplitudes. Our four selected topics demonstrate the diverse mathematical techniques involved in the study of Euler integrals, including polyhedral geometry, very affine varieties, differential equations, and computational algebra.

Matsubara-Heo, Saiei-Jaeyeong↗
DOE OSTI2021

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1977

Mathematical modelling of a human external respiratory system

A closed system of algebraic and common differential equations solved by computer is investigated. It includes equations which describe the activity pattern of the respiratory center, the phrenic nerve, the thrust produced by the diaphragm as a function of the lung volume and discharge frequency of the phrenic nerve, as well as certain relations of the lung stretch receptors and chemoreceptors on various lung and blood characteristics, equations for lung biomechanics, pulmonary blood flow, alveolar gas exchange and capillary blood composition equations to determine various air and blood flow and gas exchange parameters, and various gas mixing and arterial and venous blood composition equations, to determine other blood, air and gas mixing characteristics. Data are presented by means of graphs and tables, and some advantages of this model over others are demonstrated by test results.

Source record↗
NASA NTRS1977

Unsteady subsonic and supersonic potential aerodynamics for complex configurations

A recently developed general theory for unsteady compressible potential fluid dynamics for complex-configuration aircraft is reviewed. The method is based on a combination of the following techniques: Green's function method (to transform the differential equation into an integral differential-delay equation), finite element method (to transform the equation into a set of differential-delay equations in time), and the Laplace transform method (to transform the differential-delay equations into algebraic equations).

Morino, L.↗