Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter space exploration”

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 361 records · Page 20

Optimizing Muscle Parameters in Musculoskeletal Modeling Using Monte Carlo Simulations

Astronauts assigned to long-duration missions experience bone and muscle atrophy in the lower limbs. The use of musculoskeletal simulation software has become a useful tool for modeling joint and muscle forces during human activity in reduced gravity as access to direct experimentation is limited. Knowledge of muscle and joint loads can better inform the design of exercise protocols and exercise countermeasure equipment. In this study, the LifeModeler(TM) (San Clemente, CA) biomechanics simulation software was used to model a squat exercise. The initial model using default parameters yielded physiologically reasonable hip-joint forces. However, no activation was predicted in some large muscles such as rectus femoris, which have been shown to be active in 1-g performance of the activity. Parametric testing was conducted using Monte Carlo methods and combinatorial reduction to find a muscle parameter set that more closely matched physiologically observed activation patterns during the squat exercise. Peak hip joint force using the default parameters was 2.96 times body weight (BW) and increased to 3.21 BW in an optimized, feature-selected test case. The rectus femoris was predicted to peak at 60.1% activation following muscle recruitment optimization, compared to 19.2% activation with default parameters. These results indicate the critical role that muscle parameters play in joint force estimation and the need for exploration of the solution space to achieve physiologically realistic muscle activation.

Hanson, Andrea↗

Decoding the Mechanisms of Phase Transitions from In Situ Microscopy Observations

Abstract Analysis of the temperature‐ and stimulus‐dependent imaging data toward elucidation of the physical transformations is an ubiquitous problem in multiple fields. Here, temperature‐induced phase transition in BaTiO 3 is explored using the machine learning analysis of domain morphologies visualized via variable‐temperature scanning transmission electron microscopy (STEM) imaging data. This approach is based on the multivariate statistical analysis of the time or temperature dependence of the statistical descriptors of the system, derived in turn from the categorical classification of observed domain structures or projection on the continuous parameter space of the feature extraction‐dimensionality reduction transform. The proposed workflow offers a powerful tool for the exploration of the dynamic data based on the statistics of image representation as a function of the external control variable to visualize the transformation pathways during phase transitions and chemical reactions. This can include the mesoscopic STEM data as demonstrated here, but also optical, chemical imaging, etc., data. It can further be extended to the higher dimensional spaces, for example, analysis of the combinatorial libraries of materials compositions.

Valleti, Sai Mani Prudhvi↗

Quasi-Wireless Capacitive Power Transfer with Secure Data Acquisition for Robotic Systems in Space Infrastructure

Space exploration is dependent on robotic systems that utilize end-effectors to collect samples, probe surfaces, and manipulate objects. These systems can rarely be designed to do all three, forcing engineers to make tradeoffs based on the mission parameters - i.e. should the robotic appendage have a claw, drill, or shovel, and which would be best suited for the mission? Additionally, as more industrial and government entities partake in space exploration, data protection is needed in transit and at rest. To address these challenges, we present a first-of-its- kind robotic linkage that has no wiring between the joints. Instead, quasi-wireless capacitive (QWiC) power transfer is used to send energy over the robot’s chassis without a return wire. This enables the system to be completely modular through the use of single-contact permanent magnet connections, allowing rapid alterations in joint kinematics and/or the changing of end-effectors. For collecting sensor data from the robotic arm and to send remote commands to it, we use a Supervisory Control and Data Acquisition (SCADA) system. Data transmission relies on MQTT and OPC UA communication protocols with encryption. The SCADA server logs and archives sensor data and provides the functionality for authorized users to send remote commands from SCADA client(s) to motors. A SCADA client can be any of the web browsers that connects to a server via a secure communication channel using SSL protocol. Furthermore, as an extra data protection mechanism, we inject noise to the sensor data traffic, which obfuscates the timing of sensor data packets and adds confusion about which data packet represents which motor.

wireless sensor networks↗

Simplifying the creation of thermal decomposition mechanisms for designing mixed fibre composites

Among the most challenging aspects of simulating thermal decomposition of fibre reinforced polymers is the determination of appropriate reaction parameters. Thermogravimetric analysis is typically used to generate decomposition data, to which the reaction parameters are then fit. When designing mixed fibre materials (e.g. combinations of glass and carbon fibres), the number of TGA experiments needed to explore the entire design space may be intractable. Here, we demonstrate the creation of two candidate proxy mixed fibre mechanisms and compare them to a mechanism created by fitting parameters from TGA on the mixed fibre composite. These mechanisms are then demonstrated in a 2D axisymmetric numerical decomposition, heat transfer, and porous flow model. We find a maximum 11% uncertainty in mass and 4% in temperature difference when using a proxy mechanism.

Scott, Sarah N. [Sandia National Lab. (SNL-CA), Li↗

O 2 dosage as a descriptor of TWC performance under lean/rich dithering in stoichiometric natural gas engines

The effects of oscillations in gas composition, known as lean/rich dithering, on the performance of a commercial Pd-based three-way catalyst (TWC) for stoichiometric natural gas (NG) engines were evaluated using synthetic exhaust flow reactor experiments. Under simulated NG exhaust conditions, NO conversion was intimately correlated to CH 4 conversion at slightly fuel-rich operating conditions. CH 4 conversion significantly varied with lambda (or O 2 concentration) and depended on the direction of the lambda change. The dynamic CH 4 conversion is likely related to the change of catalyst oxidation state and structures of Pd active sites. The CH 4 -NO cross-over point was found at a rich-biased lambda rather than stoichiometry. Compared to static operation, catalyst performance was much higher under a realistic dithering condition. The impacts of dithering parameters including amplitude and frequency on CH 4 and NO conversions were explored. O 2 dosage is proposed as an effective descriptor of the TWC dithering performance. In the O 2 -dosage space, the performance of distinct dithering parameters collapse, helping to elucidate the fundamental influence of dithering parameters on TWC performance. With the help of O2 dosage, an optimal operating window leading to high CH 4 and NO conversions was identified. A higher dithering amplitude is required to keep the catalyst in a reduced state for higher CH 4 and NO conversions, while an appropriate O 2 dosage that does not exceed the breakthrough OSC is necessary to avoid excess O 2 that results in NO slip. The dithering amplitude and O 2 -dosage-metric maps provide a new method for studying dynamic TWC performance and identifying optimum operation strategies.

03 NATURAL GAS↗

Impact of recent updates to neutrino oscillation parameters on the effective Majorana neutrino mass in 0 ν β β decay

We investigate how recent updates to neutrino oscillation parameters and the sum of neutrino masses influence the sensitivity of neutrinoless double-beta ( 0 ν β β ) decay experiments. Incorporating the latest cosmological constraints on the sum of neutrino masses and laboratory measurements on oscillations, we determine the sum of neutrino masses for both the normal hierarchy (NH) and the inverted hierarchy (IH). Our analysis reveals a narrow range for the sum of neutrino masses, approximately 0.06 eV / c 2 for NH and 0.102 eV / c 2 for IH. Utilizing these constraints, we calculate the effective Majorana masses for both NH and IH scenarios, establishing the corresponding allowed regions. Importantly, we find that the minimum neutrino mass is nonzero, as constrained by the current oscillation parameters. Additionally, we estimate the half-life of 0 ν β β decay using these effective Majorana masses for both NH and IH. Our results suggest that upcoming ton-scale experiments will comprehensively explore the IH scenario, while 100-ton-scale experiments will effectively probe the parameter space for the NH scenario, provided the background index can achieve 1 event/kton-year in the region of interest. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Cone beam neutron interferometry: From modeling to applications

Phase-grating moiré interferometers (PGMIs) have emerged as promising candidates for the next generation of neutron interferometry, enabling the use of a polychromatic beam and manifesting interference patterns that can be directly imaged by existing neutron cameras. However, the modeling of the various PGMI configurations is limited to cumbersome numerical calculations and backward propagation models which often do not enable one to explore the setup parameters. Here we generalize the Fresnel scaling theorem to introduce a k -space model for PGMI setups illuminated by a cone beam, thus enabling an intuitive forward propagation model for a wide range of parameters and experimental setups. The interference manifested by a PGMI is shown to be a special case of the Talbot effect, and the optimal fringe visibility is shown to occur at the moiré location of the Talbot distances. We derive analytical expressions for the contrast and the propagating intensity profiles in various conditions and provide the first analysis of the PGMI dark-field imaging signal when considering sample characterization. The model's predictions are compared to experimental measurements and good agreement is found between them. Last, we propose and experimentally verify a method to recover contrast at typically inaccessible PGMI autocorrelation lengths. The presented work provides a toolbox for analyzing and understanding existing PGMI setups and their future applications, for example extensions to two-dimensional PGMIs and characterization of samples with nontrivial structures. Published by the American Physical Society 2024

Sarenac, D. (ORCID:0000000185753367)↗

Search for Higgsinos in final states with low-momentum lepton-track pairs at 13 TeV

We present a search for the pair production of Higgsinos in final states with large missing transverse momentum and either two reconstructed muons or a reconstructed lepton (muon or electron) and an isolated track. The analyzed data correspond to proton-proton collisions with an integrated luminosity of 137 fb −1 , collected by the CMS experiment at $\sqrt{𝑠}$ =13 TeV in 2016, 2017, and 2018. The signal scenario assumes four nearly mass degenerate Higgsino mass eigenstates: two neutralino states $\tilde{𝜒}^0_2$ and $\tilde{𝜒}^0_1$ with a small mass difference in the range 1–10 GeV and two chargino states $\tilde{𝜒}^±_1$ with an intermediate mass. The analysis focuses on the decay of the heavier neutralino into the lighter one and a virtual 𝑍 boson, which decays into two same-flavor leptons. The leptons have small transverse momentum and/or a small opening angle between the identified muons. An isolated track is used to recover events in which only one of the two leptons is identified. Multivariate discriminants are used to enhance the sensitivity by efficiently rejecting backgrounds from SM processes or misreconstructed tracks and/or leptons. The search explores a unique phase space and probes a previously unexplored region of the signal model parameter space. Mass differences between the two neutralinos are probed down to 1.5 GeV, assuming a Higgsino mass of 100 GeV. The maximum excluded Higgsino mass is 115 GeV.

Hayrapetyan, A. [Yerevan Physics Institute]↗

The near-wake structure of the Atmosphere Explorer C /AE-C/ satellite - A parametric investigation

Measurements of ion current, electron temperature, and values of space potential obtained from the cylindrical electrostatic probe on board the Atmosphere Explorer C (AE-C) satellite were used to examine, in a parametric manner, the angular distribution of charge around the satellite. Interest is focused on nighttime equatorial data in the altitude range 275-620 km, which yields a wide range for the parameter R sub D (the radius of the satellite divided by the ambient value of the Debye lengths), including R sub D greater than 100, which is of practical significance to large space platforms. The variations of normalized ion current in the wake zone of the AE-C satellite appear to display an exponential dependence on R sub D for 'constant' values of other relevant parameters. The angular variations of electron temperature and space potential in the close vicinity of the satellite's surface were examined and compared with results from the Explorer 31 satellite. The variation of the ratio of measured to computed space potential with electron temperature was examined using data from both the AE-C and Explorer 31 satellites. It was found that the ratio is greater than unity. Possible causes for the above inequality are discussed.

Samir, U.↗

Currents and Flows in Distant Magnetospheres

Space scientists have explored, described, and explained the terrestrial magnetosphere for four decades. Rarely do they point out that the planetary and solar wind parameters controlling the size, shape, and activity of Earth's magnetosphere map out only a small portion of the space of dimensionless parameters that govern magnetospheric properties. With the discovery of Ganymede's magnetosphere, the range of parameters relevant to magnetospheric studies has grown by orders of magnitude. Consider the extremes of Ganymede's and Jupiter's magnetospheres. Jupiter's magnetosphere forms within a plasma flowing at super-Alfvenic speed, whereas Ganymede's forms in a sub-Alfvenic flow. The scale sizes of these magnetospheres, characterized by distances to the magnetopause of order 7x10(exp 6) km and 5x10(exp 3) km, respectively, differ by three orders of magnitude, ranging from 100 to 0.1 times the scale of Earth's magnetosphere. The current systems that control the structure and dynamics of a magnetosphere depend on specific plasma and field properties. Magnetopause currents at Ganymede differ greatly from the forms familiar for Earth and Jupiter, principally because the Mach number of the ambient plasma flow greatly influences the shape of the magnetosphere. A magnetodisk current, present at Jupiter because of its rapid rotation, is absent at Earth and Ganymede. The ring current, extensively investigated at Earth, is probably unimportant at Ganymede because the dynamical variations of the external flow are slow. The ring current is subsumed within the magnetodisk current at Jupiter. This paper describes and contrasts aspects of these and other current systems for the three bodies.

Kivelson, Margaret Galland↗

Unitarity bounds on effective field theories at the LHC

Effective Field Theory (EFT) extensions of the Standard Model are tools to compute observables (e.g. cross sections with partonic center-of-mass energy √s^) as a systematically improvable expansion suppressed by a new physics scale M. If one is interested in EFT predictions in the parameter space where M < √s^, concerns of self-consistency emerge, which can manifest as a violation of perturbative partial-wave unitarity. However, when we search for the effects of an EFT at a hadron collider with center-of-mass energy √s using an inclusive strategy, we typically do not have access to the event-by-event value of √s^. This motivates the need for a formalism that incorporates parton distribution functions into the perturbative partial-wave unitarity analysis. Developing such a framework and initiating an exploration of its implications is the goal of this work. Our approach opens up a potentially valid region of the EFT parameter space where M << √s. We provide evidence that there exist valid EFTs in this parameter space. The perturbative unitarity bounds are sensitive to the details of a given search, an effect we investigate by varying kinematic cuts.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Enhanced Multi-Dimensional Inversion Through Target-Specific Inversion Parameter Bounds With an Application to Crosswell Electromagnetic for Sequestration Monitoring

In geophysical inversions, lower and upper model parameter bounds are a means of solution stabilization. Further, constraints that intend to let only geologically plausible inverse solutions pass are amenable to lower and upper bounds. Reliable prior information is paramount to construct such bound constraints. It is common practice to narrow and widen bound intervals for regions of, respectively, more and less certain prior information. Contrary to this practice, we experiment with widened bound intervals in zones that are poorly resolved by a given survey configuration but where prior information would suggest structural anomalies of interest. The purpose of enlarged parameter bounds that correlate spatially with predefined targets is to let the inversion explore a larger solution space, thus increasing the potential to resolve otherwise hidden anomalies. Application of the method is based on a carbon-sequestration baseline (pre-injection) crosswell electromagnetic (EM) field survey at the Containment and Monitoring Institute Field Research Station (Alberta, Canada), where impeded measurements led to generally reduced sensitivities for the interwell region. Synthetic-data proofs of concept use augmented bounds designed to boost the resolution of artificial plume targets, indicating an enhanced illumination compared to constant bounds. Comparative field data inversions with spatially variable bounds constructed from prior resistivity and velocity information highlight non-horizontal baseline structures.

3-D electrical resistivity imaging↗

VarIabiLity seLection of AstrophysIcal sources iN PTF (VILLAIN): I. Structure function fits to 71 million objects

Light-curve variability is well-suited to characterizing objects in surveys with high cadence and a long baseline. This is especially relevant in view of the large datasets to be produced by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). We aim to determine variability parameters for objects in the Palomar Transient Factory (PTF) and explore differences between quasars (QSOs), stars, and galaxies. We relate variability and colour information in preparation for future surveys. We fit joint likelihoods to structure functions (SFs) of 71 million PTF light curves with a Markov chain Monte Carlo method. For each object, we assume a power-law SF and extract two parameters: the amplitude on timescales of one year, A, and a power-law index, γ. With these parameters and colours in the optical (Pan-STARRS1) and mid-infrared (WISE), we identify regions of parameter space dominated by different types of spectroscopically confirmed objects from SDSS. Candidate QSOs, stars, and galaxies are selected to show their parameter distributions. QSOs show high-amplitude variations in the R band, and the highest γ values. Galaxies have a broader range of amplitudes and their variability shows relatively little dependency on timescale. With variability and colours, we achieve a photometric selection purity of 99.3% for QSOs. Even though hard cuts in monochromatic variability alone are not as effective as seven-band magnitude cuts, variability is useful in characterizing object subclasses. Through variability, we also find QSOs that were erroneously classified as stars in the SDSS. We discuss perspectives and computational solutions in view of the upcoming LSST.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling Study of Hatch Spacing’s Effect on Grain Morphology in Repairing Damaged Metal Parts with Welding and Hatch Spacing’s Potential Use in Lunar Exploration

The welding process is a potential way of repairing a damaged metal component, especially cavity damage caused by a harsh environment like the lunar environment, which is characterized by large temperature differences and reduced gravity. The adjustment of welding parameter (e.g., hatch spacing) can improve production efficiency in the repair process. Seen from the microstructural level, hatch spacing sensitivity affects the metallic grain evolution and morphology in the welding process, which can further influence a repaired part’s mechanical properties; however, the study of hatch spacing’s effect on microstructure is challenging. Traditional experimental procedures are costly and time-consuming, and any change in hatch spacing value needs roll-back of experimental procedure. A modeling study can address the above challenges in experimental observation. In this research, a modeling approach based on the Kinetic Monte Carlo (KMC) Potts theory was used to simulate grain evolution and morphology with three hatch spacings. Through quantifying and analyzing the predicted grain morphologies, the effect of hatch spacing on microstructure in a welding-fabricated part was investigated. The predicted grain morphologies were validated with an EBSD image of welding microstructure, which has been published before. The primary grain morphologies were columnar grains with a small amount of fine equiaxed grains formed in the scanning path centerline. When increasing the hatch spacing, the columnar grains become larger and more lengthy, while the effect of hatch spacing on the equiaxed grains is not obvious.

Welding for repairing↗

Enhanced neutrino polarizability

We point out that neutrinos can have enhanced couplings to photons, if light (pseudo)scalar mediators are present, resulting in a potentially measurable neutrino polarizability. We show that the expected suppression from small neutrino masses can be compensated by the light mediator mass, generating dimension 7 Rayleigh operators at low scales. We explore the rich phenomenology of such models, computing in detail the constraints on the viable parameter space, spanned by the couplings of the mediator to neutrinos and photons. Finally, we build several explicit models that lead to an enhanced neutrino polarizability by modifying the inverse see-saw majoron, i.e., the pseudo-Nambu-Goldstone boson of the U(1) L global lepton number responsible for generating small neutrino masses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evaluation of Potential DAS Array Geometries for the Source Physics Experiment, Phase III

Distributed Acoustic Sensing (DAS) is a rapidly developing technology that can record acoustic wavefields at high sampling rates and with dense spatial spacing difficult to achieve with seismometers. However, the geophysical community has not fully explored DAS survey parameters to characterize the ideal array design. A better quantitative understanding of DAS array behavior prior to SPE Phase III acquisition can help improve the quality of the data recorded by guiding the DAS array design. Here we use array response functions as well as beamforming and backprojection results from forward modelling calculations to assess the performance of varying DAS array geometries to record regional and local sources. A seven-sided polygon DAS array demonstrated improved capabilities for recording regional sources over segmented linear arrays, with potential improvements in recording and locating local sources. These results help reveal DAS array performance as a function of geometry.

58 GEOSCIENCES↗

Performance of a parallel code for the Euler equations on hypercube computers

The performance of hypercubes were evaluated on a computational fluid dynamics problem and the parallel environment issues were considered that must be addressed, such as algorithm changes, implementation choices, programming effort, and programming environment. The evaluation focuses on a widely used fluid dynamics code, FLO52, which solves the two dimensional steady Euler equations describing flow around the airfoil. The code development experience is described, including interacting with the operating system, utilizing the message-passing communication system, and code modifications necessary to increase parallel efficiency. Results from two hypercube parallel computers (a 16-node iPSC/2, and a 512-node NCUBE/ten) are discussed and compared. In addition, a mathematical model of the execution time was developed as a function of several machine and algorithm parameters. This model accurately predicts the actual run times obtained and is used to explore the performance of the code in interesting but yet physically realizable regions of the parameter space. Based on this model, predictions about future hypercubes are made.

Barszcz, Eric↗

Monte Carlo modeling of the resurfacing of Venus

We have developed a three-dimensional model of venusian resurfacing that employs Monte Carlo simulations of both impact cratering and volcanism. The model simulates the production of craters on Venus by using the observed mass distributions of Earth- and Venus-crossing asteroids and comets. Lava flows are modeled by an energy minimization technique to simulate the effects of local topography on the shape and extent of flows. The model is run under a wide range of assumptions regarding the scale and time evolution of volcanism on Venus. Regions of the parameter space that result in impact crater distributions and modifications that are currently observed will be explored to place limits on the possible volcanic resurfacing history of Venus.

Bullock, M. A.↗