Effective R-Matrix Parameterizations for Nuclear Data
Explore the source record for details and available documents.
Arbanas, Goran↗
SEARCH · Engineering Papers
Results for “Approximate boundary”
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.
Natively Periodic Fast Multipole Method: Approximating the Optimal Green’s Function
The fast multipole method (FMM) obeys periodic boundary conditions "natively" if it uses a periodic Green's function for computing the multipole expansion in the interaction zone of each FMM oct-tree node. One can define the "optimal" Green's function for such a method that results in the numerical solution that converges to the equivalent particle-mesh (PM) solution in the limit of sufficiently high order of multipoles. A discrete functional equation for the optimal Green's function can be derived, but is not practically useful as methods for its solution are not known. Instead, this paper presents an approximation for the optimal Green's function that is accurate to better than 10 -3 in ${L}_{\mathrm{MAX}}$ norm and 10 -4 in L 2 norm for practically useful multipole counts. Such an approximately optimal Green's function offers a practical way for implementing the FMM with periodic boundary conditions natively, without the need to compute lattice sums or to rely on hybrid FMM-PM approaches.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
A streamline coordinate analysis of a turbulent boundary layer subject to pressure gradients and curvature on the windward side of a bump
Direct numerical simulation (DNS) of a turbulent boundary layer over the Gaussian (Boeing) bump is performed. This boundary layer exhibits a series of adverse and favourable pressure gradients and convex and concave curvature effects before separating. These effects on turbulent boundary layers are characterised and compared with a lower-Reynolds-number flow over the same geometry. The momentum budgets are analysed in the streamline-aligned coordinate system upstream of the separation region. These momentum budgets allow the simplification of equations to facilitate an integral analysis. Integral-analysis-based approximations for Reynolds stresses in the inner and outer regions of the boundary layer are also formulated. The shear and wall-normal Reynolds stress profiles normalised by these approximations exhibit a better collapse compared with friction velocity and Zagarola–Smits normalisations in the strong favourable pressure gradient region and in the mild adverse pressure region that precedes it in this flow. Simplification of these Reynolds stress approximations along with results from the DNS are used to obtain semi-empirical approximations that are able to provide stress closure in terms of wall solution fields for the turbulent boundary layer under consideration.
42 ENGINEERING↗
Effects of Strain and Strain Rate on Dynamic Grain Growth and Subgrain Evolution During Plastic Deformation of an Interstitial-Free Steel at 850 ° C
Here, the effects of strain and strain rate on dynamic grain growth (DGG) and subgrain evolution are reported for an interstitial-free steel deformed at 850 ° C. Microstructures produced during tension tests at true-strain rates of 10 -4 and to 10 -3 s -1 true strains ranging from 0.02 to 0.2 were preserved following deformation. These were characterized using electron backscatter diffraction (EBSD), including the application of spherical harmonic transform indexing to produce high-angular-resolution EBSD (HR-EBSD) data. HR-EBSD data resolved the small misorientation angles of subgrain boundaries while imaging much larger data fields than possible with previously available techniques. The resulting data confirmed that steady-state flow stress is inversely proportional to the average subgrain size and that subgrain boundary misorientation angle increases with strain. The following new observations are reported. The rate of DGG increased with respect to time but decreased with respect to strain as strain rate increased. This behavior is rationalized through a simple model using separate rate parameters for the effects of time and strain. Subgrain size was not constant during steady-state deformation, but decreased slowly with increasing strain. Subgrain size distributions and subgrain boundary misorientation angle distributions were measured, and both remained approximately log-normal during steady-state deformation. Subgrain evolution demonstrated no dependence on parent grain size, crystallographic orientation, or Taylor factor. These new data suggest that steady-state flow stress is more likely controlled by the dislocation density internal to subgrains than by the spacing between subgrain boundaries.
dynamic grain growth↗
The Spatial and Temporal Variability of the Clear Convective Boundary Layer at the ARM SGP Supersite
The convective boundary layer (CBL), also known as the mixing layer, constitutes the critical lower segment of the atmosphere that significantly influences daily human activities. The growing demand for precise weather forecasts is driven by the requirements of agriculture, transportation, and routine societal functions. Here, to enhance understanding of the CBL, this study investigates the spatiotemporal variability in the CBL and its controlling factors using four-year Doppler lidar, surface flux, and profiling measurements at five ARM Southern Great Plains sites within a 100 km radius. This investigation utilizes data collected exclusively under clear-sky conditions or scattered low-cloud conditions. Results reveal significant spatial differences in CBL evolutions. Daily mixing layer heights (MLHs) vary up to 1 km (30% of the mean) in late afternoon. There is a clear east–west contrast: western sites (C1, E32, E37) exhibit higher summer MLH (1.9–2.1 km) and vertical velocity variances (1.0–1.2 m 2 s −2 ) than eastern sites (1.6–1.8 km), reversing in winter. Temporally, the MLH peaks at 70% of the sunrise–sunset interval, the lagging heat flux (HF) peaks at 50%; and the seasonal MLH maxima lag the HF by approximately one month, influenced by nighttime PBL (planetary boundary layer) properties. The HF and lower tropospheric stability are the main factors of influence for the CBL, but site-specific dependencies highlight the critical roles of local factors, underscoring the need for including them in CBL modeling.
ARM SGP site↗
Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems
In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.
97 MATHEMATICS AND COMPUTING↗
Performance Audit Report: Automated Meteorological Monitoring Stations, Lawrence Livermore National Laboratory
The Waste and Air Quality Offices Group at the Lawrence Livermore National Laboratory operates two main meteorological monitoring sites near Livermore, California. The primary met tower (Site 200) is located in the northwest corner of the lab and consists of a fifty-two (52) meter tower with vertical wind speed sensors, horizontal wind speed and direction sensors and air temperature sensors at fifty-two (52) meters, twenty-three (23) meters and ten (10) meters. In addition, there is a temperature and humidity probe at ten (10) meters. At two (2) meters there are upward and downward facing solar and infrared radiation sensors, an air temperature sensor, barometric pressure sensor and a temperature and relative humidity probe. At ground level there is a tipping bucket rain gauge located to the southwest of the tower. The datalogger, communication peripherals and the barometric pressure sensor are housed in a small building near the tower. The Site 300 met site is located approximately 10 miles east of Livermore within the boundary of the LLNL blast test site. This site also consists of a 52-meter tower which is configured the same as the 52-meter tower at Site 200. There is also a building at this site that houses the datalogger, communication peripherals and the barometric pressure sensor. At both Site 200 and Site 300, data are collected by Ethernet with telephone modem collection as a backup. The data collection interval for both sites is fifteen (15) minutes via workstations running Loggernet software located at the main lab facility. The audits of the Lawrence Livermore National Laboratory’s two automated meteorological stations located near Livermore, California were performed on October 24 & 25, 2023. These audits were performed under the guidelines of the Quality Assurance Handbook for Air Pollution Measurement Systems, Volume IV, Version 2 and the On-Site Meteorological Program Guidance for Regulatory Modeling Applications, U.S. EPA.
54 ENVIRONMENTAL SCIENCES↗
Los Alamos National Laboratory Design Basis Flood Level Modeling Floodplain Mapping Product Report
This AOI covers a land area of approximately 67 square miles. The seven sitewide watershed boundaries are shown in Figure 1, which also shows the hydro-enforced DEM rendered as a shaded relief image.
58 GEOSCIENCES↗
Aerosol size distribution properties associated with cold-air outbreaks in the Norwegian Arctic
The aerosol particles serving as cloud condensation and ice nuclei contribute to key cloud processes associated with cold-air outbreak (CAO) events but are poorly constrained in climate models due to sparse observations. Here we retrieve aerosol number size distribution modes from measurements at Andenes, Norway, during the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) and at Zeppelin Observatory, approximately 1000 km upwind from Andenes at Svalbard. During CAO events at Andenes, the sea-spray-mode number concentration is correlated with strong over-ocean winds with a mean of 8±4 cm −3 that is 71 % higher than during non-CAO conditions. Additionally, during CAO events at Andenes, the mean Hoppel minimum diameter is 6 nm smaller than during non-CAO conditions, though the estimated supersaturation is lower, and the mean number concentration of particles that likely activated in-cloud is 109±61 cm −3 with no statistically significant difference from the non-CAO mean of 99±66 cm −3 . For CAO trajectories between Zeppelin Observatory and Andenes, the upwind-to-downwind change in number concentration is the largest for the accumulation mode with a mean decrease of 93±95 cm −3 , likely attributable primarily to precipitation scavenging. These characteristic properties of aerosol number size distributions during CAO events provide guidance for evaluating CAO aerosol–cloud interaction processes in models.
54 ENVIRONMENTAL SCIENCES↗
Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination
We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.
Lopez, Graham↗
A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields
We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Ultrasonic waves from radial mode excitation of a piezoelectric disc on the surface of an elastic solid
In this work, we discuss the ultrasonic wave propagation characteristics in an isotropic elastic solid due to the radial mode excitation of a piezoelectric disc actuator bonded to its surface. Finite element simulations using coupled electromechanical modeling is employed to investigate the wave propagation behavior. We find that the radial mode vibrations on the surface of the elastic solid generates all the three types of ultrasonic waves: longitudinal, shear, and surface waves. The waves in the solid are comprised of a central lobe and multiple side lobes based on the frequency of the radial mode excitation. The central lobe is predominantly composed of longitudinal waves and the side lobes are composed of shear waves. While the longitudinal waves have a strong central lobe and weak side lobes that are not fully developed, shear waves consists of fully developed side lobes. In addition, we observe that the longitudinal waves have fewer side lobes inside the solid compared to the shear waves. A semi-analytical approach is presented to explain the above observations. The interface traction boundary condition between the piezoelectric disc and the elastic solid is approximated as a truncated 'Bessel' excitation over the area of the piezoelectric disc in contact with the solid. A Fourier–Bessel series expansion technique is used to investigate the wave propagation behavior from the above traction boundary condition. The results obtained explain the observations from the finite element simulations with regard to the number and distribution of side lobes pertaining to the longitudinal and shear waves generated from the radial modes of the piezoelectric disc. The proposed semi-analytical approach is general and can be applied to any arbitrary axisymmetric excitations on the surface of an elastic solid.
36 MATERIALS SCIENCE↗
Efficient Generalized Boundary Detection Using a Sliding Information Distance
In this work, we present a general machine learning algorithm for boundary detection within general signals based on an efficient, accurate, and robust approximation of the universal normalized information distance. Our approach uses an adaptive sliding information distance (SLID) combined with a wavelet-based approach for peak identification to locate the boundaries. Special emphasis is placed on developing an adaptive formulation of SLID to handle general signals with multiple unknown and/or drifting section lengths. Although specialized algorithms may outperform SLID when domain knowledge is available, these algorithms are limited to specific applications and do not generalize. SLID excels in these cases. We demonstrate the versatility and efficacy of SLID on a variety of signal types, including synthetically generated sequences of tokens, binary executables for reverse engineering applications, and time series of seismic events.
42 ENGINEERING↗
On the prescription of boundary conditions for nonlocal Poisson's and peridynamics models
We introduce a technique to automatically convert local boundary conditions into nonlocal volume constraints for nonlocal Poisson’s and peridynamic models. The proposed strategy is based on the approximation of nonlocal Dirichlet or Neumann data with a local solution obtained by using available boundary, local data. The corresponding nonlocal solution converges quadratically to the local solution as the nonlocal horizon vanishes, making the proposed technique asymptotically compatible. The proposed conversion method does not have any geometry or dimensionality constraints and its computational cost is negligible, compared to the numerical solution of the nonlocal equation. The consistency of the method and its quadratic convergence with respect to the horizon is illustrated by several two-dimensional numerical experiments conducted by meshfree discretization for both the Poisson’s problem and the linear peridynamic solid model.
97 MATHEMATICS AND COMPUTING↗
Open cells exhibit weaker entrainment of free-tropospheric biomass burning aerosol into the south-east Atlantic boundary layer
This work presents synergistic satellite, airborne and surface-based observations of a pocket of open cells (POC) in the remote south-east Atlantic. The observations were obtained over and upwind of Ascension Island during the CLouds and Aerosol Radiative Impacts and Forcing (CLARIFY) and the Layered Smoke Interacting with Clouds (LASIC) field experiments. A novel aspect of this case study is that an extensive free-tropospheric biomass burning aerosol plume that had been transported from the African continent was observed to be in contact with the boundary layer inversion over the POC and the surrounding closed cellular cloud regime. The in situ measurements show marked contrasts in the boundary layer thermodynamic structure, cloud properties, precipitation and aerosol conditions between the open cells and surrounding overcast cloud field. The data demonstrate that the overlying biomass burning aerosol was mixing down into the boundary layer in the stratocumulus cloud downwind of the POC, with elevated carbon monoxide, black carbon mass loadings and accumulation-mode aerosol concentrations measured beneath the trade-wind inversion. The stratocumulus cloud in this region was moderately polluted and exhibited very little precipitation falling below cloud base. A rapid transition to actively precipitating cumulus clouds and detrained stratiform remnants in the form of thin quiescent veil clouds was observed across the boundary into and deep within the POC. The subcloud layer in the POC was much cleaner than that in the stratocumulus region. The clouds in the POC formed within an ultra-clean layer (accumulation-mode aerosol concentrations of approximately a few cm -3 ) in the upper region of the boundary layer, which was likely to have been formed via efficient collision–coalescence and sedimentation processes. Enhanced Aitken-mode aerosol concentrations were also observed intermittently in this ultra-clean layer, suggesting that new particle formation was taking place. Across the boundary layer inversion and immediately above the ultra-clean layer, accumulation-mode aerosol concentrations were ~1000 cm -3 . Importantly, the air mass in the POC showed no evidence of elevated carbon monoxide over and above typical background conditions at this location and time of year. As carbon monoxide is a good tracer for biomass burning aerosol that is not readily removed by cloud processing and precipitation, it demonstrates that the open cellular convection in the POC is not able to entrain large quantities of the free-tropospheric aerosol that was sitting directly on top of the boundary layer inversion. This suggests that the structure of the mesoscale cellular convection may play an important role in regulating the transport of aerosol from the free troposphere down into the marine boundary layer. We then develop a climatology of open cellular cloud conditions in the south-east Atlantic from 19 years of September Moderate Resolution Imaging Spectroradiometer (MODIS) Terra imagery. This shows that the maxima in open cell frequency (> 0.25) occurs far offshore and in a region where subsiding biomass burning aerosol plumes may often come into contact with the underlying boundary layer cloud. If the results from the observational case study applied more broadly, then the apparent low susceptibility of open cells to free-tropospheric intrusions of additional cloud condensation nuclei could have some important consequences for aerosol–cloud interactions in the region.
54 ENVIRONMENTAL SCIENCES↗
A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics
The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.
54 ENVIRONMENTAL SCIENCES↗
Breakdown of the $v=1$ integer quantum Hall effect in the single particle and collective localization regimes
We report large signal current–voltage characteristics measurements near the $v=1$ integer quantum Hall plateau in a high mobility GaAs/AlGaAs sample that exhibits both single particle and collective localization. According to prior results, the central part of the $v=1$ plateau corresponds to single particle localization and hence an Anderson insulator. In contrast, the two regions in the flanks of this Anderson insulator were associated with collective localization and were identified with the integer quantum Hall Wigner solid. We find that, at low temperatures and under large bias currents, both the Anderson insulator and integer quantum Hall Wigner solid have a sharp onset of conduction above a threshold. Furthermore, this threshold shows an approximately monotonic dependence on the filling factor across the boundary between the two localized phases. Finally, to our surprise, we find that these large signal characteristics do not clearly distinguish between the Anderson insulator and the integer quantum Hall Wigner solid.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Complexion dictated thermal resistance with interface density in reactive metal multilayers
Multilayers composed of aluminum (A1) and platinum (Pt) exhibit a non-monotonic trend in thermal resistance with bilayer thickness as measured by time domain thermoreflectance. The thermal resistance initially increases with reduced bilayer thickness only to reach a maximum and then decrease with further shrinking of the multilayer period. These observations are attributed to the evolving impact of an intermixed amorphous complexion approximately 10 nm in thickness, which forms at each boundary between Al- and Pt-rich layers. Scanning transmission electron microscopy combined with energy dispersive X-ray spectroscopy find that the elemental composition of the complexion varies based on bilayer periodicity as does the fraction of the multilayer composed of this interlayer. These variations in complexion mitigate boundary scattering within the multilayers as shown by electronic transport calculations employing density functional theory and non-equilibrium Green's functions on amorphous structures obtained via finite temperature molecular dynamics. The lessening of boundary scattering reduces the total resistance to thermal transport leading to the observed non-monotonic trend thereby highlighting the central role of complexion on thermal transport within reactive metal multilayers.
36 MATERIALS SCIENCE↗