Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical method”

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

Elevated temperature mechanical properties of Inconel 617 surface oxide using nanoindentation

Inconel 617 is a principal candidate material for helium gas cooled very-high-temperature reactors with outlet temperatures of 700–950 °C. Recent findings showed that this alloy develops unique surface oxide layers especially at high temperature (HT) helium environment with distinctive wear, friction and contact properties. This study investigates the elevated temperature mechanical properties of Inconel 617 top surface layers aged in HT helium environment. Nanoindentation technique is used to obtain load-displacement graphs of the alloy top surface oxide in temperatures ranging from 25 °C up to 600 °C. In addition, using finite element analysis along with an iterative regression technique, a semi-numerical method is developed to further measure and quantify the material parameters and, in particular, time-independent and creep characteristics of the oxide. Although Young's modulus of the oxide is found to be relatively close to the bulk for the tested temperatures, the yield strength and hardness, in comparison to the bulk material, increase significantly as the material is oxidized after aging. The oxide exhibits significant softening as the temperature increases to 600 °C. Unlike the bulk material, diffusion through the grains is found to be the dominant creep mechanism for the oxide. Considerable difference between the mechanical properties of the oxide and the bulk material shows the need for accurate measurements of near surface mechanical properties, if reliable predictive contact and tribological models are sought at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Precise Lifetime Measurements in Light Nuclei for Benchmarking Modern Ab-initio Nuclear Structure Models

A new generation of ab-initio calculations, based on realistic two- and three-body forces, is having a profound impact on our view of how nuclei work. To improve the numerical methods, and the parameterization of 3-body forces, new precise data are needed. Electromagnetic transitions are very sensitive to the dynamics which drive mixing between configurations. We have made a series of precise (< 3%) measurements of electromagnetic transitions in the A=10 nuclei {sup 10}C and {sup 10}Be by using the Doppler Shift Attenuation method carefully. Many interesting features can be reproduced including the strong α clustering. New measurements on {sup 8}Be and {sup 12}Be highlight the interplay between the alpha clusters and their valence neutrons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Optimization of the magnetic horn for the nuSTORM non-conventional neutrino beam using the genetic algorithm

This paper describes the strategy for optimizing the magnetic horn for the neutrinos from STORed Muons (nuSTORM) facility. The nuSTORM magnetic horn is the primary collection device for the secondary particles generated by bombarding a solid target with 120 GeV protons. As a consequence of the non-conventional beamline designed for nuSTORM, the requirements on the horn are different from those for a conventional neutrino beamline. At nuSTORM, muons decay while circulating in the storage ring, and the detectors are placed downstream of the production straight so as to be exposed to the neutrinos from muon decay. nuSTORM aims at precisely measuring the neutrino cross sections, and providing a definitive statement about the existence of sterile neutrinos. The nuSTORM horn aims at focusing the pions into a certain phase space so that more muons from pion decay can be accepted by the decay ring. The paper demonstrates a numerical method that was developed to optimize the horn design to gain higher neutrino flux from the circulating muons. A Genetic Algorithm (GA) was applied to the simultaneous optimization of the two objectives in this study. In conclusion, the application of the technique discussed in this paper is not limited to either the nuSTORM facility or muon based facilities, but can be used for other neutrino facilities that use magnetic horns as collection devices.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Numerical study and test of the APS linac transverse deflecting cavity

The transverse deflecting cavity can be used to transform particle distributions in the 6D phase space, which makes it a promising component in phase space beam diagnostics and beam manipulations. In the Advanced Photon Source (APS) Linac, a LOLA-type traveling wave deflecting cavity was installed for the diagnostics of beam characteristics, such as the bunch length, bunch temporal profile, time-dependent energy spread and slice (time-correlated) transverse emittance. In this paper, this deflecting cavity is modeled and analyzed with the numerical method. The effects of the center cell and coupler cell dimensions on the performance of the whole structure are studied, which shows the coupler cell radius has a dominant effect over the coupler cell length and slot width on the global reflection coefficient and field flatness. Important RF parameters, such as S-11, field flatness and phase advance etc., of the LOLA-II cavity in tuning are calculated and discussed. After optimization, the field flatness of the cavity is 2.1%, the phase advance is 119.66 degrees degrees with a standard deviation smaller than 0.5 degrees and the bandwidth is 15.738 MHz when VSWR < 1.1 (i.e. the S-11, < -26 dB). Bead-pull measurements and RF conditioning have been also performed and discussed. Finally, the initial cavity commissioning with electron beam of 375 MeV is reported, and compared with simulation.

coupler↗
DOE OSTI2017

Simulations of single- and two-phase shock tubes and gravity-driven wave problems with the RELAP-7 nuclear reactor system analysis code

Propagation of pressure waves in single and two-phase flows is of key interest for nuclear safety analyses, reactor design, construction and life extension of nuclear reactors. In this work, we assess the capabilities of RELAP-7, the next-generation nuclear reactor system code, to solve single- and two-phase flows developing pressure waves. RELAP-7 is based on a 7-equation, two-phase flow model, with distinct phasic pressures. It is well established that these hyperbolic conservation laws can develop shocks and discontinuities and thus, require a stabilization numerical method. The all-Mach flow Entropy Viscosity Method (EVM) is now employed in RELAP-7 as a viscous stabilization approach. In this paper, we briefly describe the implementation of the EVM and use open-literature test cases that exhibit pressure waves to show that the entropy viscosity method performs adequately for these single- and two-phase cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

CTF-R: A novel residual-based thermal hydraulic solver

The traditional scientific process has been revolutionized by the advent of computational modeling, but the nuclear industry generally uses "legacy codes," which were developed early in the evolution of computers. One example of a legacy code, the thermal hydraulic subchannel code CTF, is modernized in this work through the development of a novel residual-based version, CTF-R. Unlike its predecessor, CTF-R is not limited by the strict computational limitations of the early 1980's, and can therefore be designed such that it is inherently flexible and easy to use. A case study is examined to demonstrate how the flexibility of the code can be used to improve simulation results. In this example, the coupling between the solid and liquid fields is examined. Traditionally, this coupling is modeled explicitly, which imposes numerical stability limits on the time step size. These limits are derived and it is shown that they are removed when the coupling is made implicit. Further, the development of CTF-R will enable future improvements in next generation reactor modeling, numerical methods, and coupling to other codes. Through the further development of CTF-R and other residual-based codes, state-of-the-art simulation is possible.

42 ENGINEERING↗
DOE OSTI2011

Gyrokinetic particle-in-cell optimization on emerging multi- and manycore platforms

The next decade of high-performance computing (HPC) systems will see a rapid evolution and divergence of multi- and manycore architectures as power and cooling constraints limit increases in microprocessor clock speeds. Understanding efficient optimization methodologies on diverse multicore designs in the context of demanding numerical methods is one of the greatest challenges faced today by the HPC community. In this paper, we examine the efficient multicore optimization of GTC, a petascale gyrokinetic toroidal fusion code for studying plasma microturbulence in tokamak devices. For GTC’s key computational components (charge deposition and particle push), we explore efficient parallelization strategies across a broad range of emerging multicore designs, including the recently-released Intel Nehalem-EX, the AMD Opteron Istanbul, and the highly multithreaded Sun UltraSparc T2+. We also present the first study on tuning gyrokinetic particle-in-cell (PIC) algorithms for graphics processors, using the NVIDIA C2050 (Fermi). Our work discusses several novel optimization approaches for gyrokinetic PIC, including mixed-precision computation, particle binning and decomposition strategies, grid replication, SIMDized atomic floating-point operations, and effective GPU texture memory utilization. Overall, we achieve significant performance improvements of 1.3–4.7× on these complex PIC kernels, despite the inherent challenges of data dependency and locality. Finally, our work also points to several architectural and programming features that could significantly enhance PIC performance and productivity on next-generation architectures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

The physics of irradiation of biological matter by ion beams

In radiobiology, predicting the evolution of irradiated biological matter is nowadays an active field of research to identify DNA lesions or to adapt the radiotherapeutic protocols in radiation oncology. In this context, the numerical methods, based on Monte Carlo track-structure simulations, represent the most suitable and powerful tools for understanding the radiobiological damages induced by ionizing particles. Here, we report the theoretical differential and total cross sections, computed within the quantum mechanical continuum distorted wave-eikonal initial state (CDW-EIS) method, for ion impact on water vapor and DNA nucleobases. These cross sections have been used to build up the input database for the homemade Monte Carlo track-structure TILDA-V. A comparison between the theoretical predictions and the available experimental data is presented. Micro-dosimetry results obtained with TILDA-V are also reported.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Unsteady aerodynamics of large-scale floating offshore wind turbines in surge motion

Unsteady aerodynamic loads significantly influence the design and wake flow field of floating offshore wind turbines, especially due to wave- or vibration-induced tower top motions triggering various unsteady phenomena. Recent studies show that increasing turbine sizes amplify unsteady aerodynamic effects, as their impact typically grows with rotor diameter. This work combines recent findings from experiments and simulations on model-scale FOWT aerodynamics with new numerical analyses of large-scale rotors, providing a comprehensive understanding of unsteady phenomena occurrence and impact. Numerical analyses of the IEA 15-MW and 22-MW rotors undergoing surge motions characterise the combined influence of motion-induced unsteady phenomena on rotor thrust. Results indicate that unsteady effects can reduce thrust force variations by up to 40% at realistic surge periods. These findings contrast with prior model-scale rotor investigations, attributed to the specialised design of the model rotors. Comparisons between numerical methods - the dynamic blade element momentum method in OpenFAST and free vortex wake modules in panMARE and OpenFAST - reveal persistent differences in thrust predictions under both idealised and realistic conditions. This highlights that classical blade element momentum approaches require enhancement to accurately capture unsteady loads on large-scale floating offshore wind turbines at low wind speeds.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Groundwater flowpath characteristics drive variability in per- and polyfluoroalkyl substances (PFAS) loading across a stream-wetland system

Groundwater dependent ecosystems in areas with industrial and military land use are at risk of direct exposure to a wide range of contaminants, including PFAS chemicals. Glaciated terrain often has mixed high and low permeability sediments coupled with groundwater flow-through lake features. These hydrogeologic attributes create highly complex ‘source to seep’ dynamics that make spatiotemporal contaminant transport patterns difficult to predict. We investigated one such system in detail using a suite of heat-tracing and chemical methods. Numerous (n=57) preferential groundwater discharge zones (vertical flux rates ranging 0.2 to 3.2 m/d) were identified across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, adjacent to an Air Force Base with several known PFAS source areas. Surface-water and groundwater samples were collected and analyzed (n=145) for precursors and terminal PFAS compounds between March and September 2022. Samples were collected at identified seeps along the Quashnet River (n=59), from wells upgradient from the stream-wetland system (n= 44), from contributing flow-through kettle lakes (n=8), and at multiple locations along the Quashnet River (n=34). Samples from seeps and wells had measured PFAS concentrations ranging from non-detect to approximate 3,500 ng/L (mean= 1,650 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to emergence at the discharge zones. Groundwater-lake interaction along flowpaths that sourced the sampled seeps was farther supported by significant correlations (p < 0.01) between deuterium excess and %PFAS precursors, and between %PFAS precursors and multiple terminal PFAS compounds (e.g., PFPeS, PFBS, PFHxS). However, some sampled seeps contributing groundwater to the stream-wetland system had much higher total PFAS concentrations (>1000 ng/L) than the upgradient kettle lakes, despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. PFAS compounds and water isotopic composition at sampled multilevel groundwater wells, rivers, lakes, and seeps suggest that a complex mixture of source groundwater and flowpath characteristics are responsible for diverse observed PFAS mixtures at preferential discharge zones across the stream-wetland system. Further, total PFAS loading patterns to the Quashnet River via groundwater discharge remained remarkably similar from winter to summer to fall conditions, despite a regional dry period in late summer 2022 with the upper river channel completely drying. This work addresses gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations in controlling contaminant mass loading in preferential groundwater discharge zones and presents a transferrable field toolkit for efficient characterization of spatially preferential PFAS transport dynamics.

Contaminant transport↗
DOE OSTI2010

Heat extraction from salinity-gradient solar ponds using heat pipe heat exchangers

This paper presents the results of experimental and theoretical analysis on the heat extraction process from solar pond by using the heat pipe heat exchanger. In order to conduct research work, a small scale experimental solar pond with an area of 7.0 m{sup 2} and a depth of 1.5 m was built at Khon Kaen in North-Eastern Thailand (16 27'N102 E). Heat was successfully extracted from the lower convective zone (LCZ) of the solar pond by using a heat pipe heat exchanger made from 60 copper tubes with 21 mm inside diameter and 22 mm outside diameter. The length of the evaporator and condenser section was 800 mm and 200 mm respectively. R134a was used as the heat transfer fluid in the experiment. The theoretical model was formulated for the solar pond heat extraction on the basis of the energy conservation equations and by using the solar radiation data for the above location. Numerical methods were used to solve the modeling equations. In the analysis, the performance of heat exchanger is investigated by varying the velocity of inlet air used to extract heat from the condenser end of the heat pipe heat exchanger (HPHE). Air velocity was found to have a significant influence on the effectiveness of heat pipe heat exchanger. In the present investigation, there was an increase in effectiveness by 43% as the air velocity was decreased from 5 m/s to 1 m/s. The results obtained from the theoretical model showed good agreement with the experimental data. (author)

14 SOLAR ENERGY↗
DOE OSTI2011

Effects of alkyl polyglycoside (APG) on composting of agricultural wastes

Composting is the biological degradation and transformation of organic materials under controlled conditions to promote aerobic decomposition. To find effective ways to accelerate composting and improve compost quality, numerous methods including additive addition, inoculation of microorganisms, and the use of biosurfactants have been explored. Studies have shown that biosurfactant addition provides more favorable conditions for microorganism growth, thereby accelerating the composting process. However, biosurfactants have limited applications because they are expensive and their use in composting and microbial fertilizers is prohibited. Meanwhile, alkyl polyglycoside (APG) is considered a 'green' surfactant. This study aims to determine whether APG addition into a compost reaction vessel during 28-day composting can enhance the organic matter degradation and composting process of dairy manure. Samples were periodically taken from different reactor depths at 0, 3, 5, 7, 14, 21, and 28 days. pH levels, electrical conductivity (EC), ammonium and nitrate nitrogen, seed germination indices, and microbial population were determined. Organic matter and total nitrogen were also measured. Compared with the untreated control, the sample with APG exhibited slightly increased microbial populations, such as bacteria, fungi, and actinomycetes. APG addition increased temperatures without substantially affecting compost pH and EC throughout the process. After 28 days, APG addition increased nitrate nitrogen concentrations, promoted matter degradation, and increased seed germination indices. The results of this study suggest that the addition of APG provides more favorable conditions for microorganism growth, slightly enhancing organic matter decomposition and accelerating the composting process, improving the compost quality to a certain extent.

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

Explosive X-point collapse in relativistic magnetically dominated plasma

The extreme properties of the gamma-ray flares in the Crab nebula present a clear challenge to our ideas on the nature of particle acceleration in relativistic astrophysical plasma. It seems highly unlikely that standard mechanisms of stochastic type are at work here and hence the attention of theorists has switched to linear acceleration in magnetic reconnection events. In this series of papers, we attempt to develop a theory of explosive magnetic reconnection in highly magnetized relativistic plasma which can explain the extreme parameters of the Crab flares. In the first paper, we focus on the properties of the X-point collapse. Using analytical and numerical methods (fluid and particle-in-cell simulations) we extend Syrovatsky’s classical model of such collapse to the relativistic regime. We find that the collapse can lead to the reconnection rate approaching the speed of light on macroscopic scales. During the collapse, the plasma particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\unicode[STIX]{x1D70E}$ . For sufficiently high magnetizations and vanishing guide field, the non-thermal particle spectrum consists of two components: a low-energy population with soft spectrum that dominates the number census; and a high-energy population with hard spectrum that possesses all the properties needed to explain the Crab flares.

Physics↗
DOE OSTI2017

Particle acceleration in relativistic magnetic flux-merging events

Using analytical and numerical methods (fluid and particle-in-cell simulations) we study a number of model problems involving merger of magnetic flux tubes in relativistic magnetically dominated plasma. Mergers of current-carrying flux tubes (exemplified by the two-dimensional ‘ABC’ structures) and zero-total-current magnetic flux tubes are considered. In all cases regimes of spontaneous and driven evolution are investigated. We identify two stages of particle acceleration during flux mergers: (i) fast explosive prompt X-point collapse and (ii) ensuing island merger. The fastest acceleration occurs during the initial catastrophic X-point collapse, with the reconnection electric field of the order of the magnetic field. During the X-point collapse, particles are accelerated by charge-starved electric fields, which can reach (and even exceed) values of the local magnetic field. The explosive stage of reconnection produces non-thermal power-law tails with slopes that depend on the average magnetization $\unicode[STIX]{x1D70E}$ . For plasma magnetization $\unicode[STIX]{x1D70E}\leqslant 10^{2}$ the spectrum power-law index is $p>2$ ; in this case the maximal energy depends linearly on the size of the reconnecting islands. For higher magnetization, $\unicode[STIX]{x1D70E}\geqslant 10^{2}$ , the spectra are hard, $p<2$ , yet the maximal energy $\unicode[STIX]{x1D6FE}_{\text{max}}$ can still exceed the average magnetic energy per particle, ${\sim}\unicode[STIX]{x1D70E}$ , by orders of magnitude (if $p$ is not too close to unity). The X-point collapse stage is followed by magnetic island merger that dissipates a large fraction of the initial magnetic energy in a regime of forced magnetic reconnection, further accelerating the particles, but proceeds at a slower reconnection rate.

Physics↗
DOE OSTI2023

The DESC stellarator code suite Part 3: Quasi-symmetry optimization

The DESC stellarator optimization code takes advantage of advanced numerical methods to search the full parameter space much faster than conventional tools. Only a single equilibrium solution is needed at each optimization step thanks to automatic differentiation, which efficiently provides exact derivative information. A Gauss–Newton trust-region optimization method uses second-order derivative information to take large steps in parameter space and converges rapidly. With just-in-time compilation and GPU portability, high-dimensional stellarator optimization runs take orders of magnitude less computation time with DESC compared to other approaches. This paper presents the theory of the DESC fixed-boundary local optimization algorithm along with demonstrations of how to easily implement it in the code. Example quasi-symmetry optimizations are shown and compared to results from conventional tools. Three different forms of quasi-symmetry objectives are available in DESC, and their relative advantages are discussed in detail. In the examples presented, the triple product formulation yields the best optimization results in terms of minimized computation time and particle transport. This paper concludes with an explanation of how the modular code suite can be extended to accommodate other types of optimization problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗
DOE OSTI2024

Stellarator optimization with constraints

In this work we consider the problem of optimizing a stellarator subject to hard constraints on the design variables and physics properties of the equilibrium. We survey current numerical methods for handling these constraints, and summarize a number of methods from the wider optimization community that have not been used extensively for stellarator optimization thus far. We demonstrate the utility of new methods of constrained optimization by optimizing a quasi-axisymmetric stellarator for favourable physics properties while preventing strong shaping of the plasma boundary, which can be difficult to create with external current sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗