Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “OPERATIONAL CALCULUS”

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.

65 records · Page 4

Mathematical Modeling of Diverse Phenomena

Tensor calculus is applied to the formulation of mathematical models of diverse phenomena. Aeronautics, fluid dynamics, and cosmology are among the areas of application. The feasibility of combining tensor methods and computer capability to formulate problems is demonstrated. The techniques described are an attempt to simplify the formulation of mathematical models by reducing the modeling process to a series of routine operations, which can be performed either manually or by computer.

Howard, J. C.↗

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation and Exploitation of Aggregation Abstractions for Scheduling and Resource Allocation

Our research is investigating abstraction of computational theories for scheduling and resource allocation. These theories are represented in a variant of first order predicate calculus, parameterized multisorted logic, that facilitates specification of large problems. A particular problem is conceptually stated as a set of ground sentences that are consistent with a quantified theory. We are mainly investigating the automated generation of aggregation abstractions and approximations in which detailed resource allocation constraints are replaced by constraints between aggregate demand and capacity. We are also investigating the interaction of aggregation abstractions with the more thoroughly investigated abstractions of weakening operator preconditions. The purpose of the theories for aggregated demand/capacity is threefold: first, to answer queries about aggregate properties, such as gross feasibility; second, to reduce computational costs by using the solution of aggregate problems to guide the solution of detailed problems; and third, to facilitate reformulating theories to approximate problems for which there are efficient problem solving methods. We also describe novel methods for exploiting aggregation abstractions.

Linden, Theodore A.↗

On Dirac equations for linear magnetoacoustic waves propagating in an isothermal atmosphere

A new analytical approach to study linear magnetoacoustic waves propagating in an isothermal, stratified, and uniformly magnetized atmosphere is presented. The approach is based on Dirac equations, and the theory of Sturm-Liouville operators is used to investigate spectral properties of the obtained Dirac Hamiltonians. Two cases are considered: (1) the background magnetic field is vertical, and the waves are separated into purely magnetic (transverse) and purely acoustic (longitudinal) modes; and (2) the field is tilted with respect to the vertical direction and the magnetic and acoustic modes become coupled giving magnetoacoustic waves. For the first case, the Dirac Hamiltonian possesses either a discrete spectrum, which corresponds to standing magnetic waves, or a continuous spectrum, which can be clearly identified with freely propagating acoustic waves. For the second case, the quantum mechanical perturbation calculus is used to study coupling and energy exchange between the magnetic and acoustic components of magnetoacoustic waves. It is shown that this coupling may efficiently prevent trapping of magnetoacoustic waves instellar atmospheres.

Alicki, R.↗

Generation and exploration of aggregation abstractions for scheduling and resource allocation

This paper presents research on the abstraction of computational theories for scheduling and resource allocation. The paper describes both theory and methods for the automated generation of aggregation abstractions and approximations in which detailed resource allocation constraints are replaced by constraints between aggregate demand and capacity. The interaction of aggregation abstraction generation with the more thoroughly investigated abstractions of weakening operator preconditions is briefly discussed. The purpose of generating abstract theories for aggregated demand and resources includes: answering queries about aggregate properties, such as gross feasibility; reducing computational costs by using the solution of aggregate problems to guide the solution of detailed problems; facilitating reformulating theories to approximate problems for which there are efficient problem-solving methods; and reducing computational costs of scheduling by providing more opportunities for variable and value-ordering heuristics to be effective. Experiments are being developed to characterize the properties of aggregations that make them cost effective. Both abstract and concrete theories are represented in a variant of first-order predicate calculus, which is a parameterized multi-sorted logic that facilitates specification of large problems. A particular problem is conceptually represented as a set of ground sentences that is consistent with a quantified theory.

Lowry, Michael R.↗

Numerical procedures for the calculation of the stresses in monocoques III : calculation of the bending moments in fuselage frames

This report deals with the calculation of the bending moments in and the distortions of fuselage rings upon which known concentrated and distributed loads are acting. In the procedure suggested, the ring is divided into a number of beams each having a constant radius of curvature. The forces and moments caused in the end sections of the beams by individual unit displacements of the end sections are listed in a table designated as the operations table in conformity with Southwell's nomenclature. The operations table and the external loads are equivalent to a set of linear equations. For their solution the following three procedures are presented: 1) Southwell's method of systematic relaxations. This is a step-by-step approximation procedure guided by the physical interpretation of the changes in the values of the unknown. 2) The growing unit procedure in which the individual beams are combined successively into beams of increasing length until finally the entire ring becomes a single beam. In each step of the procedure a set of not more than three simultaneous linear equations is solved. 3) Solution of the entire set of simultaneous equations by the methods of the matrix calculus. In order to demonstrate the manner in which the calculations may be carried out, the following numerical examples are worked out: 1) Curved beam with both its end sections rigidly fixed. The load is a concentrated force. 2) Egg-shape ring with symmetric concentrated loads. 3) Circular ring with antisymmetric concentrated loads and shear flow (torsion of the fuselage). 4) Same with V-braces incorporated in the ring. 5) Egg-shape ring with antisymmetric concentrated loads and shear flow (torsion of the fuselage). 6) Same with V-braces incorporated in the ring. The results of these calculations are checked, whenever possible, by calculations carried out according to known methods of analysis. The agreement is found to be good. The amount of work necessary for the solution of ring problems by the methods described in the present report is practically independent of the degree of redundancy of the structure. For this reason the methods are recommended for use particularly in problems of rings having one or more internal bracing elements.

LOADS - FUSELAGES↗

Quantum Algorithm for Linear Non-unitary Dynamics with Near-Optimal Dependence on All Parameters

We introduce a family of identities that express general linear non-unitary evolution operators as a linear combination of unitary evolution operators, each solving a Hamiltonian simulation problem. This formulation can exponentially enhance the accuracy of the recently introduced linear combination of Hamiltonian simulation (LCHS) method [An, Liu, and Lin, Physical Review Letters, 2023]. For the first time, this approach enables quantum algorithms to solve linear differential equations with both optimal state preparation cost and near-optimal scaling in matrix queries on all parameters.

Applied Dynamical Systems↗

An Influenced Future

Phobos is a vital precursor and catalyst before our next giant leap to Mars. The principle period of a Phobos' mission could be a series of robotic precursor missions for experimental perception, soil examination, ecological approval and landing site distinguishing proof. For my summer intern position at Johnson Space Center I chipped away at creating a GUNNS (General Use Nodal Network Solver) based power subsystem model for the miniATHLETE hopper, which is a conceptual (idea-based) robotic lander that will operate on Phobos. Keeping in mind the end goal to begin on my venture, I needed to comprehend my undertaking before whatever else, in which I concentrated on C++ to see how to implement the code that GUNNS generates to a Trick S_define file. Prior to coming to this internship at Johnson Space Center Dr. Edwin Zack Crues provided a class on modeling and simulation, which introduced me to the Trick simulation environment. The goal of my project was to develop a GUNNS based power subsystem model for the miniATHLETE hopper. The model needed to incorporate a solar array, battery, hopping legs, and onboard scientific instruments (sensitive measuring/recording devices). The secondary bonus goal after I completed the electrical aspect of my model was to develop a GUNNS based thermal subsystem model for the miniATHLETE hopper. Stringing the two aspects together, I would need to code up a signal aspect to make the system work as one. Accomplishing my goals would not be an easy thing, however I had successfully completed the electrical aspect model with twenty-four servos, six cameras, and multiple sensors. Venturing to complete my project has eluded me to many failures in my design to tune many things like the battery to the proper voltage and the load to the proper wattage. During this time I had touched up on advanced topics in calculus in which I implemented in the converter in my electrical model. I am currently working with my mentor Zu Qun Li to create a signal aspect to control the temperatures inside my electrical aspect model. During my time at JSC I had effectively figured out how to create subsystem models utilizing Trick and GUNNS. I obtain essential knowledge of power and thermal subsystem design for a robotic vehicle. I also learned how to work and communicate in a team effectively to accomplish a goal. Before coming to Johnson Space Center my future career and educational goals included uncertainty, however now I have a completely new look on my path to a prosperous future. My NASA experience has unquestionably impacted me to accomplish and surpass my own particular desires. After my time at Johnson Space Center I plan to apply for a coop position for NASA. This has been a dream come true that I adored each moment being at JSC realizing that I am far fit for doing things most individuals can just long for.

Owens, Christopher↗

Quasilinear theory: the lost ponderomotive effects and why they matter

Quasilinear theory (QLT) has been used for modeling wave–plasma interactions for decades but remains largely heuristic. Plasma inhomogeneity, ponderomotive effects, microscopic fluctuations, and collisions are not easily accommodated from first principles in QLT, and typically are ignored entirely, due to the limitations of the standard Fourier–Laplace global-mode approach. This results in inconsistencies, for example, violation of the action conservation for nonresonant waves. However, these issues can be avoided, and the theory can be substantially generalized and corrected, if QLT is formulated using more suitable analytical tools, particularly, the Weyl symbol calculus. Here, an attempt is made to deliver an accessible review of this modern formulation, provide intuitive calculations for special cases, and elaborate on the connection with the ‘oscillation-center QLT’ originally proposed by Dewar (Phys Fluids 16:1102, 1973). A Fokker–Planck equation for a ‘dressed’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, ponderomotive forces, and interactions with background fields for a generic Hamiltonian, so many known formulations of QLT for specific plasma models become corollaries of a single unifying theory. Also, waves are allowed to be off-shell (not constrained by a dispersion relation), which allows them to accommodate microscopic fluctuations. This leads to a collision integral of the Balescu–Lenard type that has all the usual properties but is not restricted to any specific plasma model. For on-shell waves, a generalized version of the classic oscillation-center QLT is obtained. Finally, combined with the wave-kinetic equation, this formulation not only conserves particles, momentum, and energy, like the classic QLT but also reinstates the action conservation for nonresonant waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Human Exploration Missions Study Launch Window from Earth Orbit

The determination of orbital launch window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a Delta(V) penalty. Usually, because of the Delta(V) penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: (1) One impulsive maneuver from a Low Earth Orbit (LEO), (2) Two impulsive maneuvers from LEO, (3) Three impulsive maneuvers from LEO, (4) One impulsive maneuvers from a Highly Elliptical Orbit (HEO), (5) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO) The formulation of these five different launch window modes provides a rapid means of generating realistic parametric data for space exploration studies. Also the formulation provides vector and geometrical data sufficient for use as a good starting point in detail trajectory analysis based on calculus of variations, steepest descent, or parameter optimization program techniques.

Young, Archie↗

Human Mars Mission: Launch Window from Earth Orbit

The determination of orbital window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to the earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a delta V penalty. Usually, because of the delta V penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: 1) One impulsive maneuver from a Highly Elliptical Orbit (HEO); 2) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO); 3) One impulsive maneuver from a Low Earth Orbit (LEO); 4) Two impulsive maneuvers form LEO; and 5) Three impulsive maneuvers form LEO. The formulation of these five different launch window modes provides a rapid means of generating realistic parametric data for space exploration studies. Also the formulation provides vector and geometrical data sufficient for use as a good starting point in detail trajectory analysis based on calculus of variations, steepest descent, or parameter optimization program techniques.

Young, Archie↗