Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spherical Fitting”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Resistive heating in an electrified domain with a spherical inclusion: an ALEGRA verification study

A verification study is conducted for the ALEGRA software, using the problem of an electrified medium with a spherical inclusion, paying special attention to resistive heating. We do so by extending an existing analytic solution for this problem to include both conducting and insulating inclusions, and we examine the effects of mesh resolution and mesh topology, considering both body-fitted and rectangular meshes containing mixed cells. We present observed rates of convergence with respect to mesh refinement for four electromagnetic quantities: electric potential, electric field, current density and Joule power.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Discovery of multiple p-mode pulsation frequencies in the roAp star, HD 86181

ABSTRACT We report the frequency analysis of a known roAp star, HD 86181 (TIC 469246567), with new inferences from Transiting Exoplanet Survey Satellite (TESS) data. We derive the rotation frequency to be νrot = 0.48753 ± 0.00001 d−1. The pulsation frequency spectrum is rich, consisting of two doublets and one quintuplet, which we interpret to be oblique pulsation multiplets from consecutive, high-overtone dipole, quadrupole, and dipole modes. The central frequency of the quintuplet is 232.7701 d−1 (2.694 mHz). The phases of the sidelobes, the pulsation phase modulation, and a spherical harmonic decomposition all show that the quadrupole mode is distorted. Following the oblique pulsator model, we calculate the rotation inclination, i, and magnetic obliquity, β, of this star, which provide detailed information about the pulsation geometry. The i and β derived from the best fit of the pulsation amplitude and phase modulation to a theoretical model, including the magnetic field effect, slightly differ from those calculated for a pure quadrupole, indicating the contributions from ℓ = 4, 6, 8,... are small. Non-adiabatic models with different envelope convection conditions and physics configurations were considered for this star. It is shown that models with envelope convection almost fully suppressed can explain the excitation at the observed pulsation frequencies.

Shi, Fangfei (ORCID:0000000323998216)↗
DOE OSTI2025

EMPDF : inferring the Milky Way mass with data-driven distribution function in phase space

We introduce the emPDF (empirical distribution function), a novel dynamical modelling method that infers the gravitational potential from kinematic tracers with optimal statistical efficiency under the minimal assumption of steady state. emPDF determines the best-fitting potential by maximizing the similarity between instantaneous kinematics and the time-averaged phase-space distribution function (DF), which is empirically constructed from observation upon the theoretical foundation of oPDF (Han et al. 2016). This approach eliminates the need for presumed functional forms of DFs or orbit libraries required by conventional DF- or orbit-based methods. emPDF stands out for its flexibility, efficiency, and capability in handling observational effects, making it preferable to the popular Jeans equation or other minimal assumption methods, especially for the Milky Way (MW) outer halo where tracers often have limited sample size and poor data quality. We apply emPDF to infer the MW mass profile using Gaia DR3 data of satellite galaxies and globular clusters, obtaining enclosed masses of M (,r) = 26±8, 46±8, 90±13⁠, and 149±40 x 10 10 M ⊙ at r = 30, 50, 100⁠, and 200 kpc, respectively. These are consistent with the updated constraints from simulation-informed DF fitting (Li et al. 2020). While the simulation-informed DF offers superior precision owing to the additional information extracted from simulations, emPDF is independent of such supplementary knowledge and applicable to general tracer populations. emPDF is currently implemented for tracers with complete 6D kinematics within spherical potentials, but it can potentially be extended to address more general problems.

Astrophysics of Galaxies (astro-ph.GA)↗
DOE OSTI2023

CHIPFoam Modification to Allow for Modelling SWIFT Foams

SWIFT (Silicon-water in Familiar Template) foams are silicone foams that have a continuous porosity without having any obvious spherical or ellipsoidal pore shapes. Like other foamed rubber, SWIFT foams can experience large reversible deformations. Because of their microstructure, SWIFT foams have a stiffering behavior in compression that is very abrupt compared to most other foamed rubbers. Work by Benedikt reported by Miller showed that fitting compressive test data with the CHIPFoam model in its current form is impractical and inaccurate. This report describes an extension made to CHIPFoam to allow the modeling of SWIFT foams. While the extension allows the modeling of SWIFT foams, it also extends the range of porosity of foams that can be modelled with CHIPFoam to include very high porosity (>99%) foamed rubber. The previous model was unstable for porosities greater than 70%. The changes made to CHIPFoam include the introduction of a critical porosity at which the compressible Danielsson strain energy function that couples volumetric and deviatoric mechanical behavior begins to increase with compression. This new feature also required the introduction of a third term in the multiplicative decomposition of the relative volume, J . The organization of this report is described here. The first section, SWIFT Foam Structure and Mechanical Behavior, describes the structure of SWIFT foams as contrasted to other foamed rubbers. It also contrasts the mechanical behavior of SWIFT foams to other more traditional foamed rubbers. The second section, CHIPFoam Modifications, describes the specific modifications done to be able to modify the delayed rapid compressive stiffening seen in SWIFT foams. The third section, Results, demonstrates the use of the modified model to fit compressive loading response data from a SWIFT foam sample. The fourth section, Conclusions and Further Work, puts the work in perspective and presents likely further work.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Adsorption of Mixed Micelles of Polysorbate 80 and Oleic Acid to the Air–Water Interface

The solubilization of long chain amphiphiles with limited solubility into mixed micelles composed of highly soluble surfactants plays a crucial role in modulating the stability and functionality of formulations in pharmaceutical and food systems. Subsequently, the mixed micelles adsorb from solution to the air−water surface, defining the transport mechanism of the insoluble amphiphiles to the interface. We use X-ray reflectivity to measure the composition and provide insight into the structure of these mixed monolayers, and, alongside interfacial tension measurements, we provide an understanding of how the bulk composition determines the surface composition and the dynamics of tension reduction. We use a model system consisting of an insoluble fatty acid, oleic acid (OA), and a soluble micelle-forming surfactant, polysorbate 80 (PS80). PS80 forms spherical micelles, and, above the critical micelle concentration (cmc), OA is readily solubilized inside the micelles. We show that the adsorption of PS80/OA mixed micelles to the air−water interface rapidly reduces the tension and lowers the equilibrium tension in proportion to the OA concentration. X-ray reflectivity data, fit using Parratt-slab models, quantitively demonstrates that the monolayers become enriched with OA, and we show how the OA intercalates into the PS80 monolayers.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Uncertainty-Aware Machine Learning for Small-Angle X-ray Scattering Analysis in Autonomous Experimentation

Small-angle X-ray scattering (SAXS) is a powerful high-throughput characterization tool for probing nanoscale structure in native sample environments, providing real-time morphological information such as nanoparticle size and shape during synthesis. However, automated SAXS data analysis for extracting meaningful structural parameters is non-trivial and remains a bottleneck in closed-loop experimentation towards autonomous materials discovery, which demands fast, reliable, and uncertainty-aware data analysis. Here, we develop a machine-learning approach for automated SAXS analysis tailored to closed-loop nanoparticle synthesis. A Random Forest (RF) regression model is trained on 100,000 synthetic SAXS curves generated from polydisperse spherical nanoparticles with realistic background contributions. Using normalized one-dimensional SAXS intensity profiles as input, the RF model directly predicts nanoparticle radius, size polydispersity, and background parameters, while the ensemble standard deviation across trees provides built-in uncertainty quantification (UQ). On synthetic data, we show that combining fit-quality metrics (R 2 , MAE) with thresholds on prediction uncertainty reliably identifies accurate parameter estimates without access to ground truth. We then apply the trained model to 365 experimental SAXS profiles of citrate-reduced gold nanoparticles synthesized using an automated droplet-flow microreactor with in situ SAXS at a synchrotron beamline, classifying the results into high- and low-confidence subsets based on UQ metrics. Finally, we integrate RF-based SAXS analysis into a simulated closed-loop optimization campaign using Gaussian process Bayesian optimization to minimize nanoparticle polydispersity, benchmarking against conventional automated Levenberg–Marquardt fitting. The RF-guided campaign exhibits substantially faster convergence and lower relative opportunity cost (∼0.07 vs ∼0.3), demonstrating that uncertainty-aware machine-learning SAXS analysis significantly enhances the efficiency and robustness of autonomous nanomaterials synthesis workflows.

Bayesian optimization↗
DOE OSTI2023

Effects of head modeling errors on the spatial frequency representation of MEG

Abstract Objectives. We aim to investigate the effects of head model inaccuracies on signal and source reconstruction accuracies for various sensor array distances to the head. This allows for the assessment of the importance of head modeling for next-generation magnetoencephalography (MEG) sensors, optically-pumped magnetometers (OPM). Approach. A 1-shell boundary element method (BEM) spherical head model with 642 vertices of radius 9 cm and conductivity of 0.33 S m −1 was defined. The vertices were then randomly perturbed radially up to 2%, 4%, 6%, 8% and 10% of the radius. For each head perturbation case, the forward signal was calculated for dipolar sources located at 2 cm, 4 cm, 6 cm and 8 cm from the origin (center of the sphere), and for a 324 sensor array located at 10 cm to 15 cm from the origin. Equivalent current dipole (ECD) source localization was performed for each of these forward signals. The signal for each perturbed spherical head case was then analyzed in the spatial frequency domain, and the signal and ECD errors were quantified relative to the unperturbed case. Main results. In the noiseless and high signal-to-noise ratio (SNR) case of approximately ≥6 dB, inaccuracies in our spherical BEM head conductor models lead to increased signal and ECD inaccuracies when sensor arrays are placed closer to the head. This is true especially in the case of deep and superficial sources. In the noisy case however, the higher SNR for closer sensor arrays allows for an improved ECD fit and outweighs the effects of head geometry inaccuracies. Significance. OPMs may be placed directly on the head, as opposed to the more commonly used superconducting quantum interference device sensors which must be placed a few centimeters away from the head. OPMs thus allow for signals of higher spatial resolution to be captured, resulting in potentially more accurate source localizations. Our results suggest that an increased emphasis on accurate head modeling for OPMs may be necessary to fully realize its improved source localization potential.

Yeo, Wan-Jin (ORCID:0000000241417123)↗
DOE OSTI2026

The impact of kinetic and global effects on ideal ballooning 2nd stable pedestals of conventional and low aspect-ratio tokamaks

The EPED model (Snyder et al 2011 Nucl. Fusion 51 103016) has had success in describing the pedestal structure of Type-I ELM and QH-mode plasmas in tokamaks, by combining kinetic ballooning mode (KBM) and peeling–ballooning (PB) constraints. Within EPED, the KBM constraint is often approximated by a technique using calculated ideal ballooning mode (IBM) thresholds, fit to a functional form designed to capture non-local effects. It has been noted that quantitative differences between local ideal MHD and gyro-kinetic (GK) ballooning stability can be larger at low aspect ratio. KBM critical pedestals are consistent with observations in initial studies on conventional and spherical tokamaks. In this work, the application of a reduced model for the calculation of the kinetic ballooning stability boundary is presented based on a novel and newly developed gyro-fluid system (GFS) code (Staebler et al 2023 Phys. Plasmas 30 102501). GFS is observed to capture KBMs in DIII-D as well as the NSTX pedestals, opening a route to integrating this model into EPED. Finally, high but finite $n$ global ballooning modes are observed to limit the access to the local 2nd stability and thus provide a transport mechanism that constrains the pedestal evolution with $β_{p,\textrm{ped}}$. The high $n$ global ballooning stability is approximated by its ideal MHD analog using ELITE. It is shown that nearly-local high $n$ modes with $k_yρ_s$ ~ $0.25 - 0.5$ can provide a proxy for the critical $β_{p,\textrm{ped}}$ when a 2nd stable access exists on DIII-D plasmas. The use of GFS and ELITE scaling in EPED provides improved agreement to EPED1 in an initial comparison with DIII-D pedestal data.

Anastopoulos Tzanis, Michail Savvas [Oak Ridge Nat↗
DOE OSTI2024

A magnetic diagnostic suite for the Pegasus-III experiment

Pegasus-III is an ultralow aspect ratio spherical tokamak providing a dedicated US experiment for comparative solenoid-free startup studies. A new magnetic diagnostic suite for equilibrium and low frequency (<200 kHz) magnetohydrodynamic mode analysis has been installed. These new diagnostics address the significant challenges of measuring magnetic field in a high noise environment with the majority constrained to fit in an 8 mm diagnostic gap on the high field side. Electrostatic switching noise generated by the 16 independent current feedback-controlled power supplies produces dV cm /dt ~ 1 kV/μs and volt level common mode noise on the magnetics. Immunity to this switching noise is accomplished through differential signal runs and signal processing, along with end-to-end electromagnetic interference shielding. The magnetic measurements are simultaneously digitized at 1 MHz and conditioned by precision 8 pole Butterworth filters with a corner frequency of 200 kHz to prevent aliasing down to the 16-bit level over the full passband. Ex-vessel calibrations of the B p coils were completed with a typical uncertainty of <0.5%. Stray toroidal field pickup from coil misalignment or positioning errors is corrected using a physics-based model. Comparisons of the corrected measurements to modeling agree to within 1.3% on average. This is within the 1.5% measurement uncertainty that a sensitivity analysis determined is needed for accurate fast boundary and equilibrium reconstruction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Synthetic modeling of soft x-ray emissivity for magnetic island analysis in LTX-β

We present a synthetic soft x-ray (SXR) forward-modeling framework that characterizes the emissivity structure of rotating magnetic islands in the Lithium Tokamak eXperiment-β across space and time. Magnetic islands associated with tearing modes produce modulations in line-integrated SXR brightness. Traditional tomographic methods struggle to resolve these structures in devices with limited sightlines, resulting in an under-determined inversion problem, or otherwise necessitate equilibrium reconstruction data, preventing their use in active control. Our approach removes this challenge by combining a fully three-dimensional ray-tracing model with a time-dependent emissivity prescription of a $m/n$ = 2/1 magnetic island. The model is validated through helical island geometry, rotation measured by magnetic diagnostics, and equilibrium constraints from PSI-Tri reconstructions. Synthetic brightness signals are generated for each photodiode sightline from a modeled emissivity profile and directly compared with experimental data from a tangential midplane SXR array. By fitting the synthetic diagnostic output to the observed brightness across time, we simultaneously infer key geometric parameters-such as island radial location, island width, and rotation frequency-without requiring full tomographic reconstruction, all with an average absolute deviation of less than 3%. This work demonstrates that forward modeling with a single tangential array can extract key island parameters in a spherical tokamak, providing a computationally efficient alternative to conventional SXR tomography and providing a pathway toward real-time magnetic-island characterization in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analytic Gaussian covariance matrices for galaxy N-point correlation functions

Here, we derive analytic covariance matrices for the N-point correlation functions (NPCFs) of galaxies in the Gaussian limit. Our results are given for arbitrary N and projected onto the isotropic basis functions given by spherical harmonics and Wigner 3j symbols. A numerical implementation of the 4PCF covariance is compared to the sample covariance obtained from a set of lognormal simulations, Quijote dark matter halo catalogues, and MultiDark-Patchy galaxy mocks, with the latter including realistic survey geometry. The analytic formalism gives reasonable predictions for the covariances estimated from mock simulations with a periodic-box geometry. Furthermore, fitting for an effective volume and number density by maximizing a likelihood based on Kullback-Leibler divergence is shown to partially compensate for the effects of a nonuniform window function. Our result is recently shown to facilitate NPCF analysis on a realistic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bottlebrush Block Copolymers at the Interface of Immiscible Liquids: Adsorption and Lateral Packing

Amphiphilic bottlebrush block copolymers (BBCPs), having a hydrophilic bottlebrush polymer (BP) linked covalently to a hydrophobic BP, were found to segregate to liquid-liquid interfaces to minimize the free energy of the system. The key parameter influencing the outcome of the experiments is the ratio between the degree of polymerization of the backbone (NBB) and that of the side-chain brushes (NSC). Specifically, a spherical, star-like configuration results when NBB < NSC, while a cylindrical, bottlebrush-like shape is preferred when NBB > NSC. Dynamic interfacial tension (γ) and fluorescence recovery after photobleaching (FRAP) measurements show that the BBCP configuration influences the areal density and in-plane diffusion at the fluid interface. The characteristic relaxation times associated with BBCP adsorption (τA) and reorganization (τR) were determined by fitting time-dependent interfacial tension measurements to a sum of two exponential relaxation functions. Both τA and τR initially increased with NBB up to 92 repeat units, due to the larger hydrodynamic radius in solution and slower in-plane diffusivity, attributed to a shorter cross-sectional diameter of the side-chains near the block junction. This trend reversed at NBB = 190, with shorter τA and τR attributed to increased segregation strength and exposure of the bare water/toluene interface due to tilting and/or wiggling of the backbone chains, respectively. The adsorption energy barrier decreased with higher NBB, due to a reduced BBCP packing density at the fluid interface. In conclusion, this study provides fundamental insights into macromolecular assembly at fluid interfaces, as it pertains to unique bottlebrush block architectures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Quenching factor measurements of neon nuclei in neon gas

The NEWS-G collaboration uses spherical proportional counters (SPCs) to search for weakly interacting massive particles (WIMPs). Here, in this paper, we report the first measurements of the nuclear quenching factor in neon gas at 2 bar using an SPC deployed in a neutron beam at the TUNL facility. The energy-dependence of the nuclear quenching factor is modeled using a simple power law: $αE^{β}_{nr}$; we determine its parameters by simultaneously fitting the data collected with the detector over a range of energies. We measured the following parameters in Ne : CH 4 at 2 bar: α= 0.2801±0.0050 (fit) ±0.0045 (sys) and β=0.0867±0.020 (fit) ±0.006 (sys). Our measurements do not agree with expected values from SRIM or Lindhard theory. We demonstrated the feasibility of performing quenching factor measurements at sub-keV energies in gases using SPCs and a neutron beam.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Development of an x-ray radiography platform to study laser-direct-drive energy coupling at the National Ignition Facility

A platform has been developed to study laser-direct-drive energy coupling at the National Ignition Facility (NIF) using a plastic sphere target irradiated in a polar-direct-drive geometry to launch a spherically converging shock wave. To diagnose this system evolution, eight NIF laser beams are directed onto a curved Cu foil to generate Heα line emission at a photon energy of 8.4 keV. These x rays are collected by a 100-ps gated x-ray imager in the opposing port to produce temporally gated radiographs. The platform is capable of acquiring images during and after the laser drive launches the shock wave. A backlighter profile is fit to the radiographs, and the resulting transmission images are Abel inverted to infer radial density profiles of the shock front and to track its temporal evolution. Here, the measurements provide experimental shock trajectories and radial density profiles that are compared to 2D radiation-hydrodynamic simulations using cross-beam energy transfer and nonlocal heat-transport models.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Effects of Injected Current Streams on MHD Equilibrium Reconstruction of Local Helicity Injection Plasmas in a Spherical Tokamak

Open field line currents are intrinsic to DC helicity injection plasma startup and pose a challenge for inferring the plasma equilibrium with standard reconstruction analysis. Local helicity injection (LHI) is a type of DC helicity injection which uses small, modular current sources to drive force-free current along helical field lines to produce tokamak plasmas. MHD modeling and magnetic measurements during LHI indicate the injected current streams remain coherent as helical structures on the outboard edge of a core toroidal plasma that is tokamak-like in a toroidally averaged sense. To extract core plasma equilibrium properties, external magnetic diagnostics corrected for contributions from the injected current streams are fitted by a standard Grad-Shafranov equilibrium code. An iterative approach for estimating and subtracting the stream contributions from the diagnostic signals is described and applied to a model equilibrium database to reduce systematic errors introduced by the streams. Convergence is usually attained with 2 to 4 iterations, with derived equilibrium parameters matching the prescribed axisymmetric core values to within estimated experimental uncertainties. Accurate recovery of core parameters occurs when the ratio of the net toroidal windup current from the streams to the core plasma current is less than 0.2, which is typically satisfied in most experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Main Bearing Replacement and Damage - A Field Data Study on 15 Gigawatts of Wind Energy Capacity

This study seeks to establish a comprehensive baseline of knowledge for the replacement and damage of main bearings in wind turbines. The purpose of this report is to provide a high-level summary of the data set, methodology, and results of this work. Full technical details and an extended analysis will be made available in a future publication. We collected data on main bearing replacements and reported damage from industrial partners based in Europe and the United States. In total, we obtained data for 167 wind power plants, with a combined capacity of 15.3 gigawatts (GW). Most of the data set was comprised of land-based, three-point mount, spherical roller bearings. Within this data set were 689 instances of main bearing replacement. Analysis was undertaken in two parts: first, a statistical analysis of the main bearing time-to-replacement data using survival analysis techniques; second, quantitative and qualitative analyses of the obtained damage information. Our results showed that 10% of a fixed main bearing population would be expected to have been replaced by 10.5 years. This is close to half of the 20-year design value. Fitted parametric distributions then indicated that by year 20, between 22% and 25% of main bearings are expected to have been replaced. Analysis of the damage reports revealed spalling to be the main type of damage listed. The additional presence of surface damage in the collected data indicates that at least part of the spalling cases are likely due to surface-initiated rolling contact fatigue. At this stage is not clear what proportion of spalling cases result from "wear induced", surface-initiated and subsurface-initiated rolling contact fatigue. While this work provides important insights into the current state of main bearing replacements and damage, many questions remain. An ongoing and expanding data collection and analysis effort focused on wind turbine main bearings is therefore recommended.

17 WIND ENERGY↗
DOE OSTI2023

Propellant Mass Gauging in a Spherical Tank under Micro-Gravity Conditions Using Capacitance Plate Arrays and Machine Learning

Propellant mass gauging under micro-gravity conditions is a challenging task due to the unpredictable position and shape of the fuel body inside the tank. Micro-gravity conditions are common for orbiting satellites and rockets that operate on limited fuel supplies. Capacitance sensors have been investigated for this task in recent years; however, the effect of various positions and shapes of the fuel body is not analyzed in detail. In this paper, we investigate this with various fill types, such as annular, core-annular, and stratified fills at different positions. We compare the performance among several curve-fitting-based approaches and a machine-learning-based approach, the latter of which offers superior performance in estimating the fuel content.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Overview of the Langmuir probe system on the Mega Ampere Spherical Tokamak (MAST) Upgrade

A detailed description of the Langmuir probe system on Mega Ampere Spherical Tokamak Upgrade is presented. The system features 850 tile-embedded probes and 40 bespoke electronic modules that each have the capability to drive and acquire data from up to 16 probes in a time-multiplexed manner. The system provides spatiotemporal-resolved measurements (1 cm and ~1 ms, respectively) in the divertor region of ion saturation current, electron temperature, and floating potential. The standard interpretation of current-voltage (IV) characteristics is to apply a four-parameter fit, based on unmagnetized probe theory, which includes a linear model for the ion saturation region. In conclusion, to mitigate the effect of the magnetic field, analysis is restricted to the region of the IV characteristic, which is sensitive to only the tail of the electron energy distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗