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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Gravitational Self-force Errors of Poisson Solvers on Adaptively Refined Meshes

An error in the gravitational force that the source of gravity induces on itself (a self-force error) violates both the conservation of linear momentum and the conservation of energy. If such errors are present in a self-gravitating system and are not sufficiently random to average out, the obtained numerical solution will become progressively more unphysical with time: the system will acquire or lose momentum and energy due to numerical effects. In this paper, we demonstrate how self-force errors can arise in the case where self-gravity is solved on an adaptively refined mesh when the refinement is nonuniform. Here, we provide the analytical expression for the self-force error and numerical examples that demonstrate such self-force errors in idealized settings. We also show how these errors can be corrected to an arbitrary order by straightforward addition of correction terms at the refinement boundaries.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Dynamics of a relativistic jet through magnetized media

The merger of two neutron stars (NSs) produces the emission of gravitational waves, the formation of a compact object surrounded by a dense and magnetized environment. If the binary undergoes delayed collapse a collimated and relativistic jet, which will eventually produce a short gamma-ray burst (SGRB), may be launched. The interaction of the jet with the environment has been shown to play a major role in shaping the structure of the outflow that eventually powers the gamma-ray emission. In this paper, we present a set of 2.5D RMHD simulations that follow the evolution of a relativistic non-magnetized jet through a medium with different magnetization levels, as produced after the merger of two NSs. We find that the predominant consequence of a magnetized ambient medium is that of suppressing instabilities within the jet and preventing the formation of a series of collimation shocks. One implication of this is that internal shocks lose efficiency, causing bursts with low-luminosity prompt emission. On the other hand, the jet-head velocity and the induced magnetization within the jet are fairly independent of the magnetization of the ambient medium. Future numerical studies with a larger domain are necessary to obtain light curves and spectra in order to better understand the role of magnetized media.

(magnetohydrodynamics) MHD↗
DOE OSTI2022

MAPPRAISER: A massively parallel map-making framework for multi-kilo pixel CMB experiments

Forthcoming cosmic microwave background (CMB) polarized anisotropy experiments have the potential to revolutionize our understanding of the Universe and fundamental physics. The sought-after, tale-telling signatures will be however distributed over voluminous data sets which these experiments will collect. These data sets will need to be efficiently processed and unwanted contributions due to astrophysical, environmental, and instrumental effects characterized and efficiently mitigated in order to uncover the signatures. This poses a significant challenge to data analysis methods, techniques, and software tools which will not only have to be able to cope with huge volumes of data but to do so with unprecedented precision driven by the demanding science goals posed for the new experiments. A keystone of efficient CMB data analysis is solvers of very large linear systems of equations. Such systems appear in very diverse contexts throughout CMB data analysis pipelines, however they typically display similar algebraic structures and can therefore be solved using similar numerical techniques. Linear systems arising in the so-called map-making problem are one of the most prominent and common ones. In this work we present a massively parallel, flexible and extensible framework, comprised of a numerical library, MIDAPACK, and a high level code, MAPPRAISER, which provide tools for solving efficiently such systems. Here, the framework implements iterative solvers based on conjugate gradient techniques: enlarged and preconditioned using different preconditioners. We demonstrate the framework on simulated examples reflecting basic characteristics of the forthcoming data sets issued by ground-based and satellite-borne instruments, executing it on as many as 16,384 compute cores. The software is developed as an open source project freely available to the community at: https://github.com/B3Dcmb/midapack.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Study of the Z{sub 3} symmetry in QCD at finite temperature and chemical potential using a worm algorithm

Traditional Monte Carlo simulations of QCD in the presence of a baryon chemical potential are plagued by the complex phase problem and new numerical approaches are necessary for studying the phase diagram of the theory. In this work we consider a Z{sub 3} Polyakov loop model for the deconfining phase transition in QCD and discuss how a flux representation of the model in terms of dimer and monomer variable solves the complex action problem. We present results of numerical simulations using a worm algorithm for the specific heat and two-point correlation function of Polyakov loops. Evidences of a first order deconfinement phase transition are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A new re-redistribution scheme for weighted state redistribution with adaptive mesh refinement

State redistribution (SRD) is a recently developed technique for stabilizing cut cells that result from finite-volume embedded boundary methods. SRD has been successfully applied to a variety of compressible and incompressible flow problems. When used in conjunction with adaptive mesh refinement (AMR), additional steps are needed to preserve the accuracy and conservation properties of the solution if the embedded boundary is not restricted to a single level of the mesh hierarchy. In this work, we extend the weighted state redistribution algorithm to cases where cut cells live at or near a coarse-fine interface within the domain. Here, we present numerical results that demonstrate that the algorithm is conservative when the coarse-fine interface intersects the embedded boundary. Additionally we compare the numerical solution of the Sod shock tube problem in an inclined cylinder with the analytic solution, and we compare the simulation of a shock hitting a cylindrical obstacle with experimental data. Finally we demonstrate the methodology for simulation of the multicomponent compressible Navier-Stokes equations in a piston-bowl geometry, and discuss the computational efficiency gained by not requiring the entire embedded boundary to be defined at the finest level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Robust inference of ecosystem soil water stress from eddy covariance data

Eddy covariance data are invaluable for determining ecosystem water use strategies under soil water stress. However, existing stress inference methods require numerous subjective data processing and model specification assumptions whose effect on the inferred soil water stress signal is rarely quantified. These uncertainties may confound the stress inference and the generalization of ecosystem water use strategies across multiple sites and studies. In this research, we quantify the sensitivity of soil water stress signals inferred from eddy covariance data to the prevailing data and modeling assumptions (i.e., their robustness) to compile a comprehensive list of sites with robust soil water stress signals and assess the performance of current stress inference methods. To accomplish this, we identify the most prevalent assumptions from the literature and perform a digital factorial experiment to extract probability distributions of plausible soil water stress signals and model performance at 151 FLUXNET2015 and AmeriFlux-FLUXNET sites. Here, we develop a new framework that summarizes these probability distributions to classify and rank the robustness of each site’s soil water stress signal, which we display with a user-friendly heat map. We estimate that only 5%–36% of sites exhibit a robust soil water stress signal due to deficient model performance and poorly constrained ecosystem water use parameters. We also find that the lack of robustness is site-specific, which undermines grouping stress signals by broad ecosystem categories or comparing results across studies with differing assumptions. Lastly, existing stress inference methods appear better suited for eddy covariance sites with grass/annual vegetation. Our findings call for more careful and consistent inference of ecosystem water stress from eddy covariance data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multi-fidelity microstructure-induced uncertainty quantification by advanced Monte Carlo methods

Quantifying uncertainty associated with the microstructure variation of a material can be a computationally daunting task, especially when dealing with advanced constitutive models and fine mesh resolutions in the crystal plasticity finite element method (CPFEM). Numerous studies have been conducted regarding the sensitivity of material properties and performance to the mesh resolution and choice of constitutive model. However, a unified approach that accounts for various fidelity parameters, such as mesh resolutions, integration time-steps and constitutive models simultaneously is currently lacking. Further, this paper proposes a novel uncertainty quantification (UQ) approach for computing the properties and performance of homogenized materials using CPFEM, that exploits a hierarchy of approximations with different levels of fidelity. In particular, we illustrate how multi-level sampling methods, such as multi-level Monte Carlo (MLMC) and multi-index Monte Carlo (MIMC), can be applied to assess the impact of variations in the microstructure of polycrystalline materials on the predictions of homogenized materials properties. We show that by adaptively exploiting the fidelity hierarchy, we can significantly reduce the number of microstructures required to reach a certain prescribed accuracy. Finally, we show how our approach can be extended to a multi-fidelity framework, where we allow the underlying constitutive model to be chosen from either a phenomenological plasticity model or a dislocation-density-based model.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Numerical self-consistent field theory study of the response of strong polyelectrolyte brushes to external electric fields

The response of strong polyelectrolyte (PE) brushes grafted on an electrode to electric fields generated by opposite surface charges on the PE-grafted electrode and a second parallel electrode has been numerically investigated by self-consistent field theory. The influences of grafting density, average charge fraction, salt concentration, and mobile ion size on the variation of the brush height against an applied voltage bias were investigated. In agreement with molecular dynamics simulation results, a higher grafting density requires a larger magnitude of voltage bias to achieve the same amount of relative change in the brush height. In the experimentally relevant parameter regime of the applied voltage, the brush height becomes insensitive to the voltage bias when the grafting density is high. Including the contribution of surface charges on the grafting electrode, overall charge neutrality inside the PE brushes is generally maintained, especially for PE brushes with high grafting density and high average charge fraction. Our numerical study further reveals that the electric field across the two electrodes is highly non-uniform because of the complex interplay between the surface charges on the electrodes, the charges on the grafted PE chains, and counterions.

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

Chaotic thermalization in classical gauge theories

We explore the idea that chaos concepts might be useful for understanding the thermalization in gauge theories. The SU(2) Higgs model is discussed as a prototype of system with gauge fields coupled to matter fields. Through the numerical solution of the equations of motion, we are able to characterize chaotic behavior via the corresponding Lyapunov exponent. Then it is demonstrated that the system's approach to equilibrium can be understood through direct application of the principles of Statistical Mechanics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Numerical investigation of bubble nonlinear dynamics characteristics

The complicated dynamical behaviors of bubble oscillation driven by acoustic wave can provide favorable conditions for many engineering applications. On the basis of Keller-Miksis model, the influences of control parameters, including acoustic frequency, acoustic pressure and radius of gas bubble, are discussed by utilizing various numerical analysis methods, Furthermore, the law of power spectral variation is studied. It is shown that the complicated dynamic behaviors of bubble oscillation driven by acoustic wave, such as bifurcation and chaos, further the stimulated scattering processes are revealed.

42 ENGINEERING↗
DOE OSTI2014

Development of Variational Guiding Center Algorithms for Parallel Calculations in Experimental Magnetic Equilibria

Structure-preserving algorithms obtained via discrete variational principles exhibit strong promise for the calculation of guiding center test particle trajectories. The non-canonical Hamiltonian structure of the guiding center equations forms a novel and challenging context for geometric integration. To demonstrate the practical relevance of these methods, a prototypical variational midpoint algorithm is applied to an experimental magnetic equilibrium. The stability characteristics, conservation properties, and implementation requirements associated with the variational algorithms are addressed. Furthermore, computational run time is reduced for large numbers of particles by parallelizing the calculation on GPU hardware.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Hydroelastic response of a floating runway to cnoidal waves

The hydroelastic response of mat-type Very Large Floating Structures (VLFSs) to severe sea conditions, such as tsunamis and hurricanes, must be assessed for safety and survivability. An efficient and robust nonlinear hydroelastic model is required to predict accurately the motion of and the dynamic loads on a VLFS due to such large waves. We develop a nonlinear theory to predict the hydroelastic response of a VLFS in the presence of cnoidal waves and compare the predictions with the linear theory that is also developed here. This hydroelastic problem is formulated by directly coupling the structure with the fluid, by use of the Level I Green-Naghdi theory for the fluid motion and the Kirchhoff thin plate theory for the runway. The coupled fluid structure system, together with the appropriate jump conditions are solved in two-dimensions by the finite-difference method. The numerical model is used to study the nonlinear response of a VLFS to storm waves which are modeled by use of the cnoidal-wave theory. Parametric studies show that the nonlinearity of the waves is very important in accurately predicting the dynamic bending moment and wave run-up on a VLFS in high seas.

42 ENGINEERING↗
DOE OSTI2016

Transport studies in high-performance field reversed configuration plasmas

A significant improvement of field reversed configuration (FRC) lifetime and plasma confinement times in the C-2 plasma, called High Performance FRC regime, has been observed with neutral beam injection (NBI), improved edge stability, and better wall conditioning [Binderbauer et al., Phys. Plasmas 22, 056110 (2015)]. A Quasi-1D (Q1D) fluid transport code has been developed and employed to carry out transport analysis of such C-2 plasma conditions. The Q1D code is coupled to a Monte-Carlo code to incorporate the effect of fast ions, due to NBI, on the background FRC plasma. Numerically, the Q1D transport behavior with enhanced transport coefficients (but with otherwise classical parametric dependencies) such as 5 times classical resistive diffusion, classical thermal ion conductivity, 20 times classical electron thermal conductivity, and classical fast ion behavior fit with the experimentally measured time evolution of the excluded flux radius, line-integrated density, and electron/ion temperature. The numerical study shows near sustainment of poloidal flux for nearly 1 ms in the presence of NBI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Numerical analysis of modified Central Solenoid insert design

The United States ITER Project Office (USIPO) is responsible for fabrication of the Central Solenoid (CS) for ITER project. The ITER machine is currently under construction by seven parties in Cadarache, France. The CS Insert (CSI) project should provide a verification of the conductor performance in relevant conditions of temperature, field, currents and mechanical strain. The US IPO designed the CSI that will be tested at the Central Solenoid Model Coil (CSMC) Test Facility at JAEA, Naka. To validate the modified design we performed three-dimensional numerical simulations using coupled solver for simultaneous structural, thermal and electromagnetic analysis. Thermal and electromagnetic simulations supported structural calculations providing necessary loads and strains. According to current analysis design of the modified coil satisfies ITER magnet structural design criteria for the following conditions: (1) room temperature, no current, (2) temperature 4K, no current, (3) temperature 4K, current 60 kA direct charge, and (4) temperature 4K, current 60 kA reverse charge. Fatigue life assessment analysis is performed for the alternating conditions of: temperature 4K, no current, and temperature 4K, current 45 kA direct charge. Results of fatigue analysis show that parts of the coil assembly can be qualified for up to 1 million cycles. Distributions of the Current Sharing Temperature (TCS) in the superconductor were obtained from numerical results using parameterization of the critical surface in the form similar to that proposed for ITER. Lastly, special ADPL scripts were developed for ANSYS allowing one-dimensional representation of TCS along the cable, as well as three-dimensional fields of TCS in superconductor material. Published by Elsevier B.V.

42 ENGINEERING↗
DOE OSTI2018

Two patterns of PT-symmetry breakdown in a non-numerical six-state simulation

Highlights: • Hot topic (PT symmetry breakdown in quantum theory) studied. • Next to elementary toy model considered. • Characteristic features (like quantum phase transitions) sampled. • Surprisingly, the results are non-numerical. • The method (inverse-function treatment) is brand new. - Abstract: Three-parametric family of non-Hermitian but PT-symmetric six-by-six matrix Hamiltonians H{sup (6)}(x,y,z) is considered. The PT-symmetry remains spontaneously unbroken (i.e., the spectrum of the bound-state energies remains real so that the unitary-evolution stability of the quantum system in question is shown guaranteed) in a non-empty domain D{sup (physical)} of parameters x,y,z. The construction of the exceptional-point (EP) boundary ∂D{sup (physical)} of the physical domain is preformed using an innovative non-numerical implicit-function-construction strategy. The topology of the resulting EP boundary of the spontaneous PT-symmetry breakdown (i.e., of the physical “horizon of stability”) is shown similar to its much more elementary N=4 predecessor. Again, it is shown to consist of two components, viz., of the region of the quantum phase transitions of the first kind (during which at least some of the energies become complex) and of the quantum phase transitions of the second kind (during which some of the level pairs only cross but remain real).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Numerical Study of Vortex Core in Turbulent Swirling Flow

Swirling flow has attracted a great deal of research interest. An application of swirling flow can be found in many engineering areas; such as gas cyclone and heat transfer enhancement. In nuclear engineering, swirling flow is used as heat transfer enhancement in rod bundle by installing grid spacer and mixing vane. Because of its extensive applications, there has been a great number of experimental and numerical investigations of swirling flow. McClusky et al., experimentally studied mapping of the lateral velocity in rod bundle using Laser Doppler Velocimetry (LDV) and Conner et al. predicted the flow field in rod bundle by CFD. To understand the characteristics of swirling flow such as tangential velocity and vortex core swirling flow, some researchers investigated swirling flow in a pipe instead of in rod bundle. This is due to the similarity between swirling flow in a pipe and rod bundle. In recent study, flow mapping of 2D velocity field of single-phase swirling flow had been conducted experimentally by Particle Image Velocimetry (PIV). However, this method needs calibration and complicated. The new measurement technique for flow mapping in swirling flow was done by phased array Ultrasonic Velocity Profiler (PAUVP). Hamdani et al., successfully plotted velocity field and determined the swirling core region (vortex core). This technique is nonintrusive and no need calibration. However, there was a limitation of the measurement because phased array UVP cannot capture 360 deg. velocity field inside the pipe. In addition, tangential velocity in CFD results was not agree well with experiment data due asymmetric problem. To evaluate this limitation of experimental method, the numerical simulations using Computational Fluid Dynamics (CFD) needs to be performed. The objective of this study is to use the commercial CFD code Fluent to model the developing turbulent swirling flow induced by twisted tape inside a straight pipe, and to gain an understanding of such a flow and the turbulence models involved in the simulation. Straight pipe is used instead of rod bundle for a simplification. In this work, the Renormalized Group (RNG) k-{epsilon} turbulence model and Reynolds Stress Model (RSM) turbulence model, are performed to study the phenomena of vortex core in a pipe with twisted tape inserts. Effects of the Reynolds number and asymmetric of twisted tape (inclination, {alpha} = 5 deg. and {alpha} = 20 deg.) on swirling core position are examined using water as the testing fluid. Furthermore, the Navier-Stokes equations will be solved to obtain results that can be validated with the previous experimental results. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Non-local transport in radiation-hydrodynamics codes for ICF by efficient coupling to an external Vlasov–Fokker–Planck code

Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.

Electronic transport↗