Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational approach”

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,153 records · Page 64

Enhanced beam-beam modeling to include longitudinal variation during weak-strong simulation

Beam-beam interactions pose substantial challenges in the design and operation of circular colliders, significantly affecting their performance. In particular, the weak-strong simulation approach is pivotal for investigating single-particle dynamics during the collider design phase. This paper evaluates the limitations of existing models in weak-strong simulations, noting that while they accurately account for energy changes due to slingshot effects, they fail to incorporate longitudinal coordinate changes ( z variation). To address this gap, we introduce two novel transformations that enhance Hirata’s original framework by including both z variation and slingshot effect-induced energy changes. Through rigorous mathematical analysis and extensive weak-strong simulation studies, we validate the efficacy of these enhancements in achieving a more precise simulation of beam-beam interactions. Our results reveal that although z variation constitutes a higher-order effect and does not substantially affect the emittance growth rate within the specific design parameters of the Electron-Ion Collider, the refined model offers improved accuracy, particularly in scenarios involving the interaction between beam-beam effects and other random diffusion processes, as well as in simulations incorporating realistic lattice models. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evaluation of global terrestrial evapotranspiration using state-of-the-art approaches in remote sensing, machine learning and land surface modeling

Abstract. Evapotranspiration (ET) is critical in linking global water, carbon and energy cycles. However, direct measurement of global terrestrial ET is not feasible. Here, we first reviewed the basic theory and state-of-the-art approaches for estimating global terrestrial ET, including remote-sensing-based physical models, machine-learning algorithms and land surface models (LSMs). We then utilized 4 remote-sensing-based physical models, 2 machine-learning algorithms and 14 LSMs to analyze the spatial and temporal variations in global terrestrial ET. The results showed that the ensemble means of annual global terrestrial ET estimated by these three categories of approaches agreed well, with values ranging from 589.6 mm yr-1 (6.56×104 km3 yr-1) to 617.1 mm yr-1 (6.87×10 4 km3 yr -1 ). For the period from 1982 to 2011, both the ensembles of remote-sensing-based physical models and machine-learning algorithms suggested increasing trends in global terrestrial ET (0.62 mm yr -2 with a significance level of p<0.05 and 0.38 mm yr -2 with a significance level of p<0.05, respectively). In contrast, the ensemble mean of the LSMs showed no statistically significant change (0.23 mm yr -2 , p>0.05), although many of the individual LSMs reproduced an increasing trend. Nevertheless, all 20 models used in this study showed that anthropogenic Earth greening had a positive role in increasing terrestrial ET. The concurrent small interannual variability, i.e., relative stability, found in all estimates of global terrestrial ET, suggests that a potential planetary boundary exists in regulating global terrestrial ET, with the value of this boundary being around 600 mm yr -1 . Uncertainties among approaches were identified in specific regions, particularly in the Amazon Basin and arid/semiarid regions. Improvements in parameterizing water stress and canopy dynamics, the utilization of new available satellite retrievals and deep-learning methods, and model–data fusion will advance our predictive understanding of global terrestrial ET.

58 GEOSCIENCES↗

Evaluation of global terrestrial evapotranspiration using state-of-the-art approaches in remote sensing, machine learning and land surface modeling

Evapotranspiration (ET) is critical in linking global water, carbon and energy cycles. However, direct measurement of global terrestrial ET is not feasible. Here, we first reviewed the basic theory and state-of-the-art approaches for estimating global terrestrial ET, including remote-sensing-based physical models, machine-learning algorithms and land surface models (LSMs). We then utilized 4 remote-sensing-based physical models, 2 machine-learning algorithms and 14 LSMs to analyze the spatial and temporal variations in global terrestrial ET. The results showed that the ensemble means of annual global terrestrial ET estimated by these three categories of approaches agreed well, with values ranging from 589.6 mm/yr (6.56×10^4 cu.km/yr) to 617.1 mm/yr (6.87×10^4 cu.km/yr). For the period from 1982 to 2011, both the ensembles of remote-sensing-based physical models and machine-learning algorithms suggested increasing trends in global terrestrial ET (0.62 mm/sq.yr with a significance level of p<0.05 and 0.38 mm yr−2 with a significance level of p<0.05, respectively). In contrast, the ensemble mean of the LSMs showed no statistically significant change (0.23 mm/sq.yr, p>0.05), although many of the individual LSMs reproduced an increasing trend. Nevertheless, all 20 models used in this study showed that anthropogenic Earth greening had a positive role in increasing terrestrial ET. The concurrent small interannual variability, i.e., relative stability, found in all estimates of global terrestrial ET, suggests that a potential planetary boundary exists in regulating global terrestrial ET, with the value of this boundary being around 600 mm/yr. Uncertainties among approaches were identified in specific regions, particularly in the Amazon Basin and arid/semiarid regions. Improvements in parameterizing water stress and canopy dynamics, the utilization of new available satellite retrievals and deep-learning methods, and model–data fusion will advance our predictive understanding of global terrestrial ET.

surface modeling↗

A mathematical representation of an advanced helicopter for piloted simulator investigations of control system and display variations

A mathematical model of an advanced helicopter is described. The model is suitable for use in control/display research involving piloted simulation. The general design approach for the six degree of freedom equations of motion is to use the full set of nonlinear gravitational and inertial terms of the equations and to express the aerodynamic forces and moments as the reference values and first order terms of a Taylor series expansion about a reference trajectory defined as a function of longitudinal airspeed. Provisions for several different specific and generic flight control systems are included in the model. The logic required to drive various flight control and weapon delivery symbols on a pilot's electronic display is also provided. Finally, the model includes a simplified representation of low altitude wind and turbulence effects. This model was used in a piloted simulator investigation of the effects of control system and display variations for an attack helicopter mission.

Aiken, E. W.↗

Urban Air Quality Management at Low-Cost Using Micro Air Sensors: A Case Study From Accra, Ghana

Urban air quality management is dependent on the availability of local air pollution data. In many major urban centers of Africa, there is limited to non-existent information on air quality. This is gradually changing in part due to the increasing use of micro air sensors which have the potential to enable the generation of ground-based air quality data at fine scales for understanding local emission trends. Regional literature on the application of the high-resolution data for emission source identification in this region is limited. In this study a micro air sensor was co-located at the Physics Department, University of Ghana with a reference grade instrument to evaluate its performance for estimating PM2.5 pollution accurately at fine scales and the value of this data in identification of local sources and their behavior over time. For this study 15 weeks of data at hourly resolution with approximately 2500 data pairs are generated and analyzed (June 01, 2023, to September 15, 2023). For this time period a coefficient of determination (r 2 ) of 0.83 was generated with a mean absolute error (MAE) of 5.44 μgm -3 between the pre local calibration micro air sensor (i.e. out of box) and the reference-grade instrument. Following currently accepted best practice methods (see e.g., PAS4023) a domain specific (i.e. local) calibration factor was generated using a multi-linear regression model and when this factor is applied to the micro air sensor data, a reduction i.e., improvement in MAE to 1.43 μgm -3 was found. Daily variation was calculated, a receptor model was applied, and time series plots as a function of wind direction were generated, including PM2.5/PM10 ratio scatter and count plots to explore the utility of this observational approach for local source identification. The 3 data sets were compared (out of box, domain calibrated and reference-grade) and it was found that although there were variations in the data reported, source areas highlighted based on these data were similar, with input from local sources such as traffic emissions and biomass burning. As the temporal resolution of observational data associated with these micro air sensors is higher than for reference grade instruments (primarily due to costs and logistics limitations), they have the potential to provide insight into the complex, often hyper localized sources associated with urban areas, such as those found in major African cities.

Source apportionment↗

Design and Analysis of an Advanced Three-Point Bend Test Approach for Miniature Irradiated Disk Specimens

Candidate tungsten armor materials in a magnetic confinement fusion device must be able to withstand thermal variation that leads to internal stresses caused by the impinging heat load. In addition, the thermomechanical properties of these materials are degraded by irradiation-induced defect accumulation. Fission reactor–based irradiation data are used to predict the fusion neutron damage and property change. This study examines the motivation and design of a custom-designed three-point bend test for neutron-irradiated disk specimens that are 3 mm in diameter to be able to define the flexural strength of advanced tungsten materials, alloys, and composites—and to the extent that embrittlement occurs after neutron irradiation. The theory provided shows a calculation for the flexural deflection and shear deflection due to the small-geometry constraints. A finite element deformation analysis is performed to evaluate the mechanical stress field of disk bend specimens. The stress values above 80% of the maximum stress are concentrated in 2.4 mm of the 3.0-mm length of the centerline across the tungsten disk diameter. A bend test fixture has been designed and fabricated to enable testing of these specimens with precisely engineered tolerance and minimal machine compliance. This fixture will be able to be placed inside a universal testing frame at elevated temperatures for the mechanical property evaluation of future neutron-irradiated disk specimens.

36 MATERIALS SCIENCE↗

Direct prediction of saturated neoclassical tearing modes in slab using an equilibrium approach

We demonstrate for the first time that the nonlinear saturation of neoclassical tearing modes (NTMs) can be found directly using a variational principle based on Taylor relaxation, without needing to simulate the intermediate, resistivity-dependent dynamics. As in previous investigations of classical tearing mode saturation (Loizu et al 2020 Phys. Plasmas 27 070701; Loizu and Bonfiglio 2023 J. Plasma Phys. 89 905890507), we make use of Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502), an equilibrium solver based on the variational principle of the multi-region relaxed magnetohydrodynamics (MHDs), featuring stepped pressure profiles and arbitrary magnetic topology. We work in slab geometry and employ a simple bootstrap current model J bs = C$\boldsymbol{\nabla}$p to study the bootstrap-driven tearing modes, scanning over the asymptotic matching parameter Δ' and bootstrap current strength. Saturated island widths produced by SPEC agree well with the predictions of an initial value resistive MHDs code (Huang and Bhattacharjee 2016 Astrophys. J. 818 20) while being orders of magnitude faster to calculate. Additionally, we observe good agreement with a simple analytical modified Rutherford equation, without requiring any fitting coefficients. The match is obtained for both linearly unstable classical tearing modes in the presence of bootstrap current, and NTMs, which are linearly stable but nonlinear-unstable due to the effects of the bootstrap current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A unified approach to the analysis and design of elasto-plastic structures with mechanical contact

With structural design in mind, a new unified variational model has been developed which represents the mechanics of deformation elasto-plasticity with unilateral contact conditions. For a design problem formulated as maximization of the load carrying capacity of a structure under certain constraints, the unified model allows for a simultaneous analysis and design synthesis for a whole range of mechanical behavior.

Bendsoe, Martin P.↗

Atmospheric Rotational Effects on Mars Based on the NASA Ames General Circulation Model: Angular Momentum Approach

The objective of the investigation is to determine the motion of the rotational axis of Mars as a result of mass variations in the atmosphere and condensation and sublimation of CO2 ice on the polar caps. A planet experiences this type of motion if it has an atmosphere, which is changing its mass distribution with respect to the solid body of the planet and/or it is asymmetrically changing the amount of ice at the polar caps. The physical principle involved is the conservation of angular momentum, one can get a feeling for it by sitting on a well oiled swivel chair holding a rotating wheel on a horizontal direction and then changing the rotation axis of the wheel to a vertical direction. The person holding the wheel and the chair would begin to rotate in opposite direction to the rotation of the wheel. The motions of Mars atmosphere and the ice caps variations are obtained from a mathematical model developed at the NASA Ames Research Center. The model produces outputs for a time span of one Martian year, which is equivalent to 687 Earth days. The results indicate that Mars axis of rotation moves in a spiral with respect to a reference point on the surface of the planet. It can move as far away as 35.3 cm from the initial location as a result of both mass variations in the atmosphere and asymmetric ice variations at the polar caps. Furthermore the pole performs close to two revolutions around the reference point during a Martian year. This motion is a combination of two motions, one produced by the atmospheric mass variations and another due to the variations in the ice caps. The motion due to the atmospheric variations is a spiral performing about two and a half revolutions around the reference point during which the pole can move as far as 40.9 cm. The motion due to variations in the ice caps is a spiral performing almost three revolutions during which the pole can move as far as 32.8 cm.

Sanchez, Braulio V.↗

Robust Stability Analysis of the Space Launch System Control Design: A Singular Value Approach

Classical stability analysis consists of breaking the feedback loops one at a time and determining separately how much gain or phase variations would destabilize the stable nominal feedback system. For typical launch vehicle control design, classical control techniques are generally employed. In addition to stability margins, frequency domain Monte Carlo methods are used to evaluate the robustness of the design. However, such techniques were developed for Single-Input-Single-Output (SISO) systems and do not take into consideration the off-diagonal terms in the transfer function matrix of Multi-Input-Multi-Output (MIMO) systems. Robust stability analysis techniques such as H(sub infinity) and mu are applicable to MIMO systems but have not been adopted as standard practices within the launch vehicle controls community. This paper took advantage of a simple singular-value-based MIMO stability margin evaluation method based on work done by Mukhopadhyay and Newsom and applied it to the SLS high-fidelity dynamics model. The method computes a simultaneous multi-loop gain and phase margin that could be related back to classical margins. The results presented in this paper suggest that for the SLS system, traditional SISO stability margins are similar to the MIMO margins. This additional level of verification provides confidence in the robustness of the control design.

Pei, Jing↗

A statistical approach to screening isotopic signatures in monitoring for underground nuclear explosions

The ability to differentiate between atmospheric radionuclide signatures from underground nuclear explosions (UNEs) and signals from other sources, such as medical isotope-production facilities and nuclear reactors, can be critical to the detection and monitoring of unannounced, low-yield nuclear events. Signatures having anomalously high amplitudes, compared to background levels, remain the best indicator in screening for a UNE. However, isotopic composition can further validate a suspected UNE signature, but separation from any atmospheric background composition is first necessary. To date, evaluating the challenges of performing this separation has typically involved comparing an observed background with a highly idealized deterministic model of radioxenon signature production by a UNE that does not consider the influence of post-detonation chemical/physical processes in the detonation cavity or the subsequent gas transport mechanisms that can also affect the isotopic composition of the detected gas signature. In addition, purely deterministic models, as previously employed, overlook the uncertainty inherent in estimating critical parameters characterizing the UNE and its detonation environment. In this paper, we create detailed, multi-parameter models of radionuclide evolution using the widely accepted England and Rider post-detonation radionuclide decay-chain network coupled to detailed models simulating physical production and transport processes affecting the gas signature. Because these models are governed by uncertain parameters including barometric fluctuations, realistic ranges of variation for each of the parameters influencing isotopic composition are then defined. A Latin-Hypercube sampling approach is used to obtain a random distribution of isotopic production and gas transport results associated with a given value of each parameter. We apply these results to background histories of two stations, one providing 4-isotope background measurements and the other providing two-isotope measurements associated with the 2013 DPRK announced UNE.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Orbiter altitude at ALT interface based on ALSES (Approach and Landing Shuttle Engineering Simulation) analyses. Mission planning, mission analysis and software formulation

The details of the orbiter altitude attainable at the approach and landing tests (ALT) interface determined by the approach and landing shuttle engineering simulation (manned) were documented. The analysis culminated in the verification of the trends observed in a similar analysis performed previously on the space vehicle dynamics simulation (unmanned). Altitude variations attributable to pilot steering variability ranged between 492 ft higher to 383 ft lower. The requirement for this parametric analysis is first elaborated. The specifications, assumptions, and analytical approach used to determine the orbiter altitude at the ALT interface are then presented, followed by the results of the analytical approach and the conclusions and recommendations.

Colwell, C. L.↗

A two-dimensional, time-dependent, near-earth magnetotail

A kinetic approach is presented for obtaining slowly time-dependent, self-consistent plasma, and magnetic field configurations for this region, using anisotropic particle distributions. At present, a 2D configuration is used. The time variation of the plasma/field configuration is assumed to be driven by a slowly varying, externally imposed electric field (representing, for example, a change in solar wind conditions). The numerical approach involves solving a 2D nonlinear Poisson equation for the magnetic vector potential at each time step. The local current density is calculated from particle velocity distributions under the assumption of adiabatic particle motion. The obtained equations of state are compared with those of double-adiabatic theory. Numerical examples of self-consistent responses of the near-earth plasma and field to an externally imposed electric field pulse are presented.

Whipple, Elden↗

Elevated temperature crack growth

Alloy 718 crack growth experiments were conducted to assess the ability of the selected path-independent (P-I) integrals to describe the elevated temperature crack growth behavior. These tests were performed on single edge notch (SEN) specimens under displacement control with multiple extensometers to monitor the specimen and crack mouth opening displacement (CMOD). The displacements in these tests were sufficiently high to induce bulk cyclic inelastic deformation of the specimen. Under these conditions, the linear elastic fracture mechanics (LEFM) parameter K does not correlate the crack growth data. The experimentally measured displacement gradients at the end of specimen gage length were used as the boundary conditions in elastic-plastic finite element method (FEM) analyses. These analyses were performed with a node release approach using CYANIDE, a GEAE FEM code, which included a gap element which is capable of efficiently simulating crack closure. Excellent correlation was obtained between the experimentally measured and predicted variation of stress and CMOD with crack length and the stress-CMOD loops for Alloy 718 tests conducted at 538 C. This confirmed the accuracy of the FEM crack growth simulation approach. The experimentally measured crack growth rate data correlated well the selected P-I integrals. These investigations have produced significant progress in developing P-I integrals as non-linear fracture mechanics parameters. The results suggest that this methodology has the potential of accurately describing elevated temperature crack growth behavior under the combined influence of thermal cycling and bulk elastic-inelastic deformation states.

Kim, K. S.↗

Remote Sensing of Crystal Shapes in Ice Clouds

Ice crystals in clouds exist in a virtually limitless variation of geometries. The most basic shapes of ice crystals are columnar or plate-like hexagonal prisms with aspect ratios determined by relative humidity and temperature. However, crystals in ice clouds generally display more complex structures owing to aggregation, riming and growth histories through varying temperature and humidity regimes. Crystal shape is relevant for cloud evolution as it affects microphysical properties such as fall speeds and aggregation efficiency. Furthermore, the scattering properties of ice crystals are affected by their general shape, as well as by microscopic features such as surface roughness, impurities and internal structure. To improve the representation of ice clouds in climate models, increased understanding of the global variation of crystal shape and how it relates to, e.g., location, cloud temperature and atmospheric state is crucial. Here, the remote sensing of ice crystal macroscale and microscale structure from airborne and space-based lidar depolarization observations and multi-directional measurements of total and polarized reflectances is reviewed. In addition, a brief overview is given of in situ and laboratory observations of ice crystal shape as well as the optical properties of ice crystals that serve as foundations for the remote sensing approaches. Lidar depolarization is generally found to increase with increasing cloud height and to vary with latitude. Although this variation is generally linked to the variation of ice crystal shape, the interpretation of the depolarization remains largely qualitative and more research is needed before quantitative conclusions about ice shape can be deduced. The angular variation of total and polarized reflectances of ice clouds has been analyzed by numerous studies in order to infer information about ice crystal shapes from them. From these studies it is apparent that pristine crystals with smooth surfaces are generally inconsistent with the data and thus crystal impurity, distortion or surface roughness is prevalent. However, conclusions about the dominating ice shapes are often inconclusive and contradictory and are highly dependent on the limited selection of shapes included in the investigations. Since ice crystal optical properties are mostly determined by the aspect ratios of the crystal components and their microscale structure, it is advised that remote sensing applications focus on the variation of these ice shape characteristics, rather than on the macroscale shape or habit. Recent studies use databases with nearly continuous ranges of crystal component aspect ratio and-or roughness levels to infer the variation of ice crystal shape from satellite and airborne remote sensing measurements. Here, the rationale and results of varying strategies for the remote sensing of ice crystal shape are reviewed. Observed systematic variations of ice crystal geometry with location, cloud height and atmospheric state suggested by the data are discussed. Finally, a prospective is given on the future of the remote sensing of ice cloud particle shapes.

clouds↗

Error mitigation in variational quantum eigensolvers using tailored probabilistic machine learning

Quantum computing technology has the potential to revolutionize the simulation of materials and molecules in the near future. A primary challenge in achieving near-term quantum advantage is effectively mitigating the noise effects inherent in current quantum processing units (QPUs). This challenge is also decisive in the context of quantum-classical hybrid schemes employing variational quantum eigensolvers (VQEs) that have attracted significant interest in recent years. In this paper, we present a method that employs parametric Gaussian process regression (GPR) within an active learning framework to mitigate noise in quantum computations, focusing on VQEs. Our approach, grounded in probabilistic machine learning, exploits a custom prior based on the VQE ansatz to capture the underlying correlations between VQE outputs for different variational parameters, thereby enhancing both accuracy and efficiency. We demonstrate the effectiveness of our method on a two-site Anderson impurity model and a eight-site Heisenberg model, using the IBM open-source quantum computing framework, Qiskit, showcasing substantial improvements in the accuracy of VQE outputs while reducing the number of direct QPU energy evaluations. This paper contributes to the ongoing efforts in quantum-error mitigation and optimization, bringing us a step closer to realizing the potential of quantum computing in quantum matter simulations. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗

Automated segmentation of soft X-ray tomography: Native cellular structure with submicron resolution at high-throughput for whole-cell quantitative imaging in yeast

Soft X-ray tomography (SXT) is an invaluable tool for quantitatively analyzing cellular structures at suboptical isotropic resolution. However, it has traditionally depended on manual segmentation, limiting its scalability for large datasets. Here, we leverage a deep learning-based autosegmentation pipeline to segment and label cellular structures in hundreds of cells across three Saccharomyces cerevisiae strains. This task-based pipeline uses manual iterative refinement to improve segmentation accuracy for key structures, including the cell body, nucleus, vacuole, and lipid droplets, enabling high-throughput and precise phenotypic analysis. Using this approach, we quantitatively compared the three-dimensional (3D) whole-cell morphometric characteristics of wild-type, VPH1-GFP, and vac14 strains, uncovering detailed strain-specific cell and organelle size and shape variations. We show the utility of SXT data for precise 3D curvature analysis of entire organelles and cells and detection of fine morphological features using surface meshes. Our approach facilitates comparative analyses with high spatial precision and statistical throughput, uncovering subtle morphological features at the single-cell and population level. This workflow significantly enhances our ability to characterize cell anatomy and supports scalable studies on the mesoscale, with applications in investigating cellular architecture, organelle biology, and genetic research across diverse biological contexts.

Chen, Jianhua [Lawrence Berkeley National Laborato↗

Otolith 87 Sr / 86 Sr identifies natal origin, movement and life history of Burbot Lota lota in the Kootenai River following 45 years of impoundment

Abstract Hydropower has long been an effective form of renewable energy, but the development of the resource has negatively impacted fishes throughout the world. Formulating recovery plans to mitigate these effects requires a firm understanding of species‐specific provenance, life history and life stage habitat requirements. However, the importance of this information is often only recognised after a substantial population decline, when obtaining it becomes demographically challenging. We provide a case study for this scenario and present a post hoc approach to gain insight into these critical data gaps for Burbot Lota lota in the upper Kootenai River in Montana and British Columbia. The population once supported a robust tribal and recreational harvest but declined following impoundment and these fisheries have remained closed or curtailed for nearly 20 years. The reasons for the decline are unknown and actions to assist recovery have been constrained by limited information regarding natal origin, life history and habitat use. We addressed this issue using strontium isotope ratios ( 87 Sr/ 86 Sr) in otoliths collected over four decades in Lake Koocanusa, the reservoir on the Kootenai River above Libby Dam. Due to the geologic heterogeneity of the basin and resulting variation in water 87 Sr/ 86 Sr, the assignment probability for Burbot to their region of natal origin approached unity. Otolith 87 Sr/ 86 Sr also indicated the population was comprised of mainstem‐ and tributary‐origin forms and their relative abundance changed over time. Mainstem Burbot was most abundant during the 1970s and 2000s, whereas tributary forms dominated between these periods. We suggest ecosystem changes following impoundment contributed to this shift in life history forms and the overall decline of the population. Our findings underscore the need for additional investigations of these impacts to guide restoration efforts.

Dunnigan, James L.↗