Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear operator theory”

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 253 records · Page 14

A new approach to the oxidational theory of mild wear

The wear and friction of 316 stainless steel was investigated in relation to oxidation at the real area of contact. An approach describing the oxidation characteristics of mild wear in terms of linear oxidation kinetics is presented. Linear oxidation kinetics more closely approaches the physical conditions involved where oxide growth on asperities is continuously interrupted by spalling and removal. The mild wear rate may now be described mathematically by two empirical equations in terms of the operational variables and the oxidational constants without using the critical oxide thickness, which is dependent on asperity size, shape and temperature, and is difficult to measure. The resulting activation energy for oxidational wear was significantly lower than that for static oxidation governed by diffusion controlled kinetics.

Hong, H.↗

Predictions of core plasma performance for the SPARC tokamak

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with Ion Cyclotron Range-of-Frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach Q~11, well above its mission objective of Q>2. To build confidence that SPARC will be successful, physics-based integrated modeling has also been performed. The TRANSP code coupled with the theory-based TGLF turbulence model and EPED predictions for pedestal stability find that Q~9 is attainable in standard H-mode operation and confirms Q>2 operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modeled with the Monte-Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

e + e – → ZH process in the SMEFT beyond leading order

We systematically study potential effects of Beyond the Standard Model physics in the e + e – → ZH process. To this end, we include all relevant dimension-6 Standard Model Effective Field Theory operators and work to next-to-leading order (NLO) accuracy in the electroweak coupling. We consider both polarized and unpolarized electron and positron beams and present results for $\sqrt {s}$ = 240, 365 and 500 GeV, emphasizing contributions where the NLO predictions differ significantly from the leading order results. At NLO, a sensitivity arises to operators that do not contribute at tree level, such as the Higgs tri-linear coupling, CP-violating operators, and dimension-6 operators involving the top quark, among many others. We compare the prospects of future e + e – colliders to explore these new physics effects with existing measurements from the LHC, electron EDMs (for CP violating operators), and Z pole measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Problem of Obtaining High Lift-Drag Ratios at Supersonic Speeds

The importance of the lift to drag ratio is well known to all aircraft designers since it gives, to a great extent, the aerodynamic efficiency of the airplane. Aerodynamic efficiency, however, is only one component of the grand compromise that a completed airplane represents. At subsonic speeds, lift-drag ratios of well over 200 have been measured in wind tunnels on airfoil sections; but few powered aircraft have attained (L/D) value of 20. It is invariably true that the requirements of stability and control, structure, and flight operation all contribute to reducing the design (L/D) (sub max)) considerably below those exotic values which can be predicted from unrestricted aerodynamic theory. If, however, a certain range or operating efficiency is required, there is most certainly a minimum (L/D) (sub max)) value for which the goals are just attainable. If we examine the range equation we see that range is proportional to the lift-drag ratio, the thermopropulsive efficiency, and the logarithm of the initial to final weight ratio. The appearance of the lift-drag ratio as a linear factor in the range equation indicates that every attempt should be made to increase (L/D) (sub max)); however, the search for higher (L/D)max may lead to strange and unorthodox configurations. Most frequently, such configurations are ruled out by the adverse effects of their geometry on the weight ratios. In the present paper, we will deal with the maximum lift-drag ratio problem for conventional configurations having a wing and a body in close proximity to each other. No attempt will be made to select a particular configuration as being the best. However, the promising direction to go from the aerodynamic view will be stressed with the understanding that the other factors may outweight the aerodynamics

Brown, Clinton E.↗

Goldstone bosons on celestial sphere and conformal soft theorems

In this paper, we study celestial amplitudes of Goldstone bosons and conformal soft theorems. Motivated by the success of soft bootstrap in momentum space and the important role of the soft limit behavior of tree-level amplitudes, our goal is to extend some of the methods to the celestial sphere. The crucial ingredient of the calculation is the Mellin transformation, which transforms four-dimensional scattering amplitudes to correlation functions of primary operators in the celestial CFT. The soft behavior of the amplitude is then translated to the singularities of the correlator. Only for amplitudes in “UV completed theories” (with sufficiently good high energy behavior) the Mellin integration can be properly performed. In all other cases, the celestial amplitude is only defined in a distributional sense with delta functions. We provide many examples of celestial amplitudes in UV-completed models, including linear sigma models and Z-theory, which is a certain completion of the SU(N) non-linear sigma model. We also comment on the BCFW-like and soft recursion relations for celestial amplitudes and the extension of soft bootstrap ideas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fault-tolerant resource comparison of qudit and qubit encodings for diagonal quadratic operators

Finite local Hilbert-space truncations arise naturally in quantum simulations of lattice field theories and motivate qudit encodings, but their fault-tolerant advantage over qubit encodings remains unclear. We compare the non-Clifford cost of implementing quadratic diagonal evolutions, exemplified by 𝑈 = 𝑒$^{−𝑖⁢𝑡⁢𝜙^2_𝑥}$ in a uniform field-amplitude discretization of a real scalar field, using either one logical 𝑑-level qudit or 𝑛 𝑏 = ⌈log 2⁡ 𝑑⌉ logical qubits. We analyze two standard settings: product-formula simulation and linear combination of unitaries (LCU) per block encoding, taking the resource metric to be the number of non-Clifford gates after synthesis into a discrete logical gate set. Because tight synthesis bounds for general single-qudit rotations are not known, we express the qudit constructions in terms of embedded two-level SU⁡(2) rotations and derive explicit finite-𝑑 break-even conditions for their synthesis cost; these serve as compiler targets for when qudit encodings can outperform the qubit baseline. Within the constructive models studied here, product-formula implementations would require an exponentially stronger per-primitive synthesis advantage for qudits to win asymptotically, while in the LCU setting the qubit encoding is asymptotically cheaper in 𝑑. Nevertheless, the finite-𝑑 threshold analysis identifies low-dimensional regions in which qudits can yield meaningful constant-factor savings, particularly for LCU-based implementations. As a secondary analysis of the LCU construction, we use an idealized negligible-overhead qubit-qudit code-switching model to give an absolute 𝑇-count comparison and reinterpret the savings as an allowable per-switch overhead budget.

Godwood, Samuel [Univ. of Liverpool (United Kingdo↗

Finding predictive models for singlet fission by machine learning

Singlet fission (SF), the conversion of one singlet exciton into two triplet excitons, could significantly enhance solar cell efficiency. Molecular crystals that undergo SF are scarce. Computational exploration may accelerate the discovery of SF materials. However, many-body perturbation theory (MBPT) calculations of the excitonic properties of molecular crystals are impractical for large-scale materials screening. We use the sure-independence-screening-and-sparsifying-operator (SISSO) machine-learning algorithm to generate computationally efficient models that can predict the MBPT thermodynamic driving force for SF for a dataset of 101 polycyclic aromatic hydrocarbons (PAH101). SISSO generates models by iteratively combining physical primary features. The best models are selected by linear regression with cross-validation. The SISSO models successfully predict the SF driving force with errors below 0.2 eV. Based on the cost, accuracy, and classification performance of SISSO models, we propose a hierarchical materials screening workflow. Three potential SF candidates are found in the PAH101 set.

36 MATERIALS SCIENCE↗

Statistical analysis of the modal properties of large structural systems.

A theory is developed to predict eigenvalue and eigenvector uncertainty in large dynamic models. The uncertainty is based on physical property uncertainty and should not be confused with numerical roundoff, although the method can be extended to include the latter. The theory, when implemented on a computer, is used to analyze the uncertainties in frequencies and mode shapes based on uncertainties in mass, stiffness, modulus of elasticity, etc. The method incorporates a linear statistical model which is quite adequate for handling property uncertainties of 10% or more. The model is not limited to small systems but uses certain statistical assumptions as well as selective matrix manipulations to keep the size of all matrix operations to within the number of degrees of freedom of the system. Examples are given for two longitudinal vibration problems, and the results are supported by a Monte Carlo analysis.

Collins, J. D.↗

Research in plasma physics

Three aspects of barium ion cloud dynamics are discussed. First, the effect of the ratio of ion cloud conductivity to background ionospheric conductivity on the motion of barium ion clouds is investigated and compared with observations of barium ion clouds. This study led to the suggestion that the conjugate ionosphere participates in the dynamics of barium ion clouds. Second, analytic work on the deformation of ion clouds is presented. Third, a linearized stability theory was extended to include the effect of the finite extent of an ion cloud, as well as the effect of the ratio of ion cloud to ionospheric conductivities. The stability properties of a plasma with contra-streaming ion beams parallel to a magnetic field are investigated. The results are interpreted in terms of parameters appropriate for collisionless shock waves. It is found that this particular instability can be operative only if the up-stream Alfven Mach number exceeds 5.5.

Source record↗

Sliding Mode Thermal Control System for Space Station Furnace Facility

The decoupled control of the nonlinear, multiinput-multioutput, and highly coupled space station furnace facility (SSFF) thermal control system is addressed. Sliding mode control theory, a subset of variable-structure control theory, is employed to increase the performance, robustness, and reliability of the SSFF's currently designed control system. This paper presents the nonlinear thermal control system description and develops the sliding mode controllers that cause the interconnected subsystems to operate in their local sliding modes, resulting in control system invariance to plant uncertainties and external and interaction disturbances. The desired decoupled flow-rate tracking is achieved by optimization of the local linear sliding mode equations. The controllers are implemented digitally and extensive simulation results are presented to show the flow-rate tracking robustness and invariance to plant uncertainties, nonlinearities, external disturbances, and variations of the system pressure supplied to the controlled subsystems.

Jackson Mark E.↗

Potential for Sonic Boom Reduction of the Boeing HSCT

The HSR sonic boom technology program includes a goal of reducing the objectionable aspects of sonic boom. Earlier HSCT sonic boom studies considered achieving significant sonic boom reduction by the use of arrow-wing planforms and detailed shaping of the airplane to produce shaped waveforms (non N-waves) at the ground. While these design efforts were largely successful, the added risk and cost of the airplanes were judged to be unacceptable. The objective of the current work is to explore smaller configuration refinements that could lead to reduced sonic boom impact, within design and operational constraints. A somewhat modest target of 10% reduction in sonic boom maximum overpressure was selected to minimize the effect on the configuration performance. This work was a joint NASA/Industry effort, utilizing the respective strengths of team members at Boeing, NASA Langley, and NASA Ames. The approach used was to first explore a wide range of modifications and airplane characteristics for their effects on sonic boom and drag, using classical Modified Linear Theory (MLT) methods. CFD methods were then used to verify promising, modifications and to analyze modifications for which the MLT methods were not appropriate. The tea m produced a list of configuration changes with their effects on sonic boom and, in some cases, an estimate of the drag penalty. The most promising modifications were applied to produce a boom-softened derivative of the baseline Boeing High Speed Civil Transport (HSCT) configuration. This boom-softened configuration was analyzed in detail for the reduce sonic boom impact and also for the effect of the configuration modifications on drag, weight, and overall performance relative to the baseline.

Haglund, George T.↗

A manual control theory analysis of vertical situation displays for STOL aircraft

Pilot-vehicle-display systems theory is applied to the analysis of proposed vertical situation displays for manual control in approach-to-landing of a STOL aircraft. The effects of display variables on pilot workload and on total closed-loop system performance was calculated using an optimal-control model for the human operator. The steep approach of an augmentor wing jet STOL aircraft was analyzed. Both random turbulence and mean-wind shears were considered. Linearized perturbation equations were used to describe longitudinal and lateral dynamics of the aircraft. The basic display configuration was one that abstracted the essential status information (including glide-slope and localizer errors) of an EADI display. Proposed flight director displays for both longitudinal and lateral control were also investigated.

Baron, S.↗

Protein crystal growth in low gravity

The solubility and growth of the protein canavalin, and the application of the schlieren technique to study fluid flow in protein crystal growth systems were investigated. These studies have resulted in the proposal of a model to describe protein crystal growth and the preliminary plans for a long-term space flight experiment. Canavalin, which may be crystallized from a basic solution by the addition of hydrogen (H+) ions, was shown to have normal solubility characteristics over the range of temperatures (5 to 25 C) and pH (5 to 7.5) studies. The solubility data combined with growth rate data gathered from the seeded growth of canavalin crystals indicated that the growth rate limiting step is a screw dislocation mechanism. A schlieren apparatus was constructed and flow patterns were observed in Rochelle salt (sodium potassium tartrate), lysozyme, and canavalin. The critical parameters were identified as the change in density with concentration (dp/dc) and the change in index of refraction with concentration (dn/dc). Some of these values were measured for the materials listed. The data for lyrozyme showed non-linearities in plots of optical properties and density vs. concentration. In conjunction with with W. A. Tiller, a model based on colloid stability theory was proposed to describe protein crystallization. The model was used to explain observations made by ourselves and others. The results of this research has lead to the development for a preliminary design for a long-term, low-g experiment. The proposed apparatus is univeral and capable of operation under microprocessor control.

Feigelson, Robert S.↗

Uncertainty reasoning in expert systems

Intelligent control is a very successful way to transform the expert's knowledge of the type 'if the velocity is big and the distance from the object is small, hit the brakes and decelerate as fast as possible' into an actual control. To apply this transformation, one must choose appropriate methods for reasoning with uncertainty, i.e., one must: (1) choose the representation for words like 'small', 'big'; (2) choose operations corresponding to 'and' and 'or'; (3) choose a method that transforms the resulting uncertain control recommendations into a precise control strategy. The wrong choice can drastically affect the quality of the resulting control, so the problem of choosing the right procedure is very important. From a mathematical viewpoint these choice problems correspond to non-linear optimization and are therefore extremely difficult. In this project, a new mathematical formalism (based on group theory) is developed that allows us to solve the problem of optimal choice and thus: (1) explain why the existing choices are really the best (in some situations); (2) explain a rather mysterious fact that fuzzy control (i.e., control based on the experts' knowledge) is often better than the control by these same experts; and (3) give choice recommendations for the cases when traditional choices do not work.

Kreinovich, Vladik↗

Solution of the three-dimensional Helmholtz equation with nonlocal boundary conditions

The Helmholtz equation is solved within a three-dimensional rectangular duct with a nonlocal radiation boundary condition at the duct exit plane. This condition accurately models the acoustic admittance at an arbitrarily-located computational boundary plane. A linear system of equations is constructed with second-order central differences for the Helmholtz operator and second-order backward differences for both local admittance conditions and the gradient term in the nonlocal radiation boundary condition. The resulting matrix equation is large, sparse, and non-Hermitian. The size and structure of the matrix makes direct solution techniques impractical; as a result, a nonstationary iterative technique is used for its solution. The theory behind the nonstationary technique is reviewed, and numerical results are presented for radiation from both a point source and a planar acoustic source. The solutions with the nonlocal boundary conditions are invariant to the location of the computational boundary, and the same nonlocal conditions are valid for all solutions. The nonlocal conditions thus provide a means of minimizing the size of three-dimensional computational domains.

Hodge, Steve L.↗

Model-Based Engine Control Architecture with an Extended Kalman Filter

This paper discusses the design and implementation of an extended Kalman filter (EKF) for model-based engine control (MBEC). Previously proposed MBEC architectures feature an optimal tuner Kalman Filter (OTKF) to produce estimates of both unmeasured engine parameters and estimates for the health of the engine. The success of this approach relies on the accuracy of the linear model and the ability of the optimal tuner to update its tuner estimates based on only a few sensors. Advances in computer processing are making it possible to replace the piece-wise linear model, developed off-line, with an on-board nonlinear model running in real-time. This will reduce the estimation errors associated with the linearization process, and is typically referred to as an extended Kalman filter. The non-linear extended Kalman filter approach is applied to the Commercial Modular Aero-Propulsion System Simulation 40,000 (C-MAPSS40k) and compared to the previously proposed MBEC architecture. The results show that the EKF reduces the estimation error, especially during transient operation.

engine control↗

Massively scalable workflows for quantum chemistry: BigChem and ChemCloud

Electronic structure theory, i.e., quantum chemistry, is the fundamental building block for many problems in computational chemistry. Here we present a new distributed computing framework (BigChem), which allows for an efficient solution of many quantum chemistry problems in parallel. BigChem is designed to be easily composable and leverages industry-standard middleware (e.g., Celery, RabbitMQ, and Redis) for distributed approaches to large scale problems. BigChem can harness any collection of worker nodes, including ones on cloud providers (such as AWS or Azure), local clusters, or supercomputer centers (and any mixture of these). BigChem builds upon MolSSI packages, such as QCEngine to standardize the operation of numerous computational chemistry programs, demonstrated here with Psi4, xtb, geomeTRIC, and TeraChem. BigChem delivers full utilization of compute resources at scale, offers a programable canvas for designing sophisticated quantum chemistry workflows, and is fault tolerant to node failures and network disruptions. We demonstrate linear scalability of BigChem running computational chemistry workloads on up to 125 GPUs. Finally, we present ChemCloud, a web API to BigChem and successor to TeraChem Cloud. ChemCloud delivers scalable and secure access to BigChem over the Internet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computation of Steady and Unsteady Laminar Flames: Theory

In this paper we describe the numerical analysis underlying our efforts to develop an accurate and reliable code for simulating flame propagation using complex physical and chemical models. We discuss our spatial and temporal discretization schemes, which in our current implementations range in order from two to six. In space we use staggered meshes to define discrete divergence and gradient operators, allowing us to approximate complex diffusion operators while maintaining ellipticity. Our temporal discretization is based on the use of preconditioning to produce a highly efficient linearly implicit method with good stability properties. High order for time accurate simulations is obtained through the use of extrapolation or deferred correction procedures. We also discuss our techniques for computing stationary flames. The primary issue here is the automatic generation of initial approximations for the application of Newton's method. We use a novel time-stepping procedure, which allows the dynamic updating of the flame speed and forces the flame front towards a specified location. Numerical experiments are presented, primarily for the stationary flame problem. These illustrate the reliability of our techniques, and the dependence of the results on various code parameters.

Hagstrom, Thomas↗