Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Sensitivity analysis of a layered piezoelectric system using ZFEM

The complex variable finite element method (ZFEM) is a numerical technique which aims to find the partial derivatives of the independent variables with respect to variation in dependent parameters declared in the physics. This is done by combining the complex Taylor series expansion within the weak formulation of the governing equation in a coupled system of linear equations forming a complex valued block matrix given by the Cauchy–Riemann matrix representation. In this work, two-dimensional linear first-order elements have been implemented in ZFEM to predict the design derivatives of the mechanical displacement field and the voltage potential field for a layered piezoelectric system in a steady-state study with Dirichlet boundary condition applied at the top and bottom edges of the geometry. This approach allows the standard FEM solution to quantify the sensitivity of the mechanical displacement and voltage potential fields with respect to small variations in the material properties through the information obtained from the computation of the derivatives. The domain is formed by a layered body with PZT-4 and PZT-5 stacked together. For result verification, the numerical solution obtained with ZFEM was compared to results from a commercial FEM package and the solution from the imaginary part was compared to the exact solution of a well-known benchmark problem. In conclusion, comparison of the results showed good agreement for both the real and imaginary parts of the solution and the largest sensitivities were found in PZT-5 specifically in C 13 , C 33 , and ε 33 .

42 ENGINEERING↗
NASA NTRS1999

Flutter and Forced Response Analyses of Cascades using a Two-Dimensional Linearized Euler Solver

Flutter and forced response analyses for a cascade of blades in subsonic and transonic flow is presented. The structural model for each blade is a typical section with bending and torsion degrees of freedom. The unsteady aerodynamic forces due to bending and torsion motions. and due to a vortical gust disturbance are obtained by solving unsteady linearized Euler equations. The unsteady linearized equations are obtained by linearizing the unsteady nonlinear equations about the steady flow. The predicted unsteady aerodynamic forces include the effect of steady aerodynamic loading due to airfoil shape, thickness and angle of attack. The aeroelastic equations are solved in the frequency domain by coupling the un- steady aerodynamic forces to the aeroelastic solver MISER. The present unsteady aerodynamic solver showed good correlation with published results for both flutter and forced response predictions. Further improvements are required to use the unsteady aerodynamic solver in a design cycle.

Reddy, T. S. R.↗
NASA NTRS1994

Reionization in a cold dark matter universe: The feedback of galaxy formation on the intergalactic medium

We study the coupled evolution of the intergalactic medium (IGM) and the emerging structure in the universe in the context of the cold dark matter (CDM) model, with a special focus on the consequences of imposing reionization and the Gunn-Peterson constraint as a boundary condition on the model. We have calculated the time-varying density of the IGM by coupling our detailed, numerical calculations of the thermal and ionization balance and radiative transfer in a uniform, spatially averaged IGM of H and He, including the mean opacity of an evolving distribution of gas clumps which correspond to quasar absorption line clouds, to the linearized equations for the growth of density fluctuations in both the gaseous and dark matter components in a CDM universe. We use the linear growth equations to identify the fraction of the gas which must have collapsed out at each epoch, an approach similar in spirit to the so-called Press-Schechter formalism. We identify the IGM density with the uncollapsed baryon fraction. The collapsed fraction is postulated to be a source of energy injection into the IGM, by radiation or bulk hydrodynamical heating (e.g., via shocks) or both, at a rate which is marginally enough to satisfy the Gunn-Peterson constraint at z less than 5. Our results include the following: (1) We find that the IGM in a CDM model must have contained a substantial fraction of the total baryon density of the universe both during and after its reionization epoch. (2) As a result, our previous conclusion that the observed Quasi-Stellar Objects (QSOs) at high redshift are not sufficient to ionize the IGM enough to satisfy the Gunn-Peterson constraint is confirmed. (3) We predict a detectable He II Gunn-Peterson effect at 304(1 + z) A in the spectra of quasars at a range of redshift z greater than or approx. 3, depending on the nature of the sources of IGM reionization. (4) We find, moreover, that a CDM model with high bias parameter b (i.e., b greater than or approx. 2) cannot account for the baryon content of the universe at z approximately 3 observed in quasar absorption line gas unless Omega (sub B) significantly exceeds the maximum value allowed by big bang nucleocynthesis. (5) For a CDM model with bias parameter within the allowed range of (lower) values, the lower limit to Omega(sub B) imposed by big bang nucleosynthesis (Omega(sub B) h(sup 2) greater than or equal to 0.01) combines with our results to yield the minimum IGM density for the CDM fodel. For CDM with b = 1 (Cosmic Background Explorer (COBE) normalization), we find Omega(sub IGM)(sup min) (z approximately 4) approx. equal 0.02-0.03, and Omega(sub IGM)(sup min)(z approximately 0) approx. equal 0.005-0.03, depending upon the nature of the sources of IGM reionization. (6) In general, we find that self-consistent reionization of the IGM by the collapsed baryon fraction has a strong effect on the rate of collapse. (7) As a further example, we show that the feedback effect on the IGM of energy release by the collapsed baryon fraction may explain the slow evolution of the observed comoving QSO number density between z = 5 and z = 2, followed by the sharp decline after z = 2.

Shapiro, Paul R.↗
NASA NTRS1983

Proposed solution methodology for the dynamically coupled nonlinear geared rotor mechanics equations

The equations which describe the three-dimensional motion of an unbalanced rigid disk in a shaft system are nonlinear and contain dynamic-coupling terms. Traditionally, investigators have used an order analysis to justify ignoring the nonlinear terms in the equations of motion, producing a set of linear equations. This paper will show that, when gears are included in such a rotor system, the nonlinear dynamic-coupling terms are potentially as large as the linear terms. Because of this, one must attempt to solve the nonlinear rotor mechanics equations. A solution methodology is investigated to obtain approximate steady-state solutions to these equations. As an example of the use of the technique, a simpler set of equations is solved and the results compared to numerical simulations. These equations represent the forced, steady-state response of a spring-supported pendulum. These equations were chosen because they contain the type of nonlinear terms found in the dynamically-coupled nonlinear rotor equations. The numerical simulations indicate this method is reasonably accurate even when the nonlinearities are large.

Mitchell, L. D.↗
NASA NTRS1992

A numerical method for solving systems of linear ordinary differential equations with rapidly oscillating solutions

The present numerical method for accurate and efficient solution of systems of linear equations proceeds by numerically developing a set of basis solutions characterized by slowly varying dependent variables. The solutions thus obtained are shown to have a computational overhead largely independent of the small size of the scale length which characterizes the solutions; in many cases, the technique obviates series solutions near singular points, and its known sources of error can be easily controlled without a substantial increase in computational time.

Bernstein, Ira B.↗
NASA NTRS1996

Scalable Solutions to Integral Equation and Finite Element Simulations

The application of advanced computer architecture and software to a broad range of electromagnetic problems had allowed more accurate simulations of electrically larger and more complex components and systems...This paper overviews solutions to Maxwell's equations implicitly defined through systems of linear equations.

parallel↗
NASA NTRS1978

Aeromechanical stability of helicopters with a bearingless main rotor. Part 1: Equations of motion

Equations of motion for a coupled rotor-body system were derived for the purpose of studying air and ground resonance characteristics of helicopters that have bearingless main rotors. For the fuselage, only four rigid body degrees of freedom are considered; longitudinal and lateral translations, pitch, and roll. The rotor is assumed to consist of three or more rigid blades. Each blade is joined to the hub by means of a flexible beam segment (flexbeam or strap). Pitch change is accomplished by twisting the flexbeam with the pitch-control system, the characteristics of which are variable. Thus, the analysis is capable of implicitly treating aeroelastic couplings generated by the flexbeam elastic deflections, the pitch-control system, and the angular offsets of the blade and flexbeam. The linearized equations are written in the nonrotating system retaining only the cyclic rotor modes; thus, they comprise a system of homogeneous ordinary differential equations with constant coefficients. All contributions to the linearized perturbation equations from inertia, gravity, quasi-steady aerodynamics, and the flexbeam equilibrium deflections are retained exactly.

Hodges, D. H.↗
DOE OSTI2021

Linear solvers for power grid optimization problems: A review of GPU-accelerated linear solvers

The linear equations that arise in interior methods for constrained optimization are sparse symmetric indefinite, and they become extremely ill-conditioned as the interior method converges. These linear systems present a challenge for existing solver frameworks based on sparse LU or LDL T decompositions. Here, we benchmark five well known direct linear solver packages on CPU- and GPU-based hardware, using matrices extracted from power grid optimization problems. The achieved solution accuracy varies greatly among the packages. None of the tested packages delivers significant GPU acceleration for our test cases. For completeness of the comparison we include results for MA57, which is one of the most efficient and reliable CPU solvers for this class of problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

HHL algorithm with mapping function and enhanced sampling for model predictive control in microgrids

Here, this paper presents a refined quantum Harrow Hassidim Lloyd (HHL) algorithm for microgrid control. The first novelty of the developed method is that a mapping shift function enables the original HHL algorithm to handle general linear equations with non-singular and indefinite matrix. Second, a method of Matrix Extension for Amplifying Sampling Probabilities of Intended Solution (ME-ASPI) is proposed to design the reformulated linear algebraic equations, allowing for improved sampling efficiency of the quantum tomography in the refined HHL algorithm. Then, we applied the method to solve the model predictive control (MPC) problem in nonlinear dynamical microgrids. Specifically, with the ME-ASPI method, the refined HHL algorithm can effectively obtain the intended partial optimal control inputs for MPC. The optimization of quadratic programming problem in each time step of MPC is transformed into a linear system problem, which is addressed by the proposed quantum solver through using only partial information, with the time complexity improved from $\mathscr{O}(\mathscr{N}^{2.37286})$ classically to $\mathscr{O}(\mathscr{N}^{2} log \mathscr{N}$ x $p$ log $p)$ in quantum. Numerical examples have validated the effectiveness of the refined HHL algorithm with the proposed mapping function and the ME-ASPI method. By leveraging quantum properties, the proposed method provides a hybrid quantum–classical framework for microgrid control. This generic method can also potentially tackle many other challenges in analyzing and controlling general complex engineered systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1971

Theory of the unmagnetized plasma.

The Vlasov mathematical model of a plasma, which has come to be thought more useful than any other in describing the dynamical behavior of the majority of plasmas of interest, is first examined. Macroscopic variables and moment equations; linear electrostatics solutions; plasma oscillations, ion acoustic waves, and linear instabilities are treated, as well as external fields, 'test' charges, and nonlinear Vlasov phenomena. Plasmas are statistically described, and attention is given to the kinetic theory of the stable, uniform plasma and the Balescu-Lenard equation; two-time ensemble averages and fluctuation spectra in stable plasmas; the kinetic theory of the unstable plasma; and ensembles of Vlasov plasmas. Some illustrative experiments are described. Four appendixes deal with the electrostatic approximation and transverse waves; solution of the linearized Vlasov equation in a magnetic field; estimates of correlation functions from thermal equilibrium; and equivalence of spatially uniform BBGKY and Klimontovich correlations.

Montgomery, D. C.↗