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,585 records · Page 88

Magnetocentrifugally driven flows from young stars and disks. 2: Formulation of the dynamical problem

We formulate the dynamical problem of a cool wind centrifugally driven from the magnetic interface of a young star and an adjoining Keplerian disk. We examine the situation for mildly accreting T Tauri stars that rotate slowly as well as rapidly accreting protostars that rotate near break-up. In both cases a wind can be driven from a small X-region just outside the stellar magnetopause, where the field lines assume an open geometry and are rooted to material that rotates at an angular speed equal both to the local Keplerian value and to the stellar angular speed. Assuming axial symmetry for the ideal magnetohydrodynamic flow, which requires us to postpone asking how the (lightly ionized) gas is loaded onto field lines, we can formally integrate all the governing equations analytically except for a partial equation that describes how streamlines spread in the meridional plane. Apart from the difficulty of dealing with PDEs of mixed type, finding the functional forms of the conserved quantities along streamlines - the ratio beta of magnetic field to mass flux, the specific energy H of the fluid in the rotating frame, and the total specific angular momentum J carried in the matter and the field - constitutes a standard difficulty in this kind of (Grad-Shafranov) formalism. Fortunately, because the ratio of the thermal speed of the mass-loss regions to the Keplerian speed of rotation of the interface constitutes a small parameter epsilon, we can attack the overall problem by the method of matched asymptotic expansions. This procedure leads to a natural and systematic technique for obtaining the relevant functional dependences of beta, H, and J. Moreover, we are able to solve analytically for the properties of the flow emergent from the small transsonic region driven by gas pressure without having to specify the detailed form of any of the conserved functions, beta, H, and J. This analytical solution provides inner boundary conditions for the numerical computation in a companion paper by Najita & Shu of the larger region where the main acceleration to terminal speeds occurs.

Shu, Frank H.↗

Precision Measurement of Black Hole Binary Dynamics: Analyzing the LISA Data Stream

One of the richest potential sources of insight into fundamental physics that LISA will be capable of observing is the inspiral of supermassive black hole binaries (BHBs). However, the data analysis challenge presented by the LISA data stream is quite unlike the situation for present day gravitational wave detectors. In order to make the precision measurements necessary to achieve LISA's science goals, the BHB signal must be distinguished from a data stream that not only contains instrumental noise, but potentially thousands of other signals as well, so that the "background" we wish to separate out to focus on the BHB signal is likely to be highly nonstationary and nongaussian, as well as being of scientific interest in its own right. In addition, whereas the theoretical templates that we calculate in order to ultimately estimate the parameters can afford to be somewhat inaccurate and still be effective for present day and near future detectors, this is not the case for LISA, and extremely high fidelity of the theoretical templates for high signal-to-noise signals will be required to prevent theoretical errors from dominating the parameter estimates. NVe, will describe efforts in the community of LISA data analysts to address the challenges regarding the specific issue of BHB signals. These efforts include using a Markov Chain Monte Carlo approach with the freedom to model the BHB and the other signals present in the data stream simultaneously, rather than trying to remove other signals and risk biasing the remaining data. The Mock LISA Data Challenge is a community of LISA scientists who generate rounds of simulated LISA noise with increasingly difficult signal content, and invite the LISA data analysis community to exercise their methods, or develop new methods, in an attempt to extract the parameters for the signals embedded in the mock data. In addition to practical approaches such ,is this to assess the level of parameter accuracy, one can apply the Fisher matrix formalism to assess both the statistical errors from noise and the theoretical errors

McWilliams, Sean T.↗

Axial-Flow Turbine Rotor Discharge-Flow Overexpansion and Limit-Loading Condition, Part I: Computational Fluid Dynamics (CFD) Investigation

A Computational Fluid Dynamic (CFD) investigation is conducted over a two-dimensional axial-flow turbine rotor blade row to study the phenomena of turbine rotor discharge flow overexpansion at subcritical, critical, and supercritical conditions. Quantitative data of the mean-flow Mach numbers, mean-flow angles, the tangential blade pressure forces, the mean-flow mass flux, and the flow-path total pressure loss coefficients, averaged or integrated across the two-dimensional computational domain encompassing two blade-passages, are obtained over a series of 14 inlet-total to exit-static pressure ratios, from 1.5 (un-choked; subcritical condition) to 10.0 (supercritical with excessively high pressure ratio.) Detailed flow features over the full domain-of-computation, such as the streamline patterns, Mach contours, pressure contours, blade surface pressure distributions, etc. are collected and displayed in this paper. A formal, quantitative definition of the limit loading condition based on the channel flow theory is proposed and explained. Contrary to the comments made in the historical works performed on this subject, about the deficiency of the theoretical methods applied in analyzing this phenomena, using modern CFD method for the study of this subject appears to be quite adequate and successful. This paper describes the CFD work and its findings.

Axial Flow↗

Extension of the Dytlewski-style dead time correction formalism for neutron multiplicity counting to any order

Here, neutron multiplicity counting using shift-register calculus is an established technique in the science of international nuclear safeguards for the identification, verification, and assay of special nuclear materials. Typically passive counting is used for Pu and mixed Pu-U items and active methods are used for U materials. Three measured counting rates, singles, doubles and triples are measured and, in combination with a simple analytical point-model, are used to calculate characteristics of the measurement item in terms of known detector and nuclear parameters. However, the measurement problem usually involves more than three quantities of interest, but even in cases where the next higher order count rate, quads, is statistically viable, it is not quantitatively applied because corrections for dead time losses are currently not available in the predominant analysis paradigm. In this work we overcome this limitation by extending the commonly used dead time correction method, developed by Dytlewski, to quads. We also give results for pents, which may be of interest for certain special investigations. Extension to still higher orders, may be accomplished by inspection based on the sequence presented. We discuss the foundations of the Dytlewski method, give limiting cases, and highlight the opportunities and implications that these new results expose. In particular there exist a number of ways in which the new results may be combined with other approaches to extract the correlated rates, and this leads to various practical implementations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

An efficient method for energy levels calculation using full symmetry and exact kinetic energy operator: Tetrahedral molecules

A simultaneous use of the full molecular symmetry and of an exact kinetic energy operator (KEO) is of key importance for accurate predictions of vibrational levels at a high energy range from a potential energy surface (PES). An efficient method that permits a fast convergence of variational calculations would allow iterative optimization of the PES parameters using experimental data. In this work, we propose such a method applied to tetrahedral AB{sub 4} molecules for which a use of high symmetry is crucial for vibrational calculations. A symmetry-adapted contracted angular basis set for six redundant angles is introduced. Simple formulas using this basis set for explicit calculation of the angular matrix elements of KEO and PES are reported. The symmetric form (six redundant angles) of vibrational KEO without the sin(q){sup −2} type singularity is derived. The efficient recursive algorithm based on the tensorial formalism is used for the calculation of vibrational matrix elements. A good basis set convergence for the calculations of vibrational levels of the CH{sub 4} molecule is demonstrated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficiently evaluating loop integrals in the EFTofLSS using QFT integrals with massive propagators

We develop a new way to analytically calculate loop integrals in the Effective Field Theory of Large Scale-Structure. Previous available methods show severe limitations beyond the one-loop power spectrum due to analytical challenges and computational and memory costs. Our new method is based on fitting the linear power spectrum with cosmology-independent functions that resemble integer powers of quantum field theory massive propagators with complex masses. A remarkably small number of them is sufficient to reach enough accuracy. Similarly to former approaches, the cosmology dependence is encoded in the coordinate vector of the expansion of the linear power spectrum in our basis. We first produce cosmology-independent tensors where each entry is the loop integral evaluated on a given combination of basis vectors. For each cosmology, the evaluation of a loop integral amounts to contracting this tensor with the coordinate vector of the linear power spectrum. The 3-dimensional loop integrals for our basis functions can be evaluated using techniques familiar to particle physics, such as recursion relations and Feynman parametrization. We apply our formalism to evaluate the one-loop bispectrum of galaxies in redshift space. The final analytical expressions are quite simple and can be evaluated with little computational and memory cost. We show that the same expressions resolve the integration of all one-loop N-point function in the EFTofLSS. This method, which is originally presented here, has already been applied in the first one-loop bispectrum analysis of the BOSS data to constraint ΛCDM parameters and primordial non-Gaussianities [1, 2].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Chapter 24: Strategic Energy Management (SEM) Evaluation Protocol. The Uniform Methods Project: Methods for Determining Energy Efficiency Savings for Specific Measures

Strategic energy management (SEM) focuses on achieving energy-efficiency improvements through systematic and planned changes in facility operations, maintenance, and behaviors (OM&B) and capital equipment upgrades in large energy-using facilities, including industrial buildings, commercial buildings, and multi-facility organizations such as campuses or communities. Facilities can institute a spectrum of SEM actions, ranging from a simple process for regularly identifying energy-savings actions, to establishing a formal, third-party recognized or certified SEM framework for continuous improvement of energy performance. In general, SEM programs that would be considered part of a utility program will contain a set of energy-reducing goals, principles, and practices emphasizing continuous improvements in energy performance or savings through energy management and an energy management system (EnMS).

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

On the Solution of the Continuity Equation for Precipitating Electrons in Solar Flares

Electrons accelerated in solar flares are injected into the surrounding plasma, where they are subjected to the influence of collisional (Coulomb) energy losses. Their evolution is modeled by a partial differential equation describing continuity of electron number. In a recent paper, Dobranskis & Zharkova claim to have found an "updated exact analytical solution" to this continuity equation. Their solution contains an additional term that drives an exponential decrease in electron density with depth, leading them to assert that the well-known solution derived by Brown, Syrovatskii & Shmeleva, and many others is invalid. We show that the solution of Dobranskis & Zharkova results from a fundamental error in the application of the method of characteristics and is hence incorrect. Further, their comparison of the "new" analytical solution with numerical solutions of the Fokker-Planck equation fails to lend support to their result.We conclude that Dobranskis & Zharkova's solution of the universally accepted and well-established continuity equation is incorrect, and that their criticism of the correct solution is unfounded. We also demonstrate the formal equivalence of the approaches of Syrovatskii & Shmeleva and Brown, with particular reference to the evolution of the electron flux and number density (both differential in energy) in a collisional thick target. We strongly urge use of these long-established, correct solutions in future works.

Sun: flares↗

Understanding the Awareness, Knowledge and Perceptions of Interventional Radiology Amongst Undergraduates in the UK

Introduction: The role of interventional radiology (IR) in health care is rapidly expanding in the United Kingdom (UK). However, the presence of IR in the current undergraduate curriculum is lacking. Purpose: This study investigates the awareness, knowledge and perceptions of IR in undergraduates through the clinical years of medical school to suggest future interventions. Methods: A total of 360 questionnaires, designed to test undergraduates’ knowledge of, and exposure to, IR were distributed over a 14-day period in a UK medical institution. Results: Of the 312 respondents (86.7%), the majority had never formally shadowed a radiologist (65.4%) or encountered teaching from an interventional radiologist (54.5%). Most students (88.5%) felt their knowledge of IR was lower than other specialties. Students’ knowledge of IR procedures did not statistically increase through the clinical years. Students were most cognisant of vascular procedures and least aware of IR’s role in oncology. Only 10.9% would consider a career in IR, with 39.7% citing the lack of knowledge and 18.6% citing patient contact as the top reasons for this. Conclusion: Current students have poor exposure to IR through the undergraduate curriculum. Consequently, there is no increase in knowledge of IR clinical practice through the clinical years. Multi-level interventions, incorporating universities and radiological societies are required to increase undergraduate awareness and promote IR as a career. Level of Evidence: Level 3, Local non-random sample.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Weather effects on the lifecycle of U.S. Department of Defense equipment replacement (WELDER)

Extreme weather has a direct and significant impact on buildings and infrastructure, resulting in billions of dollars of damage each year. This problem continues to grow as climate patterns change and buildings are exposed to new and different hazards than what they were designed to withstand. In order to better plan for the long-range sustainment, restoration, modernization, and eventual recapitalization of these buildings, organizations with large building portfolios, such as the U.S. Department of Defense (DoD), must have an awareness of the risks that these extreme weather events present. This research aimed to develop an approach to estimate condition loss and reduction in service life for the components of a building due to extreme weather hazards, to understand the risks that may be present in certain buildings and building systems. To achieve this objective, a damage association matrix was developed that categorizes climate hazards, the damage modes that they produce, and the individual component types impacted. This damage matrix formally links state-of-the-art climate model output, which provides projections of the probability of various climate hazards with a damage effects model that quantifies the consequence on component-level condition and service life. This method is applied to an actual portfolio of buildings in a particular geographic location and with a pre-defined component inventory that comprises the building. This approach can be aggregated to the system-, facility-, and site-level thus helping support billions of dollars in recapitalization decisions related to restoration/modernization of facilities.

54 ENVIRONMENTAL SCIENCES↗

Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface

A multilevel approach to sample the potential energy surface in a path integral formalism is proposed. The purpose is to reduce the required number of ab initio evaluations of energy and forces in ab initio path integral molecular dynamics (AI-PIMD) simulation, without compromising the overall accuracy. To validate the method, the internal energy and free energy of an Einstein crystal are calculated and compared with the analytical solutions. As a preliminary application, we assess the performance of the method in a realistic model—the FCC phase of dense atomic hydrogen, in which the calculated result shows that the acceleration rate is about 3 to 4-fold for a two-level implementation, and can be increased up to 10 times if extrapolation is used. With only 16 beads used for the ab initio potential sampling, this method gives a well converged internal energy. The residual error in pressure is just about 3 GPa, whereas it is about 20 GPa for a plain AI-PIMD calculation with the same number of beads. The vibrational free energy of the FCC phase of dense hydrogen at 300 K is also calculated with an AI-PIMD thermodynamic integration method, which gives a result of about 0.51 eV/proton at a density of r{sub s}=0.912.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Robust second-order approximation of the compressible Euler equations with an arbitrary equation of state

Here, this paper is concerned with the approximation of the compressible Euler equations supplemented with an arbitrary or tabulated equation of state. The proposed approximation technique is robust, formally second-order accurate in space, invariant-domain preserving, and works for every equation of state, tabulated or analytic, provided the pressure is nonnegative. An entropy surrogate functional that grows across shocks is proposed. The numerical method is verified with novel analytical solutions and then validated with several computational benchmarks seen in the literature including problems with composite waves.

97 MATHEMATICS AND COMPUTING↗

Temperature- and frequency-dependent dielectric properties of organic–inorganic hybrid compound: (C{sub 6}H{sub 9}N{sub 2}){sub 2}(Hg{sub 0.75}Cd{sub 0.25})Cl{sub 4}

Highlights: • We have synthesized a new hybrid compound of composition (C6H9N2)2(Hg0.75Cd0.25)Cl4. • The Ac conductivity of the title material was studied as a function of frequency and temperature. • The dielectric data have been analyzed in modulus formalism using KWW. - Abstract: The bis-2-amino-4-picolinium tetrachloromercurate-cadmate compound (C{sub 6}H{sub 9}N{sub 2}){sub 2}(Hg{sub 0.75}Cd{sub 0.25})Cl{sub 4} was prepared by hydrothermal method and characterized by X-ray diffraction (XRD) technique. The electrical properties of the compound were studied using impedance spectroscopy in the frequency and temperature range of 200 Hz–5 MHz and 308–403 K, respectively. The equivalent circuit is modeled by a combination of a parallel Rp//CPE circuit to explain the impedance results. The dielectric data were analyzed using complex electrical modulus M* at various temperatures. The activation energy responsible for the relaxation calculated from the modulus spectra is found to be almost the same as the value obtained from the temperature variation of dc conductivity. The electrical modulus and its scaling behavior are also investigated.

36 MATERIALS SCIENCE↗

Stability of a tokamak plasma with diffuse toroidal rotation

The present work considers the stability of a high- $\beta$ , large aspect ratio, circular plasma with diffuse profiles for the safety factor and the angular toroidal frequency (López & Guazzotto, Phys. Plasmas , vol. 24, 032501). An application of the Frieman–Rotenberg formalism results in a system of scalar eigenmode equations whose coupling is retained at the plasma–vacuum transition but is disregarded across the plasma column, which is a standard practice. The solution technique consists of a multidimensional shooting method for the poloidal harmonics; robust initial guesses are constructed by solving the dispersion relation in the static scenario with vanishing magnetic shear. Flow shear appears as a high- $\beta$ toroidal contribution, and we illustrate its destabilizing influence on $n=1$ external kink modes in the presence of ideal and resistive walls. Internal resonances are avoided by means of the selection of appropriate equilibrium parameters. The stabilizing influence of a finite positive average magnetic shear is also exemplified.

Physics↗

“How Should I Think about Voltage? What Is Overpotential?”: Establishing an Organic Chemistry Intuition for Electrochemistry

Redox reactions are ubiquitous in organic synthesis and intrinsic to organic electrosynthesis. The language and concepts used to describe reactions in these domains are sufficiently different to create barriers that hinder broader adoption and understanding of electrochemical methods. To bridge these gaps, this Synopsis compares chemical and electrochemical redox reactions, including concepts of free energy, voltage, kinetic barriers, and overpotential. Furthermore, this discussion is intended to increase the accessibility of electrochemistry for organic chemists lacking formal training in this area.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-Dimensional Electronic–Vibrational Spectroscopy of Coupled Molecular Complexes: A Near-Analytical Approach

This work presents theoretical calculations of the two-dimensional electronic–vibrational (2DEV) spectrum of a vibronically coupled molecular dimer using a near-analytical method. In strongly coupled dimers, where the IR mode is resonant with the electronic energy gap between the excitons, multiple infrared transitions become allowed that are forbidden in weakly coupled systems that have a nonresonant IR mode. This formalism enables the coherences and population contributions to be explored separately and allows efficient calculation of relaxation rates between the vibronic states. At short times, we find strong contributions of vibronic coherences to the 2DEV spectra. They decay fairly rapidly, giving rise to strong population signals. Although the interpretation of 2DEV spectra is considerably more complex than that for weakly coupled systems, the richness of the spectra and the necessity to consider both visible and infrared transition moments suggest that such analysis will be very valuable in characterizing the role of vibronic effects in ultrafast molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

LINKING TESTS OF GRAVITY ON ALL SCALES: FROM THE STRONG-FIELD REGIME TO COSMOLOGY

The current effort to test general relativity (GR) employs multiple disparate formalisms for different observables, obscuring the relations between laboratory, astrophysical, and cosmological constraints. To remedy this situation, we develop a parameter space for comparing tests of gravity on all scales in the universe. In particular, we present new methods for linking cosmological large-scale structure, the cosmic microwave background, and gravitational waves with classic PPN tests of gravity. Diagrams of this gravitational parameter space reveal a noticeable untested regime. The untested window, which separates small-scale systems from the troubled cosmological regime, could potentially hide the onset of corrections to GR.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-783: Using Matrixx to Determine Transit Dose Contribution Over Clinically Useful Limits of HDR Source Activity

Purpose: Most HDR brachytherapy treatment planning systems (TPS) use TG-43 formalism to calculate dose without including transit dose corrections. Historically, measurement of this contribution has required sophisticated apparatus unavailable in most hospitals. We use Matrixx to investigate several scenarios where transit dose contribution may effect a clinical treatment. Methods: Treatment plans were generated using Oncentra Brachy TPS (Version 4.3.0.410, Nucletron ) on a CT scan of a 24-catheter Freiburg applicator (Nucletron ) laid flat on the MatriXX (IBA) detector. This detector is an array of 1020 parallel plate ion chambers. All 24 catheters were digitized and dwells within a central square region of 5×5cm of the applicator were activated. Each of the active catheters had 6 dwells in increments of 1.0cm. The plans were normalized to 10mm. This places the 100% isodose line at the correct effective point of measurement, which lies half-way between the parallel plates of the ion chambers. It is also within the clinically relevant treatment depth for superficial applications. A total of 6 plans were delivered for 3 prescription doses, 1Gy, 2Gy and 4Gy using source activities of 2.9Ci and 11.2Ci. The MatriXX array was operated to capture dosimetric snaps every 500ms and yielded an integral dose at the end of treatment. Results: A comparison of integral dose from 2 different source activities shows that the transit dose contribution is larger when the source activity is higher. It is also observed that the relative transit dose contribution decreases as prescription dose increases. This is quantified by the Gamma analysis. Conclusion: We have demonstrated that the Matrixx detector can be used to evaluate the contribution for a HDR source during transit from the HDR afterloader to a dwell location, and between adjacent dwell locations.

60 APPLIED LIFE SCIENCES↗