Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Profiling”

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 163 records · Page 9

Development of a Time-Height Profile Assimilation Technique for Large-Eddy Simulation

Mesoscale-to-microscale coupling (MMC) aims to address the limited scope of traditional large-eddy simulations by driving the microscale flow with information concerning large-scale weather patterns provided by mesoscale models. We present a new offline MMC technique for horizontally homogeneous microscale flow conditions, in which internal forcing terms are computed based on mesoscale time–height profiles of mean-flow quantities. The advantage of such an approach is that it can be used to drive a microscale simulation with either mesoscale or observational data, and that it does not rely on specific terms in the mesoscale budget equations, which are typically not part of the default output of a mesoscale solver. The performance of the proposed profile assimilation technique is assessed based on the simulation of a typical diurnal cycle over the Scaled Wind Farm Technology site in west Texas. Results indicate that simple data assimilation techniques lead to unphysically high levels of shear and turbulence caused by the algorithm’s inability to cope with inaccuracies in the mesoscale time–height profiles. Modifying the algorithm to account for vertical coherence in the mesoscale source terms gives the microscale solver a greater ability to correct the provided mesoscale time–height profiles, leading to improved predictions of shear and turbulence statistics. The resulting turbulence statistics are in good agreement with meteorological tower observations and simulation results obtained with state-of-the-art coupling techniques using mesoscale budget components.

data assimilation↗

Sample environment effects on synchrotron-measured temperature profiles in an approximant of optical floating zone crystal growth

Even though the growth of crystals using optical floating zone furnaces has had an immense scientific impact, the implementation of this method remains more of an art than a science due to the difficulty of obtaining quantitative information about the sample thermal profile during crystal growth. Building on recent work demonstrating that in situ synchrotron studies can be used to map sample rod temperatures during heating, investigations were carried out to better understand how the sample environment affects the sample temperature profile. Through a combination of experimental studies and modeling efforts, it is shown that the environment in the furnace can strongly influence the sample temperature at the lamp focus, the steepness of the vertical temperature gradient, and the timescale required for sample heating and cooling - effects which can combine to produce a strong history-dependence to sample temperature profiles. It is demonstrated that the furnace effects can be effectively captured in thermal models, allowing both the steady-state and time-dependent behavior of the sample to be accurately reproduced with predictive models, and providing a launching point for improved furnace designs that can more readily deliver desired thermal profiles.

36 MATERIALS SCIENCE↗

Tutorial: Machine-Learning-Based CREASE-2D Analysis of 2D SAXS Profiles to Characterize Anisotropic Nanostructures in Soft Materials

We present a tutorial to guide users on how to extend the Computational Reverse Engineering Analysis of Scattering Experiments-2D (CREASE-2D) framework to interpret their experimental two-dimensional small-angle scattering (SAS) data from soft materials (e.g., polymers, peptide amphiphiles, biomolecular fibrils). Unlike most traditional SAS analysis approaches, which typically rely on azimuthally averaged onedimensional (1D) profiles, CREASE-2D utilizes the complete 2D scattering profile to reveal information about anisotropy in the structure. In past applications, CREASE has provided insights into complex structural features, including the cross-sectional shapes of assembled nanostructures and dispersity in these features, which are difficult to discern with existing analytical models. While (1D- ) CREASE has been applied to SANS and SAXS data, this tutorial shares the steps for implementing CREASE-2D using an example of a dipeptide solution system, for which we have SAXS data. We present details for these steps involved in using CREASE-2D to interpret SAXS profiles: how to preprocess SAXS data, define relevant structural features, generate three-dimensional real-space structures for specific values of these features, train a machine learning (ML) surrogate model to predict scattering profiles for given structural features, and optimize these features using genetic algorithms (GA). Then, we use these steps to interpret complex 2DSAXS data collected from dipeptide solutions that, in microscopy images, exhibit nanoscale structures that could be elliptical tubes/ flat tapes/cylinders or a combination of these cross sections. Open-source codes, computational hardware, and software requirements, as well as the strengths and limitations of this protocol, are also presented. We expect researchers working with (soft) biomaterials, peptide amphiphiles, amphiphilic polymer solutions, polymer nanocomposites, and blends of particles/polymers will find this CREASE-2D method and this tutorial of use.

CREASE↗

MHD stability of spherical tokamak equilibria with non-monotonic q-profiles

We use the 3D magnetohydrodynamic (MHD) code M3D-C1 [Jardin et al., Comput. Sci. Discovery 5, 014002 (2012)] to examine the MHD stability and subsequent evolution of NSTX shot 129169. This discharge had a period with a non-monotonic safety factor profile, q (reversed shear), which was terminated by a MHD event that abruptly lowered the central safety factor, q0, and greatly reduced the peakedness of the pressure profile. We show that the equilibrium just before the MHD event occurred was linearly unstable to many pressure-driven infernal modes. Modes with toroidal mode number n≥3 all had rational surfaces very close to the minimum value of q. However, a non-resonant pressure-driven (1, 1) mode was also present, and this dominated the nonlinear evolution. The final state in the simulation, after the MHD activity subsided, had a reduced and flattened pressure profile and a nearly monotonic q-profile, in qualitative agreement with experimental results. The initial state was also unstable to the resistive interchange criteria in the reversed-shear region, but the final state was stable everywhere. The “double tearing mode” (DTM) does not appear to play a role in the MHD activity of this discharge. In Appendix A, we show that in a torus, the DTM is strongly stabilized by pressure, but it is destabilized in cylindrical geometry (which has been the most extensively analyzed in the literature).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inelastic deformation of diamond single crystals shock compressed to multimegabar stresses: Wave profile calculations

As the archetypal strong solid, the response of diamond shock compressed to multimegabar stresses is important for fundamental science and for numerical simulations of wave profiles for applications in high energy density physics experiments. Previous experiments and analysis have shown that the commonly used hydrodynamic assumption is invalid for diamond shock compressed to stresses below melt and an elastic–inelastic description is needed. Here, we present a phenomenological material model for calculating wave profiles in shock compressed diamond single crystals that incorporates this description. Also, to support the modeling effort, we carried out wave profile measurements on shock compressed diamond single crystals at the Sandia Z facility to augment previous measurements. Wave profiles for [100] and [111] diamond calculated using the material model provide a good match to the elastic–inelastic response (observed two-wave structure) measured at ∼325 and ∼360 GPa. Furthermore, the calculated peak stresses for single (overdriven) waves provide a good match to the measured Hugoniot states for stresses reaching ∼700 GPa, which is near melting conditions. The present results show that the diamond single crystal response at multimegabar shock stresses is characteristic of a brittle solid—pressure-dependent strength and strength loss due to inelastic deformation.

Deformation↗

Measurements of radial neutral density profiles from Balmer- α emission in Wendelstein 7-X

Radial neutral density profiles are estimated from measurements of passive H α emission in the Wendelstein 7-X stellarator. To parametrize the generally three-dimensional distribution with a low number of degrees of freedom, the neutral density is reduced to a flux surface quantity. Accounting for emission from excitation and recombination processes, neutral density profiles are derived independently for each of the available lines of sight. Density profiles obtained from the different viewing geometries are found to vary within one order of magnitude. Toroidally oriented lines of sight predict systematically lower neutral densities when compared to poloidally oriented ones. This discrepancy is attributed to the simplifications inherent in the imposed model and significant differences in integration volumes across the viewing geometries. In line with expectations, obtained neutral densities are found to decrease with increasing plasma density. Key restrictions of the model include the reduction of the neutral density to a flux surface quantity, uncertainties in the plasma profiles and instrument function, and line integration effects outside the last closed flux surface.

Wendelstein 7-X↗

Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature (T e ) and ion temperature (T i ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts T e by 46%, with much larger variance than MMM, and underpredicts T i by 25%. As the ratio of kinetic to magnetic field pressure (β) is increased across NSTX discharges, TGLF predicts lower T e and significant flattening of the T i profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both T e and T i are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Furthermore, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.

Lestz, J. B. [General Atomics, San Diego, CA (Unit↗

Modeling electron temperature profiles in the pedestal with simple formulas for ETG transport

Abstract This paper reports on the refinement (building on (Hatch D.R. et al 2022 Phys. Plasmas 29 062501)) and application of simple formulas for electron heat transport from electron temperature gradient (ETG) driven turbulence in the pedestal. The formulas are improved by (1) improving the parameterization for certain key parameters and (2) carefully accounting for the impact of geometry and shaping in the underlying gyrokinetic simulation database. Comparisons with nonlinear gyrokinetic simulations of ETG transport in the MAST pedestal demonstrate the model’s applicability to spherical tokamaks in addition to standard aspect ratio tokamaks. We identify bounds for model applicability: the model is accurate in the steep gradient region, where the ETG turbulence is largely slab-like, but accuracy decreases as the temperature gradient becomes weaker in the pedestal top. We use the formula to model the electron temperature profile in the pedestal for four experimental scenarios while extensively varying input parameters to represent uncertainties. In all cases, the predicted electron temperature profile exhibits extreme sensitivity to separatrix temperature and density, which has implications for core-edge integration. The model reproduces the electron temperature profile for high η e = L n e / L T e scenarios but not for low η e scenarios in which microtearing modes (MTMs) have been identified. We develop a proof-of-concept model for MTM transport and explore the relative roles of ETG and MTM in setting the electron temperature profile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Gyrokinetic profile prediction and validation of a negative triangularity plasma in ASDEX Upgrade

In this work, gyrokinetic simulations are performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73–93) for a negative triangularity H-mode plasma in ASDEX Upgrade, and compared with experimental measurements. The PORTALS framework (Rodriguez-Fernandez et al 2024 Nucl. Fusion 64 076034) is used to accelerate the prediction of kinetic profiles for this plasma, using surrogate modeling and Bayesian optimization. Ion heat flux, electron heat flux, and electron particle flux are simultaneously matched across the simulated radial regime of the plasma (normalized radius $r/a = 0.35-0.90$), and the resulting ion temperature, electron temperature, and electron density profiles match well with the experimental profile data within this radial range. A synthetic Correlation Electron Cyclotron Emission diagnostic is applied to find well-matched electron temperature fluctuation properties between simulation and experiment. The flux-matched profiles provide a basis for investigation of the turbulence nature across the plasma radius, revealing the dominance of Trapped Electron Mode turbulence at $r/a = 0.35$, the dominance of Ion Temperature Gradient turbulence at $r/a = 0.55$, 0.75, and 0.83, and an instability boundary at $r/a = 0.90$.

gyrokinetic simulation↗

Nonlinear Alfvén instability simulation and EP transport for ITER reversed shear (steady-state) and monotonic q-profile regimes *

The nonlinear evolution and dynamics of energetic particle (EP) driven instabilities in both steady-state (reversed shear q-profile) and monotonic q-profile regimes of ITER are examined using the global gyro-Landau closure model FAR3d. Both neutral beam and alpha components are included, leading to synergistic effects between the two populations. In order to maintain computational feasibility, the present nonlinear simulation includes all toroidal mode numbers ranging from n = 0 to 15. It cannot be excluded that n’s above this range could introduce larger linear growth rates; however, as indicated from the results, the n = 0 to 15 range is sufficient to provide a strong linear drive and significant nonlinear transport effects. While the toroidal mode numbers n = 6,13 dominate the linear instability growth phase, nonlinear energy transfers in the saturated phase reverse this trend, leading to the dominance of lower n’s in the saturated phase of the simulation. Zonal flow structure generation and EP density flattening are further consequences of the nonlinear phase. The density profile flattening is caused by collectively driven EP radial transport fluxes, which have non-local characteristics. Instantaneous alpha particle transport fluxes are significant (Γ α ~ 2 x 10 20 m –2 s –1 ) for the reversed shear regime and lower (Γ α ~ 1.1 x 10 20 m –2 s –1 ) for the monotonic q-profile case.

AE↗

Density profile of dynamical halos

Among the most fundamental properties of a dark matter halo is its density profile. Motivated by the recent proposal by García et al . [Mon. Not. R. Astron. Soc. 521 , 2464 (2023)] to define a dynamical halo as the collection of orbiting particles in a gravitationally bound structure, we characterize the mean and scatter of the orbiting profile of dynamical halos as a function of their orbiting mass. Here, we demonstrate that the orbiting profile of individual halos at fixed mass depends on a single dynamical variable, the halo radius 𝑟 h , which characterizes the spatial extent of the profile. The scatter in halo radius at fixed orbiting mass is ≈16%. Only a small fraction of this scatter arises due to differences in halo formation time, with late-forming halos being more compact (smaller halo radii). Accounting for this additional correlation results in an ≈11% scatter in halo radius at fixed mass and halo formation time.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimation of Forest Aboveground Biomass from Derivatives of Vegetation-Structure Profiles

Several studies have found that the vertical Fourier transform of lidar, interferometric Synthetic Aperture Radar (SAR), and stereo photogrammetric profiles at empirically-determined spatial frequencies enables high-performance forest aboveground biomass (AGB) estimation. Linear combinations of real and imaginary parts of Fourier transforms of Tomographic (multi-baseline) SAR (TomoSAR) profiles, from Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) airborne data, generate ~20%-precision estimates of AGB in the Saskatchewan area of Canada. We found that this 20% precision can be improved to ~15%, a factor of 30% improvement in root mean square error (RMSE) if, in addition to using Fourier transforms of the profile itself, we use Fourier transforms of the spatial, vertical derivative of the profile. The formulation of this "derivative" algorithm is the subject of this paper.

Treuhaft, Robert↗

Weather Sensitive High Spatio-Temporal Resolution Transportation Electric Load Profiles For Multiple Decarbonization Pathways

Electrification of transport compounded with climate change will transform hourly load profiles and their response to weather. We present a novel approach to generating hourly electric load profiles that considers charging strategies and evolving sensitivity to temperature. The approach consists of downscaling annual state-scale sectoral load projections from the multisectoral Global Change Analysis Model (GCAM) into hourly electric load profiles leveraging high resolution climate and population datasets. Profiles are developed and evaluated at the Balancing Authority scale, with a 5-year increment until 2050 over the Western U.S. Interconnect for multiple decarbonization pathways and climate scenarios. The datasets are readily available for production cost model analysis. Our open source approach is transferable to other regions.

Decarbonization, tranportation, Electric vehicle c↗

Effects of Riser Diameter on Solids Holdup and Particle Velocity Profiles in Circulating Fluidized Bed Riser Systems

Tests were performed in a 0.1-m diameter small circulating fluidized bed (SCFB) and 0.3 m diameter cold flow circulating fluidized bed (CFCFB) riser systems located at the National Energy Technology Laboratory (NETL) to study the effects of riser diameter on the riser hydrodynamics. These tests were performed at solids circulation rates of G s = 20 and 75 kg/m 2 s and superficial gas velocities of U g = 5.8 and 6.5 m/s using high-density polyethylene (HDPE) pellets with a density of 0.863 g/cm 3 , particle size range of 600–1400 µm (with a Sauter mean diameter of 871 µm, placing them in the Geldart B classification). Comparisons of riser axial pressure and solids fraction profiles, radial particle velocity profiles, and radial profiles of higher statistical moments and select chaos analysis parameters were considered. The results showed that for a given U g and G s , the smaller diameter riser exhibited characteristics associated with more dilute solids flow than that observed in the larger diameter riser. Additionally, the larger diameter riser exhibited a downward flow of solids near the wall under all test conditions, whereas the smaller diameter riser data exhibited little or no indications of solids downflow near the wall. Overall, these findings suggest that, from an industrial standpoint, a direct scaleup of small-scale tests cannot readily be accomplished as the solids holdup and the solids velocity profiles in small units (those normally tested in the laboratory) are not similar to those of large units and the performance of large units can therefore not be predicted from small-scale tests.

42 ENGINEERING↗

Space-Charge and Other Effects in Fermilab Booster and IOTA Rings’ Ionization Profile Monitors

Ionization profile monitors (IPMs) are widely used in accelerators for non-destructive and fast diagnostics of high energy particle beams. At low beam intensities, initial velocities of the secondaries to collect (ions or electrons) result in the IPM profile smearing. At high beam intensities, the space-charge forces make the measured IPM profiles significantly different from those of the beams. We analyze dynamics of the secondaries in IPMs, describe an effective algorithm to reconstruct the beam sizes from the measured IPM profiles and apply it to the Fermilab 8 GeV proton Booster and 70 MeV/c IOTA ring IPMs.

Shiltsev, Vladimir D.↗

Beam Profile Measurements Utilizing an Amplitude Modulated Pulsed Fiber Laser at PIP2IT

Fermilab is undertaking the development of a new 800 MeV superconducting RF linac to replace its present normal conducting 400 MeV linac. The PIP-II linac consists of a warm front-end generating 2 mA of 2.1 MeV H⁻ followed immediately by a series of superconducting RF cryomodules to 800 MeV. To limit the potential damage to the superconducting RF cavities, PIP-II will utilize laser-based monitors to obtain beam profiles via photoionization. This paper will present the results of transverse and longitudinal beam profile measurements using a prototype profile monitor that was tested with 2.1 MeV H⁻ beam at the PIP-II Injector Test (PIP2IT) accelerator. This prototype profile monitor utilizes a high repetition rate fiber laser and fiber optic transport into the PIP2IT enclosure. In addition, results will be shown of narrow-band electron detection from amplitude modulated laser pulses.

43 PARTICLE ACCELERATORS↗

Simcenter TestLab Error in uploading reference profiles

Our team exports data from Matlab to a Simcenter TestLab formatted .mat file. The .mat file is then uploaded as a reference profile for Environmental Random Control vibration testing. Operations are applied to the data before export that introduce floating point errors. We are experiencing errors uploading the generated reference profiles if the frequency values have a certain level of floating point error (1 Hz vs 1.000000000000001 Hz). I generated some test cases of reference profiles that were giving us trouble, and outlined the errors we were seeing. Two different errors, depending on how the control frequency range is defined. Error is fixed by copy/pasting the frequency column in the Profile Editor.

97 MATHEMATICS AND COMPUTING↗

Retrieved Number Concentration of Cloud Condensation Nuclei (RNCCN) Profile Value-Added Product Report

The cloud condensation nuclei (CCN) concentration at cloud base is the most relevant measure of the aerosol that influences droplet formation in clouds. Since the CCN concentration depends on supersaturation, a more general measure of the CCN concentration is the CCN spectrum (values at multiple supersaturations). The CCN spectrum is now measured at the surface at several U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility observatories and by the ARM Mobile Facility (AMF) but is not measured at the cloud base. Rather than rely on expensive aircraft measurements for all studies of aerosol effects on clouds, a way to project CCN observations made at the surface to cloud base is needed. Remote sensing of aerosol extinction provides information about the vertical profile of the aerosol but cannot be directly related to the CCN concentration because the aerosol extinction is strongly influenced by humidification, particularly near cloud base. Ghan and Collins (2004) and Ghan et al. (2006) propose a method to remove the influence of humidification from the extinction profiles and tie the “dry extinction” retrieval to the surface CCN concentration, thus estimating the CCN profile. This methodology has been implemented as ARM’s Retrieved Number Concentration of CCN (RNCCN) Profile Value-Added Product (VAP).

54 ENVIRONMENTAL SCIENCES↗