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

THREE-DIMENSIONAL BOLTZMANN HYDRO CODE FOR CORE COLLAPSE IN MASSIVE STARS. I. SPECIAL RELATIVISTIC TREATMENTS

We propose a novel numerical method for solving multi-dimensional, special relativistic Boltzmann equations for neutrinos coupled with hydrodynamics equations. This method is meant to be applied to simulations of core-collapse supernovae. We handle special relativity in a non-conventional way, taking account of all orders of v/c. Consistent treatment of the advection and collision terms in the Boltzmann equations has been a challenge, which we overcome by employing two different energy grids: Lagrangian remapped and laboratory fixed grids. We conduct a series of basic tests and perform a one-dimensional simulation of core-collapse, bounce, and shock-stall for a 15 M {sub ☉} progenitor model with a minimum but essential set of microphysics. We demonstrate in the latter simulation that our new code is capable of handling all phases in core-collapse supernova. For comparison, a non-relativistic simulation is also conducted with the same code, and we show that they produce qualitatively wrong results in neutrino transfer. Finally, we discuss a possible incorporation of general relativistic effects into our method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

AN OPEN-SOURCE NEUTRINO RADIATION HYDRODYNAMICS CODE FOR CORE-COLLAPSE SUPERNOVAE

We present an open-source update to the spherically symmetric, general-relativistic hydrodynamics, core-collapse supernova (CCSN) code GR1D. The source code is available at http://www.GR1Dcode.org. We extend its capabilities to include a general-relativistic treatment of neutrino transport based on the moment formalisms of Shibata et al. and Cardall et al. We pay special attention to implementing and testing numerical methods and approximations that lessen the computational demand of the transport scheme by removing the need to invert large matrices. This is especially important for the implementation and development of moment-like transport methods in two and three dimensions. A critical component of neutrino transport calculations is the neutrino–matter interaction coefficients that describe the production, absorption, scattering, and annihilation of neutrinos. In this article we also describe our open-source neutrino interaction library NuLib (available at http://www.nulib.org). We believe that an open-source approach to describing these interactions is one of the major steps needed to progress toward robust models of CCSNe and robust predictions of the neutrino signal. We show, via comparisons to full Boltzmann neutrino-transport simulations of CCSNe, that our neutrino transport code performs remarkably well. Furthermore, we show that the methods and approximations we employ to increase efficiency do not decrease the fidelity of our results. We also test the ability of our general-relativistic transport code to model failed CCSNe by evolving a 40-solar-mass progenitor to the onset of collapse to a black hole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

An implementation of the phase-field model based on coupled thermomechanical finite element solvers for large-strain twinning, explicit dynamic fracture and the classical Stefan problem

The implementation of a phase-field model in finite elements usually requires significant expertise and involves the development of a user element with additional degrees of freedom. An alternative implementation of the phase-field model within a thermo-mechanical finite element simulation package was presented in (Cho et al 2012 Int. J. Solids Struct. 49 1973–1992), where the phase-field variable is treated as the temperature degree of freedom. However, this approach has only been used for small strain phase-field modelling of martensitic transformations and quasistatic phase-field modelling of fracture. Here, we present a phase-field finite element implementation via the temperature degree of freedom for several additional cases from the literature: (i) the large-strain phase-field description of deformation twinning presented in (Clayton and Knap 2011 Physica D 240 841–858), (ii) phase-field description of brittle fracture with inertial effects based on the theory from (Molnár and Gravouil 2017, Finite Elem. Anal. Des. 130 27–38) and (Miehe et al 2010 Int. J. Numer. Methods Eng. 83 1273–1311) and (iii) the classical Stefan problem of solidification presented in (Mackenzie and Robertson 2002 J. Comput. Phys. 181 526–544). The last problem involves the temperature and phase-field variables as unknowns.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Linearised Fokker–Planck collision model for gyrokinetic simulations

We introduce a gyrokinetic, linearised Fokker–Planck collision model that satisfies conservation laws and is accurate at arbitrary collisionalities. The differential test-particle component of the operator is exact; the integral field-particle component is approximated using a spherical harmonic and a modified Laguerre polynomial expansion developed by Hirshman and Sigmar (1976 Phys. Fluids 19 1532). The numerical methods of the implementation in the δf-gyrokinetic code stella (Barnes et al 2019 J. Comput. Phys. 391 365–80) are discussed, and conservation properties of the operator are demonstrated. The collision model is then benchmarked against the collision model of the gyrokinetic solver GS2 in the limiting cases of a reduced test-particle collision operator and energy- and momentum-conserving operator. The accuracy of the full collision model is investigated by solving the parallel Spitzer-Härm problem for the transport coefficients. It is shown that retaining collisional energy flux and higher-order terms in the field-particle operator reduces errors in the transport coefficients from 10%–25% for a simple momentum- and energy-conserving model to under 1%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

An overview of discharge plasma modeling for Hall effect thrusters

The current status and challenges of the discharge plasma modeling capabilities in Hall effect thrusters (HETs) are reviewed. Partially magnetized plasmas in HETs present unique physics in that collisions with neutral atoms (a feature of low-temperature plasmas) and instabilities and oscillations (a feature of high-temperature plasmas) coexist. A discharge plasma model is constructed by choosing the numerical methods, e.g. fluid or kinetic, for individual plasma constituents, including neutral atoms, ions, and electrons. The governing equations, boundary conditions, and collision models of a xenon discharge plasma in HETs are discussed. Various models are used to illustrate device-scale phenomena, including the breathing mode and azimuthally rotating spokes, as well as small-scale physics, such as the electron cyclotron drift instability and current-carrying instability that may occur in a hollow cathode plume.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Prediction of strain, inter-layer interaction and critical current in CORC ® wires under axial strain by T-A modeling

Superconducting conductors on round core (CORC ® ) cables and wires can meet the needs of large high-field magnets, such as particle accelerators and compact nuclear fusion machines, due to their simple cabling process, high current-carrying capacity and reliable operation under high mechanical stresses. Many high-field magnets require CORC ® cables to carry a current of thousands of amperes in a background magnetic field exceeding 20 T. As a result, the large electromagnetic forces will deform the cable in the axial direction due to hoop stress and in the transverse direction by compressive stress. Therefore, it is essential to determine the irreversible deformation limit of the CORC ® cable under axial tensile load and optimize the cabling parameters to potentially extend this limit. Analytical and numerical methods are developed to assess the performance degradation of CORC ® wires under axial tensile load. The strain level, interlayer contact pressure and friction and their impact on the critical current are calculated by combining the mechanical response and the T-A method. Analyzing the results shows that the winding angle of the tape and the Poisson’s ratio of the inner core are key factors affecting the irreversible tensile strain limit of CORC ® wires. The smaller the winding angle and the higher the Poisson’s ratio of the inner core, the higher the irreversible tensile strain limit. For multi-layer CORC ® wires, the initial contact pressure caused by the cabling process must also be considered. The inter-layer interaction is coupled with the tape strain of each layer. The results of this research can serve as a basis for optimizing and designing CORC ® wires with extended irreversible strain limits.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Spatial correlations in driven-dissipative photonic lattices

In this paper, we study the nonequilibrium steady-state of interacting photons in cavity arrays as described by the driven-dissipative Bose–Hubbard and spin-1/2 XY model. For this purpose, we develop a self-consistent expansion in the inverse coordination number of the array ($\sim 1/z$) to solve the Lindblad master equation of these systems beyond the mean-field approximation. Our formalism is compared and benchmarked with exact numerical methods for small systems based on an exact diagonalization of the Liouvillian and a recently developed corner-space renormalization technique. We then apply this method to obtain insights beyond mean-field in two particular settings: (i) we show that the gas–liquid transition in the driven-dissipative Bose–Hubbard model is characterized by large density fluctuations and bunched photon statistics. (ii) We study the antibunching–bunching transition of the nearest-neighbor correlator in the driven-dissipative spin-1/2 XY model and provide a simple explanation of this phenomenon.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Cosmological perturbations during the Bose-Einstein condensation of dark matter

In the present work, we analyze the evolution of the scalar and tensorial perturbations and the quantities relevant for the physical description of the Universe, as the density contrast of the scalar perturbations and the gravitational waves energy density during the Bose-Einstein condensation of dark matter. The behavior of these parameters during the Bose-Einstein phase transition of dark matter is analyzed in details. To study the cosmological dynamics and evolution of scalar and tensorial perturbations in a Universe with and without cosmological constant we use both analytical and numerical methods. The Bose-Einstein phase transition modifies the evolution of gravitational waves of cosmological origin, as well as the process of large-scale structure formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Stability of BEC galactic dark matter halos

In this paper we show that spherically symmetric BEC dark matter halos, with the sin r/r density profile, that accurately fit galactic rotation curves and represent a potential solution to the cusp-core problem are unstable. We do this by introducing back the density profiles into the fully time-dependent Gross-Pitaevskii-Poisson system of equations. Using numerical methods to track the evolution of the system, we found that these galactic halos lose mass at an approximate rate of half of its mass in a time scale of dozens of Myr. We consider this time scale is enough as to consider these halos are unstable and unlikely to be formed. We provide some arguments to show that this behavior is general and discuss some other drawbacks of the model that restrict its viability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Inflation in a two 3-form fields scenario

A setting constituted by N 3-form fields, without any direct interaction between them, minimally coupled to gravity, is introduced in this paper as a framework to study the early evolution of the universe. We focus particularly on the two 3-forms case. An inflationary scenario is found, emerging from the coupling to gravity. More concretely, the fields coupled in this manner exhibit a complex interaction, mediated by the time derivative of the Hubble parameter. Our investigation is supported by means of a suitable choice of potentials, employing numerical methods and analytical approximations. In more detail, the oscillations on the small field limit become correlated, and one field is intertwined with the other. In this type of solution, a varying sound speed is present, together with the generation of isocurvature perturbations. The mentioned features allow to consider an interesting model, to test against observation. It is subsequently shown how our results are consistent with current CMB data (viz.Planck and BICEP2)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Caustics of 1/r{sup n} binary gravitational lenses: from galactic haloes to exotic matter

We investigate the caustic topologies for binary gravitational lenses made up of two objects whose gravitational potential declines as 1/r{sup n}. With n<1 this corresponds to power-law dust distributions like the singular isothermal sphere. The n>1 regime can be obtained with some violations of the energy conditions, one famous example being the Ellis wormhole. Gravitational lensing provides a natural arena to distinguish and identify such exotic objects in our Universe. We find that there are still three topologies for caustics as in the standard Schwarzschild binary lens, with the main novelty coming from the secondary caustics of the close topology, which become huge at higher n. After drawing caustics by numerical methods, we derive a large amount of analytical formulae in all limits that are useful to provide deeper insight in the mathematics of the problem. Our study is useful to better understand the phenomenology of galaxy lensing in clusters as well as the distinct signatures of exotic matter in complex systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A perturbative approach to the redshift space power spectrum: beyond the Standard Model

We develop a code to produce the power spectrum in redshift space based on standard perturbation theory (SPT) at 1-loop order. The code can be applied to a wide range of modified gravity and dark energy models using a recently proposed numerical method by A.Taruya to find the SPT kernels. This includes Horndeski's theory with a general potential, which accommodates both chameleon and Vainshtein screening mechanisms and provides a non-linear extension of the effective theory of dark energy up to the third order. Focus is on a recent non-linear model of the redshift space power spectrum which has been shown to model the anisotropy very well at relevant scales for the SPT framework, as well as capturing relevant non-linear effects typical of modified gravity theories. We provide consistency checks of the code against established results and elucidate its application within the light of upcoming high precision RSD data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Anisotropic deformations of spatially open cosmology in massive gravity theory

We combine analytical and numerical methods to study anisotropic deformations of the spatially open homogeneous and isotropic cosmology in the ghost free massive gravity theory with flat reference metric. We find that if the initial perturbations are not too strong then the physical metric relaxes back to the isotropic de Sitter state. However, the dumping of the anisotropies is achieved at the expense of exciting the Stueckelberg fields in such a way that the reference metric changes and does not share anymore with the physical metric the same rotational and translational symmetries. As a result, the universe evolves towards a fixed point which does not coincide with the original solution, but for which the physical metric is still de Sitter. If the initial perturbation is strong, then its evolution generically leads to a singular anisotropic state or, for some parameter values, to a decay into flat spacetime. We also present an infinite dimensional family of new homogeneous and isotropic cosmologies in the theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Scalar geons in Born-Infeld gravity

The existence of static, spherically symmetric, self-gravitating scalar field solutions in the context of Born-Infeld gravity is explored. Upon a combination of analytical approximations and numerical methods, the equations for a free scalar field (without a potential term) are solved, verifying that the solutions recover the predictions of General Relativity far from the center but finding important new effects in the central regions. We find two classes of objects depending on the ratio between the Schwarzschild radius and a length scale associated to the Born-Infeld theory: massive solutions have a wormhole structure, with their throat at r ≈ 2 M , while for the lighter configurations the topology is Euclidean. The total energy density of these solutions exhibits a solitonic profile with a maximum peaked away from the center, and located at the throat whenever a wormhole exists. The geodesic structure and curvature invariants are analyzed for the various configurations considered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

A digital CDS technique and its performance testing

Readout noise is a critical parameter for characterizing the performance of charge-coupled devices (CCDs), which can be greatly reduced by the correlated double sampling (CDS) circuit. However, a conventional CDS circuit inevitably introduces new noise since it consists of several active analog components such as operational amplifiers. This article proposes a digital CDS circuit technique, which transforms the pre-amplified CCD signal into a train of digital presentations by a high-speed data acquisition card directly without the noisy CDS circuit, then implements the digital CDS algorithm through a numerical method. A readout noise of 3.3 e - and an energy resolution of 121 eV@5.9 keV can be achieved via the digital CDS technique.

42 ENGINEERING↗
DOE OSTI2020

Construction of quasisymmetric stellarators using a direct coordinate approach

In this work, optimized stellarator configurations and their analytical properties are obtained using a near-axis expansion approach. Such configurations are associated with good confinement as the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak. This makes them appealing as fusion reactor candidates. Using a direct coordinate approach, where the magnetic field and flux surface functions are found explicitly in terms of the position vector at successive orders in the distance to the axis, the set of ordinary differential equations for first and second order quasisymmetry is derived. Examples of quasi-axisymmetric shapes are constructed using a pseudospectral numerical method. Finally, the direct coordinate approach is used to independently verify two hypotheses commonly associated with quasisymmetric magnetic fields, namely that the number of equations exceeds the number of parameters at third order in the expansion and that the near-axis expansion does not prohibit exact quasisymmetry from being achieved on a single flux surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Bifurcation-like transition of divertor conditions induced by X-point radiation in KSTAR L-mode plasmas *

Abstract Density ramps with ion grad B drift directed into lower single null KSTAR L-mode plasmas are associated with a simultaneous and abrupt reduction of the divertor particle flux on both low- and high-field-side targets when the mid-plane line averaged electron density reaches a given level. Target embedded Langmuir probe signals show a clear ‘cliff edge’ behavior similar to that observed in the divertor target electron temperature in DIII-D H-mode plasmas (Eldon et al 2017 Nucl. Fusion 57 066039; McLean et al 2015 J. Nucl. Mater. 463 533–6). The collapse of the particle flux is observed along the whole divertor target area (from private flux region to the far scrape-off layer (SOL)). The critical upstream density of this target flux cliff is invariant under fuel gas throughput modulation. The transition along the cliff occurs in tens of milliseconds. With the cliff, carbon impurities and deuterium neutrals transported through the X-point to the core produce a strong radiation spot near the X-point, seen on bolometric signals, and increase the upstream density. The experimental observations are consistent with time-dependent SOLPS-ITER simulations, which also demonstrate an abrupt transition of the target flux and upstream density with the increase in X-point radiation. The timescale of the cliff predicted by SOLPS-ITER is consistent with the experiment, although, it is influenced by gas throughput or time-dependent numerical methods. In the L-mode phase space of separatrix electron density and temperature, branches are divided based on target temperature, because the latter is strongly coupled to the radiation front and ionization front due to the monotonic characteristic of the parallel electron temperature distribution. Since the H-mode condition operates at a much higher upstream density and electron temperature in phase space, dissipation from sputtered carbon alone leads to the density limit before reaching the X-point radiation condition. This is therefore consistent with the fact that cliffs have never been observed in H-mode KSTAR experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗