Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Methods”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Computational Fluid Dynamics Prediction of a Modified Savonius Wind Turbine with Novel Blade Shapes

The Savonius wind turbine is a type of vertical axis wind turbine (VAWTs) that is simply composed of two or three arc-type blades which can generate power even under poor wind conditions. A modified Savonius wind turbine with novel blade shapes is introduced with the aim of increasing the power coefficient of the turbine. The effect of blade fullness, which is a main shape parameter of the blade, on the power production of a two-bladed Savonius wind turbine is investigated using transient computational fluid dynamics (CFD). Simulations are based on the Reynolds Averaged Navier-Stokes (RANS) equations with a renormalization group turbulent model. This numerical method is validated with existing experimental data and then utilized to quantify the performance of design variants. Results quantify the relationship between blade fullness and turbine performance with a blade fullness of 1 resulting in the highest coefficient of power, 0.2573. This power coefficient is 10.98% higher than a conventional Savonius turbine.

17 WIND ENERGY↗
DOE OSTI2024

The Effect of Interlayer Delay on the Heat Accumulation, Microstructures, and Properties in Laser Hot Wire Directed Energy Deposition of Ti-6Al-4V Single-Wall

Laser hot wire directed energy deposition (LHW-DED) is a layer-by-layer additive manufacturing technique that permits the fabrication of large-scale Ti-6Al-4V (Ti64) components with a high deposition rate and has gained traction in the aerospace sector in recent years. However, one of the major challenges in LHW-DED Ti64 is heat accumulation, which affects the part quality, microstructure, and properties of as-built specimens. These issues require a comprehensive understanding of the layerwise heat-accumulation-driven process–structure–property relationship in as-deposited samples. In this study, a systematic investigation was performed by fabricating three Ti-6Al-4V single-wall specimens with distinct interlayer delays, i.e., 0, 120, and 300 s. The real-time acquisition of high-fidelity thermal data and high-resolution melt pool images were utilized to demonstrate a direct correlation between layerwise heat accumulation and melt pool dimensions. The results revealed that the maximum heat buildup temperature of the topmost layer decreased from 660 °C to 263 °C with an increase to a 300 s interlayer delay, allowing for better control of the melt pool dimensions, which then resulted in improved part accuracy. Furthermore, the investigation of the location-specific composition, microstructure, and mechanical properties demonstrated that heat buildup resulted in the coarsening of microstructures and, consequently, the reduction of micro-hardness with increasing height. Extending the delay by 120 s resulted in a 5% improvement in the mechanical properties, including an increase in the yield strength from 817 MPa to 859 MPa and the ultimate tensile strength from 914 MPa to 959 MPa. Cooling rates estimated at 900 °C using a one-dimensional thermal model based on a numerical method allowed us to establish the process–structure–property relationship for the wall specimens. The study provides deeper insight into the effect of heat buildup in LHW-DED and serves as a guide for tailoring the properties of as-deposited specimens by regulating interlayer delay.

36 MATERIALS SCIENCE↗
DOE OSTI2016

POLARIZED SYNCHROTRON EMISSIVITIES AND ABSORPTIVITIES FOR RELATIVISTIC THERMAL, POWER-LAW, AND KAPPA DISTRIBUTION FUNCTIONS

Synchrotron emission and absorption determine the observational appearances of many astronomical systems. In this paper, we describe a numerical scheme for calculating synchrotron emissivities and absorptivities in all four Stokes parameters for arbitrary gyrotropic electron distribution functions, building on earlier work by Leung, Gammie, and Noble. We use this technique to evaluate the emissivities and the absorptivities for a thermal (Maxwell–Jüttner), isotropic power-law, and an isotropic kappa distribution function. The latter contains a power-law tail at high particle energies that smoothly merges with a thermal core at low energies, as is characteristic of observed particle spectra in collisionless plasmas. We provide fitting formulae and error bounds on the fitting formulae for use in codes that solve the radiative transfer equation. The numerical method and the fitting formulae are implemented in a compact C library called symphony. We find that the kappa distribution has a source function that is indistinguishable from a thermal spectrum at low frequency and transitions to the characteristic self-absorbed synchrotron spectrum, ∝ν{sup 5/2}, at high frequency; the linear polarization fraction for a thermal spectrum is near unity at high frequency; and all distributions produce O(10%) circular polarization at low frequency for lines of sight sufficiently close to the magnetic field vector.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Formal Solutions for Polarized Radiative Transfer. I. The DELO Family

The discussion regarding the numerical integration of the polarized radiative transfer equation is still open and the comparison between the different numerical schemes proposed by different authors in the past is not fully clear. Aiming at facilitating the comprehension of the advantages and drawbacks of the different formal solvers, this work presents a reference paradigm for their characterization based on the concepts of order of accuracy , stability , and computational cost . Special attention is paid to understand the numerical methods belonging to the Diagonal Element Lambda Operator family, in an attempt to highlight their specificities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Emulating Simulations of Cosmic Dawn for 21 cm Power Spectrum Constraints on Cosmology, Reionization, and X-Ray Heating

Current and upcoming radio interferometric experiments are aiming to make a statistical characterization of the high-redshift 21 cm fluctuation signal spanning the hydrogen reionization and X-ray heating epochs of the universe. However, connecting 21 cm statistics to the underlying physical parameters is complicated by the theoretical challenge of modeling the relevant physics at computational speeds quick enough to enable exploration of the high-dimensional and weakly constrained parameter space. In this work, we use machine learning algorithms to build a fast emulator that can accurately mimic an expensive simulation of the 21 cm signal across a wide parameter space. We embed our emulator within a Markov Chain Monte Carlo framework in order to perform Bayesian parameter constraints over a large number of model parameters, including those that govern the Epoch of Reionization, the Epoch of X-ray Heating, and cosmology. As a worked example, we use our emulator to present an updated parameter constraint forecast for the Hydrogen Epoch of Reionization Array experiment, showing that its characterization of a fiducial 21 cm power spectrum will considerably narrow the allowed parameter space of reionization and heating parameters, and could help strengthen Planck's constraints on σ{sub 8}. We provide both our generalized emulator code and its implementation specifically for 21 cm parameter constraints as publicly available software.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Constructing Current Singularity in a 3D Line-tied Plasma

We revisit Parker's conjecture of current singularity formation in 3D line-tied plasmas using a recently developed numerical method, variational integration for ideal magnetohydrodynamics in Lagrangian labeling. With the frozen-in equation built-in, the method is free of artificial reconnection, and hence it is arguably an optimal tool for studying current singularity formation. Using this method, the formation of current singularity has previously been confirmed in the Hahm–Kulsrud–Taylor problem in 2D. In this paper, we extend this problem to 3D line-tied geometry. The linear solution, which is singular in 2D, is found to be smooth for arbitrary system length. However, with finite amplitude, the linear solution can become pathological when the system is sufficiently long. The nonlinear solutions turn out to be smooth for short systems. Nonetheless, the scaling of peak current density versus system length suggests that the nonlinear solution may become singular at finite length. Finally, with the results in hand, we can neither confirm nor rule out this possibility conclusively, since we cannot obtain solutions with system length near the extrapolated critical value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

The Aemulus Project. III. Emulation of the Galaxy Correlation Function

Using the N-body simulations of the Aemulus Project, we construct an emulator for the nonlinear clustering of galaxies in real and redshift space. We construct our model of galaxy bias using the halo occupation framework, accounting for possible velocity bias. The model includes 15 parameters, including both cosmological and galaxy bias parameters. We demonstrate that our emulator achieves ~1% precision at the scales of interest, 0.1 h -1 Mpc < r < 10 h -1 Mpc, and recovers the true cosmology when tested against independent simulations. Our primary parameters of interest are related to the growth rate of structure, f, and its degenerate combination, fσ 8 . Using this emulator, we show that the constraining power on these parameters monotonically increases as smaller scales are included in the analysis, all the way down to 0.1 h -1 Mpc. For a BOSS-like survey, the constraints on fσ 8 from r < 30 h -1 Mpc scales alone are nearly a factor of two tighter than those from the fiducial BOSS analysis of redshift-space clustering using perturbation theory at larger scales. The combination of real- and redshift-space clustering allows us to break the degeneracy between f and σ 8, yielding an 11% constraint on f alone for a BOSS-like analysis. The current Aemulus simulations limit this model to surveys of massive galaxies. Future simulations will allow this framework to be extended to all galaxy target types, including emission-line galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

On the Growth and Saturation of the Gyroresonant Streaming Instabilities

The self-regulation of cosmic-ray (CR) transport in the interstellar and intracluster media has long been viewed through the lenses of linear and quasi-linear kinetic plasma physics. Such theories are believed to capture the essence of CR behavior in the presence of self-generated turbulence but cannot describe potentially critical details arising from the nonlinearities of the problem. We utilize the particle-in-cell numerical method to study the time-dependent nonlinear behavior of the gyroresonant streaming instabilities, self-consistently following the combined evolution of particle distributions and self-generated wave spectra in one-dimensional periodic simulations. We demonstrate that the early growth of instability conforms to the predictions from linear physics, but that the late-time behavior can vary depending on the properties of the initial CR distribution. We emphasize that the nonlinear stages of instability depend strongly on the initial anisotropy of CRs—highly anisotropic CR distributions do not efficiently reduce to Alfvénic drift velocities, owing to reduced production of left-handed resonant modes. We derive estimates for the wave amplitudes at saturation and the timescales for nonlinear relaxation of the CR distribution and then demonstrate the applicability of these estimates to our simulations. Bulk flows of the background plasma due to the presence of resonant waves are observed in our simulations, confirming the microphysical basis of CR-driven winds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Multi-frequency General Relativistic Radiation-hydrodynamics with M 1 Closure

We report on recent upgrades to our general relativistic radiation-magnetohydrodynamics code, Cosmos++ , which expands the two-moment, M 1 , radiation treatment from gray to multi-frequency transport, including Doppler and gravitational frequency shifts. The solver accommodates either photon (Bose–Einstein) or neutrino (Fermi–Dirac) statistical distribution functions with absorption, emission, and elastic scattering processes. An implicit scheme is implemented to simultaneously solve the primitive inversion problem together with the radiation–matter coupling source terms, providing stability over a broad range of opacities and optical depths where the interaction terms can be stiff. We discuss our formulations and numerical methods, and validate our methods against a wide variety of test problems spanning optically thin to thick regimes in flat, weakly curved, and strongly curved spacetimes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Mira-Titan Universe. III. Emulation of the Halo Mass Function

We construct an emulator for the halo mass function over group and cluster mass scales for a range of cosmologies, including the effects of dynamical dark energy and massive neutrinos. The emulator is based on the recently completed Mira-Titan Universe suite of cosmological N-body simulations. The main set of simulations spans 111 cosmological models with 2.1 Gpc boxes. We extract halo catalogs in the redshift range z = [0.0, 2.0] and for masses M-200c >= 10(13)M(circle dot)/h. The emulator covers an eight-dimensional hypercube spanned by {Omega(m)h(2), Omega(b)h(2), Omega(nu)h(2), sigma(8), h, n(s), w(0), w(a)}; spatial flatness is assumed. We obtain smooth halo mass functions by fitting piecewise second-order polynomials to the halo catalogs and employ Gaussian process regression to construct the emulator while keeping track of the statistical noise in the input halo catalogs and uncertainties in the regression process. For redshifts z less than or similar to 1, the typical emulator precision is better than 2% for 10(13)-10(14)M(circle dot)/h and <10<^> M similar or equal to 101(circle dot)(54)/h. For comparison, fitting functions using the traditional universal form for the halo mass function can be biased at up to 30% at M similar or equal to 10(144)M(circle dot)/h for z = 0. Our emulator is publicly available at https://github.com/SebastianBocquet/MiraTitanHMFemulator.

cosmology: theory↗
DOE OSTI2021

Fermi-type particle acceleration from magnetic reconnection at the termination shock of a relativistic striped wind

An oblique-rotating pulsar generates a relativistic striped wind in a pulsar wind nebula (PWN). The termination shock of the PWN compresses the Poynting-flux-dominated flow and drives magnetic reconnection. By carrying out particle-in-cell simulations of the termination shock of the PWN, we study the shock structure as well as the energy conversion processes and particle acceleration mechanisms. With the recent advances in the numerical methods, we extend the simulations to the ultrarelativistic regime with a bulk Lorentz factor of up to γ 0 = 10 6 . Magnetic reconnection at the termination shock is highly efficient at converting magnetic energy to particle kinetic energy and accelerating particles to high energies. Similar to earlier studies, we find that the resulting energy spectra crucially depend on λ/d e (λ is the wavelength of the striped wind and de is the relativistic plasma skin depth). When λ/d e is large (λ ≳ 40d e ), the downstream particle spectra form a power-law distribution in the magnetically dominated relativistic wind regime. By analyzing particle trajectories and statistical quantities relevant to particle energization, we find that Fermi-type mechanism dominates the particle acceleration and power-law formation. We find that the results for particle acceleration are scalable as γ 0 and σ 0 increase to large values. The maximum energy for electrons and positrons can reach hundreds of TeV if the wind has a bulk Lorentz factor of γ0 ≈ 106 and a magnetization parameter of σ 0 = 10, which can explain the recent observations of high-energy gamma rays from PWNe

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Improved Models of R Coronae Borealis Stars

We present an improved numerical method to model subsolar He+CO-WD merger progenitors of R Corona Borealis stars that builds on our previous work. These improvements include a smooth entropy transition from the core to the envelope of the post-merger, inclusion of single-zone nucleosynthesis to mimic the effects of burning during the merger event, and post-processing the models with a larger nuclear network for analysis of s-process nucleosynthesis. We perform a parameter study to understand the effects of the entropy transition, peak temperature, and overshooting on our models. The models that best agree with observations of R Corona Borealis stars are processed with a much larger nuclear network to investigate s-process nucleosynthesis and the dredge-up of s-process products into the outer envelope in detail. We present a model with a significant enhancement in s-process elements, which also agrees with observed surface abundances and isotopic ratios of 16 O/ 16 O and C/O between 1 and 10. Finally, we find that the neutron exposure and initial neutron densities this model requires to obtain such an enhancement are much more consistent with i-process nucleosynthesis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Energetic Nonthermal Electrons within the Above-the-looptop Regions in Solar Flares: Acceleration, Feedback, and Quasiperiodic Pulsations

Solar flares are among the most dramatic events in the solar system, releasing substantial magnetic energy and accelerating a large number of electrons to high energies. Notably, in certain events, the above-the-looptop region may contain a significant population of nonthermal electrons, both in number and energy. For the first time, we adopt a novel numerical method that combines magnetohydrodynamics with energetic particles incorporating feedback from nonthermal electrons to investigate electron acceleration and transport in solar flares. We find that a large fraction of energetic electrons are accelerated via the current sheet and termination shock regions. Most energetic electrons are concentrated in the above-the-looptop region, carrying a sizable amount of the released energy. We observe that greater feedback of nonthermal electrons leads to steeper energy spectra. The energy density of the nonthermal electrons oscillates due to the periodic impact of magnetic islands into the above-the-looptop region, which may help explain the observed quasiperiodic pulsations. Our simulations provide new insights into the origin of nonthermal electrons and associated emissions in the above-the-looptop region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

CosmosDG: An hp -adaptive Discontinuous Galerkin Code for Hyper-resolved Relativistic MHD

We have extended Cosmos++, a multidimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and arbitrarily high-order finite element structures. The new finite element implementation, called CosmosDG, is based on a discontinuous Galerkin (DG) formulation, using both entropy-based artificial viscosity and slope limiting procedures for the regularization of shocks. High-order multistage forward Euler and strong-stability preserving Runge–Kutta time integration options complement high-order spatial discretization. We have also added flexibility in the code infrastructure allowing for both adaptive mesh and adaptive basis order refinement to be performed separately or simultaneously in a local (cell-by-cell) manner. We discuss in this report the DG formulation and present tests demonstrating the robustness, accuracy, and convergence of our numerical methods applied to special and general relativistic MHD, although we note that an equivalent capability currently also exists in CosmosDG for Newtonian systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The LSST DESC DC2 Simulated Sky Survey

Here, we describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep area of approximately 300 deg 2 , as well as a deep drilling field of approximately 1 deg 2 . We simulate 5 yr of the planned 10 yr survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the data set to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic test bed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time domain cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. II. General Relativistic Module

Black hole (BH) environments often host plasmas that are fully collisionless or contain intrinsically collisionless regions, including relativistic jets and coronae, where particle energization is ubiquitous. Capturing the physics of these systems requires numerical methods capable of modeling relativistic, magnetized, collisionless plasmas in strong gravitational fields. In this work, we introduce the general relativistic module for Entity—the first open-source, coordinate-agnostic, performance-portable, particle-in-cell code. The code enables fast axisymmetric simulations of collisionless plasmas around BHs on any modern high-performance computing architecture (both GPUs and CPUs).

Galishnikova, Alisa [Flatiron Institute, New York,↗
DOE OSTI2023

Evaluation of a cloudy cold-air pool in the Columbia River basin in different versions of the High-Resolution Rapid Refresh (HRRR) model

The accurate forecast of persistent orographic cold-air pools in numerical weather prediction models is essential for the optimal integration of wind energy into the electrical grid during these events. Model development efforts during the second Wind Forecast Improvement Project (WFIP2) aimed to address the challenges related to this. We evaluated three versions of the National Oceanic and Atmospheric Administration (NOAA) High-Resolution Rapid Refresh model with two different horizontal grid spacings against in situ and remote sensing observations to investigate how developments in physical parameterizations and numerical methods targeted during WFIP2 impacted the simulation of a persistent cold-air pool in the Columbia River basin. Differences amongst model versions were most apparent in simulated temperature and low-level cloud fields during the persistent phase of the cold-air pool. The model developments led to an enhanced low-level cloud cover, resulting in better agreement with the observations. This removed a diurnal cycle in the near-surface temperature bias at stations throughout the basin by reducing a cold bias during the night and a warm bias during the day. However, low-level clouds did not clear sufficiently during daytime in the newest model version, which leaves room for further model developments. The model developments also led to a better representation of the decay of the cold-air pool by slowing down its erosion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Albany/FELIX: A parallel, scalable and robust, finite element, first-order Stokes approximation ice sheet solver built for advanced analysis

This paper describes a new parallel, scalable and robust finite element based solver for the first-order Stokes momentum balance equations for ice flow. The solver, known as Albany/FELIX , is constructed using the component-based approach to building application codes, in which mature, modular libraries developed as a part of the Trilinos project are combined using abstract interfaces and template-based generic programming, resulting in a final code with access to dozens of algorithmic and advanced analysis capabilities. Following an overview of the relevant partial differential equations and boundary conditions, the numerical methods chosen to discretize the ice flow equations are described, along with their implementation. The results of several verification studies of the model accuracy are presented using (1) new test cases for simplified two-dimensional (2-D) versions of the governing equations derived using the method of manufactured solutions, and (2) canonical ice sheet modeling benchmarks. Model accuracy and convergence with respect to mesh resolution are then studied on problems involving a realistic Greenland ice sheet geometry discretized using hexahedral and tetrahedral meshes. Also explored as a part of this study is the effect of vertical mesh resolution on the solution accuracy and solver performance. The robustness and scalability of our solver on these problems is demonstrated. Lastly, we show that good scalability can be achieved by preconditioning the iterative linear solver using a new algebraic multilevel preconditioner, constructed based on the idea of semi-coarsening.

54 ENVIRONMENTAL SCIENCES↗