Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “algebraic methods”

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 307 records · Page 17

Application of boundary integral method to elastoplastic analysis of V-notched beams

The boundary integral equation method was applied in the solution of the plane elastoplastic problem. The use of this method was illustrated by obtaining stress and strain distributions for a number of specimens with a single-edge notch and subjected to pure bending. The boundary integral equation method reduced the inhomogeneous biharmonic equation to two coupled Fredholm-type integral equations. These integral equations were replaced by a system of simultaneous algebraic equations and solved numerically in conjunction with a method of successive elastic solutions.

Rzasnicki, W.↗

Application of boundary integral method to elastoplastic analysis of V-notched beams

The boundary integral equation method was applied in the solution of the plane elastoplastic problems. The use of this method was illustrated by obtaining stress and strain distributions for a number of specimens with a single edge notch and subjected to pure bending. The boundary integral equation method reduced the nonhomogeneous biharmonic equation to two coupled Fredholm-type integral equations. These integral equations were replaced by a system of simultaneous algebraic equations and solved numerically in conjunction with the method of successive elastic solutions.

Rzasnicki, W.↗

Elliptic generation of composite three-dimensional grids about realistic aircraft

An elliptic method for generating composite grids about realistic aircraft is presented. A body-conforming grid is first generated about the entire aircraft by the solution of Poisson's differential equation. This grid has relatively coarse spacing, and it covers the entire physical domain. At boundary surfaces, cell size is controlled and cell skewness is nearly eliminated by inhomogeneous terms, which are found automatically by the program. Certain regions of the grid in which high gradients are expected, and which map into rectangular solids in the computational domain, are then designated for zonal refinement. Spacing in the zonal grids is reduced by adding points with a simple, algebraic scheme. Details of the grid generation method are presented along with results of the present application, a wing-body configuration based on the F-16 fighter aircraft.

Sorenson, R. L.↗

Elliptic surface grid generation on minimal and parametrized surfaces

An elliptic grid generation method, which generates boundary conforming grids in a two dimensional physical space, is presented. The method is based on the composition of an algebraic and elliptic transformation. The composite mapping obeys the Poisson grid generation system with control functions specified by the algebraic transformation. It is shown that the grid generation on a minimal surface in a three dimensional space is equivalent to the grid generation in a two dimensional domain in physical space. A second elliptic grid generation method, which generates boundary conforming grids on smooth surfaces, is presented. Concerning surface modeling, it is shown that bicubic Hermit interpolation is an excellent method to generate a smooth surface crossing a discrete set of control points.

Spekreijse, S. P.↗

Explicit, implicit, and hybrid methods

Time integration methods can be separated into two groups: explicit and implicit. Methods which do not involve the solution of any algebraic equations are called explicit, while those that require the solution of equations are called implicit. The relative advantages and disadvantages of explicit and implicit methods are summarized. The major trend in the past decade of research was to use hybrization methods to take advantage of the complementary nature of the positive attributes of explicit and implicit integration. These trends are briefly discussed.

Belytschko, T.↗

Application of boundary integral method to elastoplastic analysis of v-notched beams

The boundary integral equation method was applied in the solution of the plane elastoplastic problems. The use of this method was illustrated by obtaining stress and strain distributions for a number of specimens with a single edge notch and subjected to pure bending. The boundary integral equation method reduced the non-homogeneous biharmonic equation to two coupled Fredholm-type integral equations. These integral equations were replaced by a system of simultaneous algebraic equations and solved numerically in conjunction with the method of successive elastic solutions.

Rzasnicki, W.↗

A least-squares finite element method for 3D incompressible Navier-Stokes equations

The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations, and results in symmetric, positive definite algebraic system. An additional compatibility equation, i.e., the divergence of vorticity vector should be zero, is included to make the first-order system elliptic. The Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. The flow in a half of 3D cubic cavity is calculated at Re = 100, 400, and 1,000 with 50 x 52 x 25 trilinear elements. The Taylor-Gortler-like vortices are observed at Re = 1,000.

Jiang, Bo-Nan↗

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↗

A Game-Theoretic Quantum Algorithm for Solving Magic Squares

Variational quantum algorithms (VQAs) offer a promising near-term approach to finding optimal quantum strategies for playing non-local games. These games test quantum correlations beyond classical limits and enable entanglement verification. In this work, we present a variational framework for the Magic Square Game (MSG), a two-player non-local game with perfect quantum advantage. We construct a value Hamiltonian that encodes the game’s parity and consistency constraints, then optimize parameterize quantum circuits to minimize this cost. Our approach build on the stabilizer formalism, leverages commutation structure for circuit design, and is hardware-efficient. Compared to existing work, our contribution emphasizes algebraic structure an interpretability. We validate our method through numerical experiments and outline generalizations to larger games.

Chehade, Sarah [ORNL]↗

Quasi-periodic orbits about the translunar libration point.

Analytical solutions for quasi-periodic orbits about the translunar libration point are obtained by using the method of Lindstedt-Poincare and computerized algebraic manipulations. The solutions include the effects of nonlinearities, lunar orbital eccentricity, and the sun's gravitational field. For a small-amplitude orbit, the orbital path as viewed from the earth traces out a Lissajous figure. This is due to a small difference in the fundamental frequencies of the in-plane and out-of-plane oscillations. However, when the amplitude of the in-plane oscillation is greater than 32,379 km, there is a corresponding value of the out-of-plane amplitude that will produce a path where the fundamental frequencies are equal. This synchronized trajectory describes a 'halo orbit' of the moon.

Farquhar, R. W.↗

Quasi-periodic orbits about the translunar libration point.

Analytical solutions for quasi-periodic orbits about the translunar libration point are obtained by using the method of Lindstedt-Poincare and computerized algebraic manipulations. The solutions include the effects of nonlinearities, lunar orbital eccentricity, and the sun's gravitational field. For a small-amplitude orbit, the orbital path as viewed from the earth traces out a Lissajous figure. This is due to a small difference in the fundamental frequencies of the in-plane and out-of-plane oscillations. However, when the amplitude of the in-plane oscillation is greater than 32,379 km, there is a corresponding value of the out-of-plane amplitude that will produce a path where the fundamental frequencies are equal. This synchronized trajectory describes a 'halo orbit' of the moon.

Farquhar, R. W.↗

Initial values for the integration scheme to compute the eigenvalues for propagation in ducts

A scheme for the calculation of eigenvalues in the problem of acoustic propagation in a two-dimensional duct is described. The computation method involves changing the coupled transcendental nonlinear algebraic equations into an initial value problem involving a nonlinear ordinary differential equation. The simplest approach is to use as initial values the hardwall eigenvalues and to integrate away from these values as the admittance varies from zero to its actual value with a linear variation. The approach leads to a powerful root finding routine capable of computing the transverse and axial wave numbers for two-dimensional ducts for any frequency, lining, admittance and Mach number without requiring initial guesses or starting points.

Eversman, W.↗

The axisymmetric elasticity problem for a laminated plate containing a circular hole

The elasticity problem for a laminated thick plate which consists of two bonded dissimilar layers and which contains a circular hole is considered. The problem is formulated for arbitrary axisymmetric tractions on the hole surface by using the Love strain function. Through the expansion of the boundary conditions into Fourier series the problem is reduced to an infinite system of algebraic equations which is solved by the method of reduction. Of particular interest in the problem are the stresses along the interface as they relate to the question of delamination failure of the composite plate. These stresses are calculated and are observed to become unbounded at the hole boundary. An approximate treatment of the singular behavior of the stress state is presented and the stress intensity factors are calculated.

Delale, F.↗

Three-dimensional elasticity solution of an infinite plate with a circular hole

The elasticity problem for a thick plate with a circular hole is formulated in a systematic fashion by using the z-component of the Galerkin vector and that of Muki's harmonic vector function. The problem was originally solved by Alblas. The reasons for reconsidering it are to develop a technique which may be used in solving the elasticity problem for a multilayered plate and to verify and extend the results given by Alblas. The problem is reduced to an infinite system of algebraic equations which is solved by the method of reduction. Various stress components are tabulated as functions of a/h, z/h, r/a, and nu, a and 2h being the radius of the hole and the plate thickness and nu, the Poisson's ratio. The significant effect of the Poisson's ratio on the behavior and the magnitude of the stresses is discussed.

Delale, F.↗

Grid flexibility and patching techniques

The numerical determination of combustor flowfields is of great value to the combustor designer. An a priori knowledge of the flow behavior can speed the combustor design process and reduce the number of experimental test rigs required to arrive at an optimal design. Even 2-D steady incompressible isothermal flow predictions are of use; many codes of this kind are available, each employing different techniques to surmount the difficulties arising from the nonlinearity of the governing equations and from typically irregular combustor geometries. Mapping techniques (algebraic and elliptic PDE), and adaptive grid methods (both multi-grid and grid embedding) as applied to axisymmetric combustors are discussed.

Keith, T. G.↗

A Newton/upwind method and numerical study of shock wave/boundary layer interactions

The objective of the paper is two-fold. First, an upwind/central differencing method for solving the steady Navier-Stokes equations is described. The symmetric line relation method is used to solve the resulting algebraic system to achieve high computational efficiency. The grid spacings used in the calculations are determined from the triple-deck theory, in terms of Mach and Reynolds numbers and other flow parameters. Thus the accuracy of the numerical solutions is improved by comparing them with experimental, analytical, and other computational results. Secondly, the shock wave/boundary layer interactions are studied numerically, with special attention given to the flow separation. The concept of free interaction is confirmed. Although the separated region varies with Mach and Reynolds numbers, it is found that the transverse velocity component behind the incident shock, which has not been identified heretofore, is also an important parameter. A small change of this quantity is sufficient to eliminate the flow separation entirely.

Liou, Meng-Sing↗

Active vibration mitigation of distributed parameter, smart-type structures using Pseudo-Feedback Optimal Control (PFOC)

A new, near-optimal feedback control technique is introduced that is shown to provide excellent vibration attenuation for those distributed parameter systems that are often encountered in the areas of aeroservoelasticity and large space systems. The technique relies on a novel solution methodology for the classical optimal control problem. Specifically, the quadratic regulator control problem for a flexible vibrating structure is first cast in a weak functional form that admits an approximate solution. The necessary conditions (first-order) are then solved via a time finite-element method. The procedure produces a low dimensional, algebraic parameterization of the optimal control problem that provides a rigorous basis for a discrete controller with a first-order like hold output. Simulation has shown that the algorithm can successfully control a wide variety of plant forms including multi-input/multi-output systems and systems exhibiting significant nonlinearities. In order to firmly establish the efficacy of the algorithm, a laboratory control experiment was implemented to provide planar (bending) vibration attenuation of a highly flexible beam (with a first clamped-free mode of approximately 0.5 Hz).

Patten, W. N.↗

An approximate solution for interlaminar stresses in laminated composites: Applied mechanics program

An approximate solution for interlaminar stresses in finite width, laminated composites subjected to uniform extensional, and bending loads is presented. The solution is based upon the principle of minimum complementary energy and an assumed, statically admissible stress state, derived by considering local material mismatch effects and global equilibrium requirements. The stresses in each layer are approximated by polynomial functions of the thickness coordinate, multiplied by combinations of exponential functions of the in-plane coordinate, expressed in terms of fourteen unknown decay parameters. Imposing the stationary condition of the laminate complementary energy with respect to the unknown variables yields a system of fourteen non-linear algebraic equations for the parameters. Newton's method is implemented to solve this system. Once the parameters are known, the stresses can be easily determined at any point in the laminate. Results are presented for through-thickness and interlaminar stress distributions for angle-ply, cross-ply (symmetric and unsymmetric laminates), and quasi-isotropic laminates subjected to uniform extension and bending. It is shown that the solution compares well with existing finite element solutions and represents an improved approximate solution for interlaminar stresses, primarily at interfaces where global equilibrium is satisfied by the in-plane stresses, but large local mismatch in properties requires the presence of interlaminar stresses.

Rose, Cheryl A.↗