Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physics Regularization”

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 595 records · Page 33

Nonlocal chiral contributions to generalized parton distributions of the proton at nonzero skewness

We compute the one-loop contributions to spin-averaged generalized parton distributions (GPDs) in the proton from pseudoscalar mesons with intermediate octet and decuplet baryon states at nonzero skewness. Our framework is based on nonlocal covariant chiral effective theory, with ultraviolet divergences regularized by introducing a relativistic regulator derived consistently from the nonlocal Lagrangian. Using the splitting functions calculated from the nonlocal Lagrangian, we find the nonzero skewness GPDs from meson loops by convoluting with the phenomenological pion GPD and the generalized distribution amplitude, and verify that these satisfy the correct polynomiality properties. We also compute the lowest two moments of GPDs to quantify the meson loop effects on the Dirac, Pauli, and gravitational form factors of the proton.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An Evaluation of Antarctica as a Calibration Target for Passive Microwave Satellite Missions

Passive microwave remote sensing at L-band (1.4 GHz) is sensitive to soil moisture and sea surface salinity, both important climate variables. Science studies involving these variables can now take advantage of new satellite L-band observations. The first mission with regular global passive microwave observations at L-band is the European Space Agency's Soil Moisture and Ocean Salinity (SMOS), launched November, 2009. A second mission, NASA's Aquarius, was launched June, 201l. A third mission, NASA's Soil Moisture Active Passive (SMAP) is scheduled to launch in 2014. Together, these three missions may provide a decade-long data record -- provided that they are intercalibrated. The intercalibration is best performed at the radiance (brightness temperature) level, and Antarctica is proving to be a key calibration target. However, Antarctica has thus far not been fully characterized as a potential target. This paper will present evaluations of Antarctica as a microwave calibration target for the above satellite missions. Preliminary analyses have identified likely target areas, such as the vicinity of Dome-C and larger areas within East Antarctica. Physical sources of temporal and spatial variability of polar firn are key to assessing calibration uncertainty. These sources include spatial variability of accumulation rate, compaction, surface characteristics (dunes, micro-topography), wind patterns, and vertical profiles of density and temperature. Using primarily SMOS data, variability is being empirically characterized and attempts are being made to attribute observed variability to physical sources. One expected outcome of these studies is the potential discovery of techniques for remotely sensing--over all of Antarctica--parameters such as surface temperature.

Kim, Edward↗

QED medium effects in (anti)neutrino-nucleus and electron-nucleus scattering: Elastic scattering on nucleons

Interpretation of current and future neutrino oscillation and electron scattering experiments requires knowledge of lepton-nucleon and lepton-nucleus interactions at the percent level. We study the exchange of photons between charged particles and the nuclear medium for (anti)neutrino-, electron-, and muon-induced reactions inside a large nucleus. While quantum electrodynamics (QED)-medium contributions are formally suppressed by two powers of the electromagnetic coupling constant α when compared to the leading-order cross sections, low-energy modes and the nuclear size enhance the effect by orders of magnitude. They require a proper infrared regularization, which we implement as a screening of the electromagnetic interactions at atomic length scales or above. We provide approximate analytic expressions for the distortion of (anti)neutrino-nucleus and charged lepton-nucleus cross sections and evaluate the QED-medium effects for realistic values of the screening scale on the example of elastic scattering with nucleons inside the nucleus. We find new permille- to percent-level effects, which were not considered in either (anti)neutrino-nucleus or electron-nucleus scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Studying excited-state-specific perturbation theory on the Thiel set

Here, we explore the performance of a recently introduced N 5 -scaling excited-state-specific second order perturbation theory (ESMP2) on the singlet excitations of the Thiel benchmarking set. We find that, without regularization, ESMP2 is quite sensitive to π system size, performing well in molecules with small π systems but poorly in those with larger π systems. With regularization, ESMP2 is far less sensitive to π system size and shows a higher overall accuracy on the Thiel set than CC2, equation of motion-coupled cluster with singles and doubles, CC3, and a wide variety of time-dependent density functional approaches. Unsurprisingly, even regularized ESMP2 is less accurate than multi-reference perturbation theory on this test set, which can, in part, be explained by the set’s inclusion of some doubly excited states but none of the strong charge transfer states that often pose challenges for state-averaging. Beyond energetics, we find that the ESMP2 doubles norm offers a relatively low-cost way to test for doubly excited character without the need to define an active space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

LaserNetUS. Final report

In this proposal we have formed a consortium in response to a request from the DoE Office of Fusion Energy Sciences to establish a national network of Petawatt-class laser systems to act as domestic user facilities that can enable a broad range of frontier scientific research. This new entity is called LaserNet US. The network formed here is directly responsive to recommendations made in a recently released National Academy of Sciences Report with regard to US strategy for high intensity laser research, “Opportunities in Intense Ultrafast Lasers: Reaching for the Brightest Light”. As detailed in this report, the research in this area has the potential to transform science in a number of research fields and to open up new areas of fundamental research. High field science was initiated in the US in the 1990’s, and subsequently Europe and Asia have embraced this research field, investing more than a billion dollars in this area over the past several years. Our network includes five academic and two national lab-based high intensity laser facilities. These facilities are distributed geographically throughout the US and the network provides access to complementary facilities which have unique world-class laser and experimental capabilities. LaserNet US will make these existing Petawatt-class laser facilities available to users from around the country who until now have not had regular access to such machines. Consequently, the network, which will develop over time, will also become a key element in driving forward national research in high field and high energy density plasma science. In particular, the University of Michigan provides access to the Hercules laser facility and the T-cubed laser facility for collaborative experiments as part of the first year operation. In the second year - to begin in summer 2019 - the experiments will be chosen via proposal submission and review by an outside panel formed by DOE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synthesis, Crystal and Electronic Structure of the New Ternary Compound Ca3InAs3

Crystals of a new ternary compound in the Ca-In-As family, Ca3InAs3, have been successfully synthesized via flux growth techniques. This is only the third known compound between the respective elements. As elucidated by single-crystal X-ray diffraction measurements, Ca3InAs3 crystallizes in the orthorhombic space group Pnma (No. 62, Pearson symbol oP28) with unit cell parameters a = 12.296(2) Å, b = 4.2553(7) Å, and c = 13.735(2) Å. The smallest building motifs of the structure are InAs4 tetrahedra, which are connected to one another by shared As corners, forming infinite [InAs2As2/2] chains. The latter propagate along the crystallographic b-axis. The As-In-As bond angles within the InAs4 tetrahedra deviate from the ideal 109.5° value and range from 98.12(2)° to 116.53(2)°, attesting to a small distortion from the regular tetrahedral geometry. Electronic structure calculations indicate the opening of a bandgap, consistent with the expected (Ca2+)3(In3+)(As3–)3 formula breakdown based on conventional oxidation numbers. The calculations also show that the Ca–As interactions are an intermediate between covalent and ionic, while providing evidence of strong covalent features of the In–As interactions. Weak s-p hybridization of In states was observed, supporting the experimentally found deviation of the InAs4 moiety from the ideal tetrahedral symmetry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interpolation of compound semiconductor alloy parameters from those of their constituents

Several methods have been proposed for interpolation of the value of physical parameters of quaternary alloys from those of their constituent ternary and binary sub-alloys. These expressions agree when non-linear bowing terms are not required; they differ in how the bowing terms of the bounding ternaries should be utilized. Common interpolation expressions for quaternaries can be generalized into two groups: (1) those that use a linear interpolation of the nearest ternary parameter values and (2) those that interpolate over binary values with a bowing term derived from the bounding ternaries. The second group of methods is equivalent to a polynomial expansion over the alloy’s interpolation space. For compound semiconductor alloys, the geometry of the composition space is the direct sum of the group-III and group-V mixture sub-spaces. The mixture sub-spaces are best described using barycentric coordinates on a regular simplex. A general polynomial expansion of the value of an alloy parameter using barycentric coordinates for the group-III and group-V simplex spaces is described along with an algorithm to generate interpolation expressions for alloys with arbitrary numbers of elements, including quinary and senary alloys. It is shown that a polynomial expansion produces values in closer agreement with the direct gap of quaternaries lattice-matched to common substrates than do approaches using an interpolation of the ternary values, despite a prominent recommendation to the contrary. Finally, a quaternary correction term is described that improves the predicted direct bandgap energies of GaInAsSb for compositions near those lattice matched to InP, InAs, and GaSb.

Olesberg, Jonathon T. [Sandia National Laboratorie↗

The ATLAS experiment software on ARM

With an increased dataset obtained during the Run 3 of the LHC at CERN and the even larger expected increase of the dataset by more than one order of magnitude for the HL-LHC, the ATLAS experiment is reaching the limits of the current data processing model in terms of traditional CPU resources based on x86_64 architectures and an extensive program for software upgrades towards the HL-LHC has been set up. The ARM architecture is becoming a competitive and energy efficient alternative. Some surveys indicate its increased presence in HPCs and commercial clouds, and some WLCG sites have expressed their interest. Chip makers are also developing their next generation solutions on ARM architectures, sometimes combining ARM and GPU processors in the same chip. Consequently it is important that the ATLAS software embraces the change and is able to successfully exploit this architecture. We report on the successful porting to ARM of the Athena software framework, which is used by ATLAS for both online and offline computing operations. Furthermore we report on the successful validation of simulation workflows running on ARM resources. For this we have set up an ATLAS Grid site using ARM compatible middleware and containers on Amazon Web Services (AWS) ARM resources. The ARM version of Athena is fully integrated in the regular software build system and distributed in the same way as other software releases. In addition, the workflows have been integrated into the HEPscore benchmark suite which is the planned WLCG wide replacement of the HepSpec06 benchmark used for Grid site pledges. In the overall porting process we have used resources on AWS, Google Cloud Platform (GCP) and CERN. A performance comparison of different architectures and resources will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The CMOS Pseudo-Thyristor: a zero-static current discriminator circuit

We report a very low power discriminator circuit for pixelized detectors, called the Pseudo-Thyristor is described in this document. It is a positive feedback topology using regular PMOS and NMOS field-effect transistors (FET's) with zero static current. When a small charge is injected into the circuit, it flips rapidly due to the positive feedback and outputs a logic transition for further digitization. Simulation shows that in a 65 nm process, it is possible to achieve a detecting threshold below 5 fC while maintain the average power consumption below 10 micro-Watts when the hit occupancy is <10% for 40 MHz operation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An Evaluation of Antarctica as a Calibration Target for Passive Microwave Satellite Missions with Climate Data Record Applications

Passive microwave remote sensing at L-band (1.4 GHz) is sensitive to soil moisture and sea surface salinity, both important climate variables. Science studies involving these variables can now take advantage of new satellite L-band observations. The first mission with regular global passive microwave observations at L-band is the European Space Agency's Soil Moisture and Ocean Salinity (SMOS), launched November, 2009. A second mission, NASA's Aquarius, was launched June, 201 I. A third mission, NASA's Soil Moisture Active Passive (SMAP) is scheduled to launch in 2014. Together, these three missions may provide a decade-long data record-provided that they are intercalibrated. The intercalibration is best performed at the radiance (brightness temperature) level, and Antarctica is proving to be a key calibration target. However, Antarctica has thus far not been fully characterized as a potential target. This paper will present evaluations of Antarctica as a microwave calibration target for the above satellite missions. Preliminary analyses have identified likely target areas, such as the vicinity of Dome-C and larger areas within East Antarctica. Physical sources of temporal and spatial variability of polar firn are key to assessing calibration uncertainty. These sources include spatial variability of accumulation rate, compaction, surface characteristics (dunes, micro-topography), wind patterns, and vertical profiles of density and temperature. Using primarily SMOS data, variability is being empirically characterized and attempts are being made to attribute observed variability to physical sources. One expected outcome of these studies is the potential discovery of techniques for remotely sensing--over all of Antarctica-parameters such as surface temperature.

Kim, Edward↗

In-Canal Assay of High Specific Activity 60 Co at the Advanced Test Reactor

Production of high specific activity (HSA) 60 Co has recently resumed at the Advanced Test Reactor (ATR) at the Idaho National Laboratory. The technical steps of performing in-canal assay of HSA 60 Co targets at the ATR are described herein. The HSA targets are assayed on a regular basis, between cycles at the ATR, in order to assess the progress of activation. The targets are also assayed at the conclusion of the irradiation in order to provide activity estimates for the distributor and to be used for safety and shipping evaluations. These target assay activities must take place in the ATR canal to provide sufficient radiological shielding. Furthermore, a specialized assay fixture is used in conjunction with custom 60 Co radiation standards to assess the irradiated target. The specific activity of each irradiated target was determined by first measuring the 60 Co standards to determine a sensitivity factor. Then, each irradiated target was measured, and the integrated measurement was multiplied by the sensitivity factor to determine the total activity of the irradiated target. Finally, a correction factor was determined to adjust the reported activity, accounting for differences in the physical geometries of the standards and targets. These methods were used to accurately assay the gross activity of nine irradiated 60 Co targets, two of which were delivered to the distributor where confirmatory hot-cell assays were performed verifying the accuracy of the in-canal assay method.

07 ISOTOPE AND RADIATION SOURCES↗

Density functional based methods for light harvesting molecules

The goal of the proposal was to computationally study the photoinduced charge transfer excited states of organic molecules. The systems of interest are the large organic supramolecules that have potential for photovoltaic applications. The simulation of the charge transfer excited states in large supramolecular assemblies require quantum mechanical methods that can be applied to systems with hundreds of atoms. Through this project we have implemented a perturbative delta-self-consistent field method that can easily describe the charge transfer excited states with regular density functionals without range separation. The cost of the calculations of an excited state is similar to that of the ground state. The method also allows for calculation of the energy gradients of the excited state. The developed method was applied to study a number of different supramolecular systems containing hundreds of atoms with excellent agreement with experiment. The absorption characteristics and the influence of density functional approximations were also examined. These led to several studies on the variety of organic electron donor and acceptor molecules. Apart from these developments through this project we also improved the code capabilities by implementing new meta-GGA functionals. This was a major step since the code is purely density functional and earlier the functional capabilities were limited only to generalized gradient approximations. A total of 23 peer-reviewed publications in top chemical physics journals resulted from this grant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pulse-Timing Studies of X Ray Pulsars

The pulse-timing projects supported by NASA Grant NAG8-695 were motivated in large part by our long-standing interest in the physics of rotating neutron stars and particularly the variations in rotation due to internal and external fluctuating torques. For accretion-powered pulsars, observed as compact galactic X-ray sources, our work has been motivated by questions regarding the physics of matter accretion, mass transfer, and mass loss in these X-ray binary systems. For rotation-powered pulsars using X-ray, optical, and radio observations, we have examined the internal structure of neutron stars and mechanisms for secular spin-down. These issues are still central to the continuing process of understanding the complex behavior of these fascinating systems. The work supported by this grant is based mainly on Ginga observations of three pulsating X-ray sources. We observed two of these sources (Her X-1 and SMC X-1) within the framework of a NASA-ISAS program for U.S.-Japan collaborations to obtain and analyze Ginga observations. In addition, we joined with members of the Ginga team to apply our pulse-timing methods to PSR 0540-69 data obtained as a collateral benefit of the regular monitoring by Ginga of SN 1987A. In addition to these Ginga observations, there exist other relevant data on all three of these sources. The direction of our investigations has been affected by the necessity of including these supporting data in our analyses, and setting up the necessary collaborations has sometimes entailed extra work that was not anticipated in our proposals. Indeed, we have often taken the initiative in establishing these joint projects.

Deeter, John E.↗

Physically-based parameterization of spatially variable soil and vegetation using satellite multispectral data

A stochastic-geometric landsurface reflectance model is formulated and tested for the parameterization of spatially variable vegetation and soil at subpixel scales using satellite multispectral images without ground truth. Landscapes are conceptualized as 3-D Lambertian reflecting surfaces consisting of plant canopies, represented by solid geometric figures, superposed on a flat soil background. A computer simulation program is developed to investigate image characteristics at various spatial aggregations representative of satellite observational scales, or pixels. The evolution of the shape and structure of the red-infrared space, or scattergram, of typical semivegetated scenes is investigated by sequentially introducing model variables into the simulation. The analytical moments of the total pixel reflectance, including the mean, variance, spatial covariance, and cross-spectral covariance, are derived in terms of the moments of the individual fractional cover and reflectance components. The moments are applied to the solution of the inverse problem: The estimation of subpixel landscape properties on a pixel-by-pixel basis, given only one multispectral image and limited assumptions on the structure of the landscape. The landsurface reflectance model and inversion technique are tested using actual aerial radiometric data collected over regularly spaced pecan trees, and using both aerial and LANDSAT Thematic Mapper data obtained over discontinuous, randomly spaced conifer canopies in a natural forested watershed. Different amounts of solar backscattered diffuse radiation are assumed and the sensitivity of the estimated landsurface parameters to those amounts is examined.

Jasinski, Michael F.↗

A mountain-induced moist baroclinic wave test case for the dynamical cores of atmospheric general circulation models

Idealized test cases for the dynamical cores of atmospheric general circulation models are informative tools to assess the accuracy of the numerical designs and investigate the general characteristics of atmospheric motions. A new test case is introduced that is built upon a baroclinically unstable base state with an added orographic barrier. The topography is analytically prescribed and acts as a trigger of both baroclinic Rossby waves and inertia–gravity waves on a rotating, regular-sized planet. Both dry and idealized moist configurations are suggested. The latter utilizes the Kessler warm-rain precipitation scheme. The test case enhances the complexity of the existing test suite hierarchy and focuses on the impacts of two midlatitudinal mountain ridges on the circulation. Selected simulation examples from four dynamical cores are shown. These are the Spectral Element and Finite Volume dynamical cores, which are part of the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM), versions 2.1.3 and 2.2, and the Cubed-Sphere Finite Volume dynamical cores, which is new to CESM version 2.2. In addition, the Model for Prediction Across Scales (MPAS) is tested. The overall flow patterns agree well in the four dynamical cores, but the details can vary greatly. The examples highlight the broad palette of use cases for the test case and reveal physics–dynamics coupling issues.

58 GEOSCIENCES↗

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗

Topological terms with qubit regularization and relativistic quantum circuits

Qubit regularization provides a rich framework to explore quantum field theories. The freedom to choose how the important symmetries of the theory are embedded in the qubit regularization scheme allows us to construct new lattice models with rich phase diagrams. Some of the phases can contain topological terms which lead to critical phases. In this work we introduce and study the SU(3)-F qubit regularization scheme to embed the SO(3) spin-symmetry. We argue that qubit models in this regularization scheme contain several phases including a critical phase which describes the k = 1 Wess-Zumino-Witten (WZW) conformal field theory (CFT) at long distances, and two massive phases one of which is trvially gapped and the other which breaks the lattice translation symmetry. We construct a simple space-time Euclidean lattice model with a single coupling U and study it using the Monte Carlo method. We show the model has a critical phase at small U and a trivially massive phase at large U with a first order transition separating the two. Another feature of our model is that it is symmetric under space-time rotations, which means the temporal and spatial lattice spacing are connected to each other. The unitary time evolution operator obtained by a Wick rotation of the transfer matrix of our model can help us compute the physics of the k = 1 WZW CFT in real time without the need for tuning the temporal lattice spacing to zero. We use this idea to introduce the concept of a relativistic quantum circuit on a discrete space-time lattice.

Bhattacharya, Tanmoy↗