Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal modes”

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 127 records · Page 7

Application of Monte Carlo code SHIFT for MSBR dose rate calculations

This paper presents a study that applied the Monte Carlo code SHIFT to calculate the radiological environment within a detailed Molten Salt Breeder Reactor (MSBR) model. It represents one of the first applications of the SHIFT code in stand-alone simulations for non-light water reactors to help demonstrate the code's potential uses in the design, licensing, and operation of advanced reactors. The radiological conditions of the MSBR were modeled when the reactor is at two different operation modes: normal full power and a drained state. The Forward Weighted-Consistent Adjoint Driven Importance Sampling (FW-CADIS) hybrid method in SHIFT was applied successfully to calculate the ex-core neutron and gamma dose rates for the MSBR at full power. Neutron and gamma dose rates within the drain cell were also calculated for the MSBR at the drained state by integrating the source terms obtained from an ORIGEN-S depletion calculation into the SHIFT simulation. The results indicate that in the drained state, the delayed gammas from the depleted fuel salt are the main contributors to the dose rates. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Accelerating full-waveform inversion using source stacking: synthetic experiments at the global scale in a realistic 3-D earth model

SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.

Geochemistry & Geophysics↗

Isotope effects on transport in LHD

Abstract Isotope effects are one of the most important issues for predicting future reactor operations. Large helical device (LHD) is the presently working largest stellarator/helical device using super conducting helical coils. In LHD, deuterium experiments started in 2017. Extensive studies regarding isotope effects on transport have been carried out. In this paper, the results of isotope effect studies in LHD are reported. The systematic studies were performed adjusting operational parameters and nondimensional parameters. In L mode like normal confinement plasma, where internal and edge transport barriers are not formed, the scaling of global energy confinement time ( τ E ) with operational parameters shows positive mass dependence ( M 0.27 ; where M is effective ion mass) in electron cyclotron heating plasma and no mass dependence ( M 0.0 ) in neutral beam injection heating plasma. The non-negative ion mass dependence is anti-gyro-Bohm scaling. The role of the turbulence in isotope effects was also found by turbulence measurements and gyrokinetic simulation. Better accessibility to electron and ion internal transport barrier (ITB) plasma is found in deuterium (D) plasma than in hydrogen (H). Gyro kinetic non-linear simulation shows reduced ion heat flux due to the larger generation of zonal flow in deuterium plasma. Peaked carbon density profile plays a prominent role in reducing ion energy transport in ITB plasma. This is evident only in plasma with deuterium ions. New findings on the mixing and non-mixing states of D and H particle transports are reported. In the mixing state, ion particle diffusivities are higher than electron particle diffusivities and D and H ion density profiles are almost identical. In the non-mixing state, ion particle diffusivity is much lower than electron diffusivity. Deuterium and hydrogen ion profiles are clearly different. Different turbulence structures were found in the mixing and non-mixing states suggesting different turbulence modes play a role.

Physics↗

A cosine-based correlation information entropy approach for building automatic fault detection baseline construction

Building automatic fault detection and diagnosis (AFDD) technologies have shown great potential for energy savings. To enable AFDD, a baseline depicting the normal operation mode is needed to detect whether the building operation deviates from normality. Existing research using physics-based knowledge and models for AFDD has mainly taken a trial-and-error approach to determine if a given baseline is sufficient via empirical experiments. A mechanism to support decisions such as how many samples and what samples should be included in the baseline is currently lacking. In this study, a data-driven method for AFDD baseline construction based on information entropy is developed. The entropy is derived based on cosine similarity among typical building automation system measurements in conjunction with outdoor weather information. The performance of the proposed method is evaluated using real building data. Evaluation results indicate that the fault detection strategy adopting the proposed method has similar or better accuracy in detecting faults compared to the same fault detection strategy using the baseline construction method from the literature. Additionally, the use of entropy enables the proposed method to automatically construct and assess the baseline consisting of information-rich samples.

42 ENGINEERING↗

A Stable and Ultrafast Control for Multiparallel Grid-Tied SiC Inverters With LVRT Capability Considering Cable Impedance

The stability of multiparalleled grid-tied inverters remains a challenging issue especially during low-voltage ride-through (LVRT) in the presence of grid impedance. SiC device-based grid-tied inverters can achieve advantages such as higher efficiency, reduced size of interface filters, and fast dynamics; however, they can exacerbate this stability issue particularly when LVRT occurs along with high inrush current due to their low inertia characteristics. In this article, a direct deadbeat control method without any PI control loop is proposed for multiparallel grid-tied SiC inverters. The proposed control achieves instantaneous current regulations within switching cycles. In addition, it mitigates high inrush currents and maintains stability during LVRT transients without requiring mode transitions between normal and LVRT operation modes, thereby achieving an ultrafast response and seamless operation through this unified control method. The Lyapunov method is applied to analyze the multiparallel inverter stability in the presence of grid impedance with proposed deadbeat control method. The inverter V-I trajectories during LVRT transients are derived to illustrate the ultrafast and stable feature of proposed control. Finally, the experimental results of two grid-tied SiC inverters operating at LVRT transients are provided to verify the advantages of the proposed control method.

14 SOLAR ENERGY↗

High Free-Tropospheric Aitken-Mode Aerosol Concentrations Buffer Cloud Droplet Concentrations in Large-Eddy Simulations of Precipitating Stratocumulus

A new Aitken mode aerosol microphysics scheme is developed for a large eddy simulation model in order to better investigate cloud-aerosol interactions in the marine boundary layer and to study the Aitken buffering hypothesis of McCoy et al. (2021), https://doi.org/10.1029/2020jd033529. This scheme extends the single-mode two-moment prognostic aerosol scheme of Berner et al. (2013), https://doi.org/10.5194/acp-13-12549-2013. Seven prognostic variables represent accumulation and Aitken log-normal aerosol modes in air and droplets as well as 3 gas species. Scavenging of interstitial and other unactivated aerosol by cloud and rain drops are treated using the scheme described in Berner et al. (2013), https://doi.org/10.5194/acp-13-12549-2013. The scheme includes coagulation of unactivated aerosol and a simple chemistry model with gas phase H 2 SO 4 , SO 2 , and DMS as prognostic variables to capture basic influences of sulfur chemistry on the model aerosols. Nucleation of H 2 SO 4 aerosol particles from gas-phase H 2 SO 4 is neglected. A deep, precipitating stratocumulus case (VAMOS Ocean Cloud Atmosphere Land Study RF06) is used to test the new scheme. The presence of the Aitken mode aerosol increases the cloud droplet concentration through activation of the larger Aitken particles and delays the creation of an ultraclean, strongly precipitating cumulus state. Scavenging of unactivated accumulation and Aitken particles by cloud and precipitation droplets accelerates the collapse. Increasing either the above-inversion Aitken concentration or the surface Aitken flux increases the Aitken population in the boundary layer and prevents the transition to an ultraclean state.

54 ENVIRONMENTAL SCIENCES↗

Drop impact tester development for spent nuclear fuel vibration integrity study

The accumulated damage of SNF rod is defined as the integrated damage induced by the SNF assembly harmonic vibrations and SNF assembly & cask system contact interactions induced transient shocks, experienced by SNF rods during normal conditions of transportation (NCT). The CIRFT technology has provided means to resolve/understand the SNF rod harmonic vibration reliability issues associated with maximum axial clad tubing stress induced by SNF inertia related dynamic bending loading. However, the SNF rod integrity in association with the transient shock loads during NCT through dynamic contact interactions among SNF assembly cannot be investigated systematically through CIRFT technology alone. Thus, the proposed drop impact tester development become a viable tool for SNF system transient shock effect investigation. Furthermore, the transient shocks induced loading mode is normally perpendicular to the SNF rod axial orientation, i.e., in a transverse orientation to the SNF rod. To proceed with such loading mode investigation, a new device, “SNF Drop Impact Tester” that provides the transverse impact shock load onto a SNF rod through a free-drop projectile was developed in FY20 for SNF System Vibration Reliability Investigation; where the SNF system contact impact intensity can be calibrated accordingly with projectile weight and projectile head geometry. This newly developed research tool can further provide a detailed understanding about the effect of dynamic contact-interaction loadings; and combined harmonic vibration and transient shock loading modes on the fatigue damage evolution of the HBU SNF under NCT.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

HIPED: Machine learning framework for spherical tokamak pedestal prediction and optimization

We introduce a Machine Learning framework, HIPED (HeIght and width Predictor for Edge Dynamics), for predicting and optimizing pedestal and core performance in spherical tokamak plasmas. Trained on pedestal and core datasets from the third MAST-U campaign, HIPED provides accurate estimates of pedestal height and width. The results reveal notable differences compared with conventional aspect-ratio studies; for instance, a simple power-law relation between pedestal width and height has very low accuracy. Instead, additional parameters such as normalized plasma pressure, elongation, and Greenwald fraction significantly improve accuracy. HIPED can also be trained only on `control room parameters' to inform experimentalists of which controllable parameters to adjust for improving core-integrated performance. The framework further includes a multi-objective optimization scheme that helps guide experimental planning and optimization. We find Pareto-optimal discharges with respect to various features, including distance from edge-localized modes and normalized plasma pressure, track their parameter trajectories over time, and identify the control room parameters required for these Pareto-optimal discharges. This provides a framework for systematically optimizing core and edge performance according to different experimental priorities.

Parisi, Jason F. [Princeton Plasma Physics Laborat↗

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quasiclassical sampling and Wigner sampling of initial vibrational coordinates and momenta for polyatomic molecules in Monte Carlo molecular dynamics simulations

In a quasiclassical trajectory simulation, the vibrational modes are initialised with quantised vibrational energies, but vibrational phases are sampled by Monte Carlo. This requires an algorithm to assign coordinates and momenta to the various atoms. In this work, we present two methods for implementing this for nonrotating polyatomic molecules, namely, fixed-energy vibrational-state-selected initial conditions and thermal initial conditions. We also present a method for initiating classical trajectories with a ground-state Wigner distribution. These vibrational treatments are sufficient to initialise trajectories for unimolecular processes, and we also show how they can be applied to simulate bimolecular collision processes. The treatments of unimolecular and bimolecular collision processes are available in two Python codes called wigner_state_selected.py and bimolecular_collision.py, respectively, which will generate initial condition files that are recognisable by the SHARC and SHARC-MN computer programs for dynamics calculations. Both codes are available as standalone programs, as well as being included in SHARC-MN, and they will be included in future versions of SHARC. Here, the methods implemented in these codes are mostly also available in the ANT computer program, and those that are not available in ANT will be incorporated in future versions of ANT.

Wigner distribution↗

Experimental study of a thermal energy storage-integrated heat pump system for load shifting in space cooling and heating

This study experimentally investigates an all-season thermal energy storage-integrated heat pump (TES-HP) system developed to enhance building energy efficiency and support grid-interactive operation through load shifting in both cooling and heating. A 4-ton commercial rooftop air-source heat pump was modified by integrating a hydronic thermal energy storage (TES) unit containing a phase change material (PCM) with a melting temperature of 22.0 °C. The system employs two reversing valves and three electronic expansion valves to enable six operating modes, including normal, TES charging, and TES discharging in both cooling and heating seasons. A subcooling-controlled electronic expansion valve was implemented to mitigate refrigerant maldistribution caused by unequal internal volumes among the heat exchangers. Compared with a baseline heat pump, the TES-HP reduced power consumption by 30–50% during cooling at an ambient temperature of 40.6 °C, and by up to 60% during heating at −15.0 °C, while maintaining high performance under cold-climate conditions. Additionaly to experimental evaluation, a simplified annualized on-peak analysis was conducted to compare the TES-HP with the baseline system, indicating an on-peak electricity reduction of approximately 876 kWhₑ per unit per year and associated on-peak cost savings under time-of-use pricing. These results demonstrate that the TES-HP effectively decouples HP operation from adverse ambient conditions, improving flexibility, resilience, and energy efficiency. This work advances prior research by experimentally proving the integration of PCM-based TES for a year-round, load-flexible HVAC operation.

42 ENGINEERING↗

Effects of Intercalating Molecules on the Polymer Properties of DNA

Atomic force microscopy (AFM) enables determination of physical properties from single DNA molecules. Insertion of aromatic molecules into the structure of DNA results in morphological changes. However, the accompanying changes to elastic properties due to this insertion are not fully understood. AFM was used to examine the morphological effects of intercalator binding and report changes in the elastic properties of intrinsically straight DNA molecules. The persistence length and polymer extension were characterized in the presence of three intercalating molecules: ethidium bromide and the less well studied chloroquine and acridine. It was found that all three intercalators significantly increased the bending persistence length. In addition, an analysis of the normal bending modes of the static molecules corroborated these results. Furthermore, this approach of measuring binding effects of intercalators on DNA physical properties using a model system of intrinsically straight DNA is applicable to other DNA binding ligands and other modes of DNA interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The dynamical tides of spinning Newtonian stars

ABSTRACT We carefully develop the framework required to model the dynamical tidal response of a spinning neutron star in an inspiralling binary system, in the context of Newtonian gravity, making sure to include all relevant details and connections to the existing literature. The tidal perturbation is decomposed in terms of the normal oscillation modes, used to derive an expression for the effective Love number which is valid for any rotation rate. In contrast to previous work on the problem, our analysis highlights subtle issues relating to the orthogonality condition required for the mode-sum representation of the dynamical tide and shows how the prograde and retrograde modes combine to provide the overall tidal response. Utilizing a slow-rotation expansion, we show that the dynamical tide (the effective Love number) is corrected at first order in rotation, whereas in the case of the static tide (the static Love number) the rotational corrections do not enter until second order.

Pnigouras, P. (ORCID:0000000318959431)↗

Eigenvalue-Based Micromagnetic Analysis of Switching in Spin-Torque-Driven Structures

We present an eigenvalue-based approach for studying the magnetization dynamics in magnetic nanostructures driven by spintronic excitations, such as spin-transfer torque and spin-orbit torque. The approach represents the system dynamics in terms of normal oscillation modes (eigenstates) with corresponding complex eigenfrequencies. The dynamics is driven by a small number of active eigenstates and often considering just a single eigenstate is sufficient. We develop a perturbation theory that provides semianalytical dynamic solutions by using eigenstates for the case in the absence of damping and spintronic excitations as a basis. The approach provides useful insights into dynamics in such systems and allows solving several difficulties in their modeling, such as extracting the switching current in magnetic random-access memories and understanding switching mechanisms. We show that the presented approach directly predicts the critical switching current, i.e., switching current for an infinite time. The approach also provides solutions for the switching dynamics allowing the switching current to be obtained for a finite switching time, provided that the system symmetry is broken, e.g., by tilting the polarizer, so that switching by a finite pulse is possible.

42 ENGINEERING↗

Evidence of Completion of Milestone 3: Experimental Testbed

Construction and deployment of the experimental testbed encountered scheduling delays, and completion was subsequently further delayed by curtailment of LANL operations due to the COVID-19 pandemic. Between March and October, LANL was in a state of “Limited Operations” such that only certain “mission-critical” work was performed with approximately 25% on-site staffing. In spite of this, much of the lab space preparation was executed and the five initial model house” units were constructed and installed during this period. Still, the pace of construction was impacted by difficulty in coordinating personnel and minimizing contact between workers. Currently, LANL is in the mode of “Normal Operations with Maximum Telework,” which continues to limit the availability of on-site personnel. Nevertheless, Milestone 3 is now complete. This document describes the experiment as it is deployed and provides status for each task for this milestone.

99 GENERAL AND MISCELLANEOUS↗