Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal 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 1,189 records · Page 66

Online Primal-Dual Methods With Measurement Feedback for Time-Varying Convex Optimization

This paper addresses the design and analysis of feedback-based online algorithms to control systems or networked systems based on performance objectives and engineering constraints that may evolve over time. The emerging timevarying convex optimization formalism is leveraged to model optimal operational trajectories of the systems, as well as explicit local and network-level operational constraints. Departing from existing batch and feed-forward optimization approaches, the design of the algorithms capitalizes on an online implementation of primal-dual projected-gradient methods; the gradient steps are, however, suitably modified to accommodate feedback from the system in the form of measurements - hence, the term 'online optimization with feedback.' By virtue of this approach, the resultant algorithms can cope with model mismatches in the algebraic representation of the system states and outputs, they avoid pervasive measurements of exogenous inputs, and they naturally lend themselves to a distributed implementation. Under suitable assumptions, analytical convergence claims are established in terms of dynamic regret. Furthermore, when the synthesis of the feedback-based online algorithms is based on a regularized Lagrangian function, Q-linear convergence to solutions of the time-varying optimization problem is shown.

97 MATHEMATICS AND COMPUTING↗

Aerospace engineering design by systematic decomposition and multilevel optimization

A method for systematic analysis and optimization of large engineering systems, by decomposition of a large task into a set of smaller subtasks that is solved concurrently is described. The subtasks may be arranged in hierarchical levels. Analyses are carried out in each subtask using inputs received from other subtasks, and are followed by optimizations carried out from the bottom up. Each optimization at the lower levels is augmented by analysis of its sensitivity to the inputs received from other subtasks to account for the couplings among the subtasks in a formal manner. The analysis and optimization operations alternate iteratively until they converge to a system design whose performance is maximized with all constraints satisfied. The method, which is still under development, is tentatively validated by test cases in structural applications and an aircraft configuration optimization.

Sobieszczanski-Sobieski, J.↗

SU-E-T-182: Clinical Implementation of TG71-Based Electron MU Calculation and Comparison with a Commercial Secondary Calculation

Purpose: The TG-71 report was published in 2014 to present standardized methodologies for MU calculations and determination of dosimetric quantities. This work explores the clinical implementation of a TG71-based electron MU calculation algorithm and compares it with a recently released commercial secondary calculation program–Mobius3D (Mobius Medical System, LP). Methods: TG-71 electron dosimetry data were acquired, and MU calculations were performed based on the recently published TG-71 report. The formalism in the report for extended SSD using air-gap corrections was used. The dosimetric quantities, such PDD, output factor, and f-air factors were incorporated into an organized databook that facilitates data access and subsequent computation. The Mobius3D program utilizes a pencil beam redefinition algorithm. To verify the accuracy of calculations, five customized rectangular cutouts of different sizes–6×12, 4×12, 6×8, 4×8, 3×6 cm{sup 2}–were made. Calculations were compared to each other and to point dose measurements for electron beams of energy 6, 9, 12, 16, 20 MeV. Each calculation / measurement point was at the depth of maximum dose for each cutout in a 10×10 cm{sup 2} or 15×15cm{sup 2} applicator with SSDs 100cm and 110cm. Validation measurements were made with a CC04 ion chamber in a solid water phantom for electron beams of energy 9 and 16 MeV. Results: Differences between the TG-71 and the commercial system relative to measurements were within 3% for most combinations of electron energy, cutout size, and SSD. A 5.6% difference between the two calculation methods was found only for the 6MeV electron beam with 3×6 cm{sup 2}cutout in the 10×10{sup 2}cm applicator at 110cm SSD. Both the TG-71 and the commercial calculations show good consistency with chamber measurements: for 5 cutouts, <1% difference for 100cm SSD, and 0.5–2.7% for 110cm SSD. Conclusions: Based on comparisons with measurements, a TG71-based computation method and a Mobius3D program produce reasonably accurate MU calculations for electron-beam therapy.

60 APPLIED LIFE SCIENCES↗

Progress on GPS standardization

It has been clear for some time that a desirable and necessary step for improvement of the accuracy of GPS time comparisons is to establish GPS standards which may be adopted by receiver designers and users. For this reason, a formal body, the CCDS Group on GPS Time Transfer Standards (CGGTTS), was created in 1991. It operates under the auspices of the permanent CCDS Working Group on TAI, with the objective of recommending procedures and models for operational time transfer by the GPS common-view method. It works in close cooperation with the Subcommittee on Time of the Civil GPS Service Interface Committee. The members of the CGGTTS have met in December 1991 and in June 1992. Following these two formal meetings, a number of decisions were taken for unifying the treatment of GPS short-term data and for standardizing the format of GPS data files. A formal CGGTTS Recommendation is now being written concerning these points. This paper relates on the work carried out by the CGGTTS.

Thomas, C.↗

How to use the Standard Model effective field theory

Here, we present a practical three-step procedure of using the Standard Model effective field theory (SM EFT) to connect ultraviolet (UV) models of new physics with weak scale precision observables. With this procedure, one can interpret precision measurements as constraints on a given UV model. We give a detailed explanation for calculating the effective action up to one-loop order in a manifestly gauge covariant fashion. This covariant derivative expansion method dramatically simplifies the process of matching a UV model with the SM EFT, and also makes available a universal formalism that is easy to use for a variety of UV models. A few general aspects of RG running effects and choosing operator bases are discussed. Finally, we provide mapping results between the bosonic sector of the SM EFT and a complete set of precision electroweak and Higgs observables to which present and near future experiments are sensitive. Many results and tools which should prove useful to those wishing to use the SM EFT are detailed in several appendices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thermodynamically consistent algorithms for models of incompressible multiphase polymer solutions with a variable mobility

Here we present a general strategy for developing structure and property preserving numerical algorithms for thermodynamically consistent models of incompressible multiphase polymer solutions with a variable mobility. We first present a formalism to derive thermodynamically consistent, incompressible, multiphase polymer models. Then, we develop the general strategy, known as the supplementary variable method, to devise thermodynamically consistent numerical approximations to the models. We illustrate the numerical strategy using newly developed models of incompressible diblock copolymer solutions coupled with an electric and a magnetic field, respectively. Mesh refinement is conducted to verify convergence rates of the developed schemes. Some numerical examples are given to exhibit underlying dynamics absent from and driven by the external fields, respectively, highlighting differences between models with the variable and constant mobilities.

97 MATHEMATICS AND COMPUTING↗

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

74 ATOMIC AND MOLECULAR PHYSICS↗

Resurgence and Lefschetz thimble in three-dimensional $\mathcal{N}=2$ supersymmetric Chern–Simons matter theories

We study a certain class of supersymmetric (SUSY) observables in three-dimensional N=2 SUSY Chern–Simons (CS) matter theories and investigate how their exact results are related to the perturbative series with respect to coupling constants given by inverse CS levels. We show that the observables have non-trivial resurgent structures by expressing the exact results as a full transseries consisting of perturbative and non-perturbative parts. As real mass parameters are varied, we encounter Stokes phenomena at an infinite number of points, where the perturbative series becomes non-Borel-summable due to singularities on the positive real axis of the Borel plane. We also investigate the Stokes phenomena when the phase of the coupling constant is varied. For these cases, we find that the Borel ambiguities in the perturbative sector are canceled by those in non-perturbative sectors and end up with an unambiguous result which agrees with the exact result even on the Stokes lines. We also decompose the Coulomb branch localization formula, which is an integral representation for the exact results, into Lefschetz thimble contributions and study how they are related to the resurgent transseries. We interpret the non-perturbative effects appearing in the transseries as contributions of complexified SUSY solutions which formally satisfy the SUSY conditions but are not on the original path integral contour.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamics of a qubit in a high-impedance transmission line from a bath perspective

We investigate the quantum dynamics of a generic model of light-matter interaction in the context of high-impedance waveguides, focusing on the behavior of the photonic states generated in the waveguide. The model treated consists simply of a two-level system coupled to a bosonic bath (the Ohmic spin-boson model). Quantum quenches as well as scattering of an incident coherent pulse are studied using two complementary methods. First, we develop an approximate ansatz for the electromagnetic waves based on a single multimode coherent state wave function; formally, this approach combines in a single framework ideas from adiabatic renormalization, the Born-Markov approximation, and input-output theory. Second, we present numerically exact results for scattering of a weak intensity pulse by using numerical renormalization group (NRG) calculations. NRG provides a benchmark for any linear response property throughout the ultrastrong-coupling regime. We find that in a sudden quantum quench, the coherent state approach produces physical artifacts, such as improper relaxation to the steady state. These previously unnoticed problems are related to the simplified form of the ansatz that generates spurious correlations within the bath. In the scattering problem, NRG is used to find the transmission and reflection of a single photon, as well as the inelastic scattering of that single photon. Simple analytical formulas are established and tested against the NRG data that predict quantitatively the transport coefficients for up to moderate environmental impedance. These formulas resolve pending issues regarding the presence of inelastic losses in the spin-boson model near absorption resonances, and could be used for comparison to experiments in Josephson waveguide quantum electrodynamics. Finally, the scattering results using the coherent state wave-function approach are compared favorably to the NRG results for very weak incident intensity. Finally, we end our study by presenting results at higher power where the response of the system is nonlinear.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Phonon thermal conductance across GaN-AlN interfaces from first principles

The vibrational thermal conductances (G) across GaN-AlN interfaces are computed using a nonequilibrium Green's function formalism in the harmonic limit with bulk and interfacial interatomic force constants (IFCs) fully from density functional theory. Several numerical methods and supercell configurations are employed to examine the sensitivity of G to variances of IFCs. In particular, the effects of supercell size, the enforcement of symmetry constraints, and truncation of IFCs near the interface, and atomic relaxation on phonon transmission and conductance are explored. Our fully first-principles calculations are compared with common approximations and measured G values inferred from thermal conductivity measurements for GaN-AlN superlattices. Here, our calculated value, G~300MWm -2 K -1 , is nearly half that from measurements. This discrepancy is critically analyzed in terms of the physical assumptions of the calculations and the derivation of the experimental values. This work provides guidelines to determine “physically correct” sets of interfacial IFCs from first principles for thermal conductance calculations using minimal computational resources. It also contributes toward developing predictive calculations and a more complete picture of thermal conduction across interfaces, a step toward first-principles multiscale thermal transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DELTACON: A Principled Massive-Graph Similarity Function with Attribution

How much did a network change since yesterday? How different is the wiring between Bob's brain (a left-handed male) and Alice's brain (a right-handed female)? Graph similarity with known node correspondence, i.e. the detection of changes in the connectivity of graphs, arises in numerous settings. In this work, we formally state the axioms and desired properties of the graph similarity functions, and evaluate when state-of-the-art methods fail to detect crucial connectivity changes in graphs. We propose DeltaCon, a principled, intuitive, and scalable algorithm that assesses the similarity between two graphs on the same nodes (e.g. employees of a company, customers of a mobile carrier). In our experiments on various synthetic and real graphs we showcase the advantages of our method over existing similarity measures. We also employ DeltaCon to real applications: (a) we classify people to groups of high and low creativity based on their brain connectivity graphs, and (b) do temporal anomaly detection in the who-emails-whom Enron graph.

97 MATHEMATICS AND COMPUTING↗

Vortex equations: Singularities, numerical solution, and axisymmetric vortex breakdown

A method of weighted residuals for the computation of rotationally symmetric quasi-cylindrical viscous incompressible vortex flow is presented and used to compute a wide variety of vortex flows. The method approximates the axial velocity and circulation profiles by series of exponentials having (N + 1) and N free parameters, respectively. Formal integration results in a set of (2N + 1) ordinary differential equations for the free parameters. The governing equations are shown to have an infinite number of discrete singularities corresponding to critical values of the swirl parameters. The computations point to the controlling influence of the inner core flow on vortex behavior. They also confirm the existence of two particular critical swirl parameter values: one separates vortex flow which decays smoothly from vortex flow which eventually breaks down, and the second is the first singularity of the quasi-cylindrical system, at which point physical vortex breakdown is thought to occur.

Bossel, H. H.↗

Pilot interaction with automated airborne decision making systems

The role of the pilot and crew for future aircraft is discussed. Fifteen formal experimental studies and the development of a variety of models of human behavior based on queueing history, pattern recognition methods, control theory, fuzzy set theory, and artificial intelligence concepts are presented. L.F.M.

Rouse, W. B.↗

Solution influence on biomolecular equilibria - Nucleic acid base associations

Various attempts to construct an understanding of the influence of solution environment on biomolecular equilibria at the molecular level using computer simulation are discussed. First, the application of the formal statistical thermodynamic program for investigating biomolecular equilibria in solution is presented, addressing modeling and conceptual simplications such as perturbative methods, long-range interaction approximations, surface thermodynamics, and hydration shell. Then, Monte Carlo calculations on the associations of nucleic acid bases in both polar and nonpolar solvents such as water and carbon tetrachloride are carried out. The solvent contribution to the enthalpy of base association is positive (destabilizing) in both polar and nonpolar solvents while negative enthalpies for stacked complexes are obtained only when the solute-solute in vacuo energy is added to the total energy. The release upon association of solvent molecules from the first hydration layer around a solute to the bulk is accompanied by an increase in solute-solvent energy and decrease in solvent-solvent energy. The techniques presented are expectd to displace less molecular and more heuristic modeling of biomolecular equilibria in solution.

Pohorille, A.↗

Orbital stability of the unseen solar companion linked to periodic extinction events

Evidence from three-dimensional numerical modelling is presented that only cometary orbits with a limited range in inclination with respect to the galactic plane are formally stable for the length of time required to cause periodic extinction events. The calculations were done using Cowell's method employing a fourth-order Runge-Kutta integration scheme in an inertial reference frame in orbit about the Galaxy. Tidal perturbations in the radial direction due to the Galaxy and the Coriolis forces are included. The vertical component of the gravitational field of the galactic disk is superimposed on these forces. The results indicate that orbits for Nemesis that are inclined at more than 30 deg to the galactic plane are not allowed and suggests that the search for Nemesis should be concentrated toward the plane of the Galaxy. Perturbations by passing stars or molecular clouds may make even the low-inclination orbits unstable.

Torbett, M. V.↗

A mathematical theory of learning control for linear discrete multivariable systems

When tracking control systems are used in repetitive operations such as robots in various manufacturing processes, the controller will make the same errors repeatedly. Here consideration is given to learning controllers that look at the tracking errors in each repetition of the process and adjust the control to decrease these errors in the next repetition. A general formalism is developed for learning control of discrete-time (time-varying or time-invariant) linear multivariable systems. Methods of specifying a desired trajectory (such that the trajectory can actually be performed by the discrete system) are discussed, and learning controllers are developed. Stability criteria are obtained which are relatively easy to use to insure convergence of the learning process, and proper gain settings are discussed in light of measurement noise and system uncertainties.

Phan, Minh↗

Periodic Comet Showers, Mass Extinctions, and the Galaxy

Geologic data on mass extinctions of life and evidence of large impacts on the Earth are thus far consistent with a quasi-periodic modulation of the flux of Oort cloud comets. Impacts of large comets and asteroids are capable of causing mass extinction of species, and the records of large impact craters and mass show a correlation. Impacts and extinctions display periods in the range of approximately 31 +/- 5 m.y., depending on dating methods, published time scales, length of record, and number of events analyzed. Statistical studies show that observed differences in the formal periodicity of extinctions and craters are to be expected, taking into consideration problems in dating and the likelihood that both records would be mixtures of periodic and random events. These results could be explained by quasi-periodic showers of Oort Cloud comets with a similar cycle. The best candidate for a pacemaker for comet showers is the Sun's vertical oscillation through the plane of the Galaxy, with a half-period over the last 250 million years in the same range. We originally suggested that the probability of encounters with molecular clouds that could perturb the Oort comet cloud and cause comet showers is modulated by the Sun's vertical motion through the galactic disk. Tidal forces produced by the overall gravitational field of the Galaxy can also cause perturbations of cometary orbits. Since these forces vary with the changing position of the solar system in the Galaxy, they provide a mechanism for the periodic variation in the flux of Oort cloud comets into the inner solar system. The cycle time and degree of modulation depend critically on the mass distribution in the galactic disk. Additional information is contained in the original extended abstract.

Rampino, M. R.↗

NASA Astronauts on Soyuz: Experience and Lessons for the Future

The U. S., Russia, and, China have each addressed the question of human-rating spacecraft. NASA's operational experience with human-rating primarily resides with Mercury, Gemini, Apollo, Space Shuttle, and International Space Station. NASA s latest developmental experience includes Constellation, X38, X33, and the Orbital Space Plane. If domestic commercial crew vehicles are used to transport astronauts to and from space, Soyuz is another example of methods that could be used to human-rate a spacecraft and to work with commercial spacecraft providers. For Soyuz, NASA's normal assurance practices were adapted. Building on NASA's Soyuz experience, this report contends all past, present, and future vehicles rely on a range of methods and techniques for human-rating assurance, the components of which include: requirements, conceptual development, prototype evaluations, configuration management, formal development reviews (safety, design, operations), component/system ground-testing, integrated flight tests, independent assessments, and launch readiness reviews. When constraints (cost, schedule, international) limit the depth/breadth of one or more preferred assurance means, ways are found to bolster the remaining areas. This report provides information exemplifying the above safety assurance model for consideration with commercial or foreign-government-designed spacecraft. Topics addressed include: U.S./Soviet-Russian government/agency agreements and engineering/safety assessments performed with lessons learned in historic U.S./Russian joint space ventures

Source record↗