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,495 records · Page 83

Anisotropic fluid dynamical simulations of heavy-ion collisions

We present VAH, a (3+1)–dimensional simulation that evolves the far-from-equilibrium quark-gluon plasma produced in ultrarelativistic heavy-ion collisions with anisotropic fluid dynamics. We solve the hydrodynamic equations on an Eulerian grid using the Kurganov–Tadmor algorithm in combination with a new adaptive Runge–Kutta method. Our numerical scheme allows us to start the simulation soon after the nuclear collision, largely avoiding the need to integrate it with a separate pre-equilibrium dynamics module. We test the code's performance by simulating on the Eulerian grid conformal and non-conformal Bjorken flow as well as conformal Gubser flow, whose (0+1)–dimensional solutions are precisely known. Finally, we compare non-conformal anisotropic hydrodynamics to second-order viscous hydrodynamics in central Pb+Pb collisions and find that the former's longitudinal flow profile responds more consistently to the fluid's gradients along the spacetime rapidity direction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Causal explicit algorithm for heat conduction in a plasma

Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Conjugate heat transfer↗

Leveraging large language models to address data scarcity in machine learning for graphene synthesis

Machine learning in experimental materials science faces significant challenges due to the scarcity of data, which are costly and time-consuming to generate, particularly when relying on in-house experiments. Literature data mining offers a potential solution but introduces issues like mixed data quality, inconsistent formats, and non-uniform reporting of synthesis parameters, resulting in partially missing and heterogeneous features across the dataset. Here, we propose data imputation and feature engineering methods that employ pre-trained large language models (LLMs) to enhance machine learning performance on scarce, heterogeneous datasets, demonstrated on graphene CVD synthesis data and the ML-HydPARK hydrogen storage dataset. GPT models perform data imputation via tailored prompting and semantic normalization of inconsistently reported features through embeddings, for example, to harmonize the complex nomenclature of CVD substrates. Beyond yielding more diverse and richer feature representations than traditional methods such as K-nearest neighbors (KNN) and Multivariate Imputation by Chained Equations (MICE), LLM-based data imputation is evaluated against dataset characteristics and prompting strategies. We vary the level of autonomy granted to the LLM, from generic prompting that leverages pre-trained knowledge for autonomous data generation to data-informed prompting that constrains outputs using target-specific information, and demonstrate which level of autonomy yields superior imputation performance across datasets and feature types. The proposed data engineering methods markedly improve downstream performance; for example, in graphene layer number classification using a support vector machine (SVM), binary accuracy increases from 39% to 65% and ternary accuracy from 52% to 72%. Fine-tuning experiments on both datasets show that combining our proposed LLM-based data imputation and feature encoding methods with numerical machine learning predictors outperforms standalone fine-tuned LLM predictors in data-scarce settings. The proposed strategies emphasize data enhancement techniques rather than refining learning architectures or regularizing loss functions, offering a broadly applicable framework for improving machine learning performance on scarce, inhomogeneous datasets.

Chemical vapor deposition↗

Comparing Poynting flux dominated magnetic tower jets with kinetic-energy dominated jets

Magnetic towers represent one of two fundamental forms of MHD out flows. Driven by magnetic pressure gradients, these flows have been less well studied than magneto-centrifugally launched jets even though magnetic towers may well be as common. We present new results exploring the behavior and evolution of magnetic tower outflows and demonstrate their connection with pulsed power experimental studies and purely hydrodynamic jets which might represent the asymptotic propagation regimes of magneto-centrifugally launched jets. High-resolution AMR MHD simulations (using the AstroBEAR code) provide insights into the underlying physics of magnetic towers and help us constrain models of their propagation. Our simulations have been designed to explore the effects of thermal energy losses and rotation on both towerflows and their hydro counterparts. We find these parameters have significant effects on the stability of magnetic towers, but mild effects on the stability of hydro jets. Current-driven perturbations in the Poynting Flux Dominated (PDF) towers are shown to be amplified in both the cooling and rotating cases. Our studies of the long term evolution of the towers show that the formation of weakly magnetized central jets within the tower are broken up by these instabilities becoming a series of collimated clumps which magnetization properties vary over time. In addition to discussing these results in light of laboratory experiments, we address their relevance to astrophysical observations of young star jets and outflow from highly evolved solar type stars.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Green's function quantum average atom model

A quantum average atom model is reformulated using Green's functions. This allows integrals along the real energy axis to be deformed into the complex plane. The advantage being that sharp features such as resonances and bound states are broadened by a Lorentzian with a half-width chosen for numerical convenience. An implementation of this method therefore avoids numerically challenging resonance tracking and the search for weakly bound states, without changing the physical content or results of the model. A straightforward implementation results in up to a factor of 5 speed-up relative to an optimized orbital based code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Using distributed temperature sensing to detect CO 2 leakage along the injection well casing

The objective of this study is to evaluate the sensitivity of Distributed Temperature Sensing (DTS) data to detect CO 2 leakage along an injection well casing. This paper describes the relationship between the CO 2 leakage rate and temperature response at DTS locations, and the method and numerical model used for understanding such a relationship. The uncertainties in the parameters are propagated to the interpretation of DTS measurements, which may lead to false positive and false negative identifications of CO 2 leakage. We propose to identify CO 2 leakage by analyzing both temperature history plots at selected vertical locations as well as the vertical temperature profiles at different times. The analysis should be combined with numerical simulations for the estimation of leakage rates. In addition, leakage needs to be confirmed by data from a warm-up test (temperature recovery after injection of cold CO 2 is stopped) to minimize the probability of false identifications.

42 ENGINEERING↗

A wavelet-based computational method for solving stochastic Itô–Volterra integral equations

This paper presents a computational method based on the Chebyshev wavelets for solving stochastic Itô–Volterra integral equations. First, a stochastic operational matrix for the Chebyshev wavelets is presented and a general procedure for forming this matrix is given. Then, the Chebyshev wavelets basis along with this stochastic operational matrix are applied for solving stochastic Itô–Volterra integral equations. Convergence and error analysis of the Chebyshev wavelets basis are investigated. To reveal the accuracy and efficiency of the proposed method some numerical examples are included.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A stochastic Galerkin method for the Boltzmann equation with uncertainty

We develop a stochastic Galerkin method for the Boltzmann equation with uncertainty. The method is based on the generalized polynomial chaos (gPC) approximation in the stochastic Galerkin framework, and can handle random inputs from collision kernel, initial data or boundary data. We show that a simple singular value decomposition of gPC related coefficients combined with the fast Fourier-spectral method (in velocity space) allows one to compute the high-dimensional collision operator very efficiently. In the spatially homogeneous case, we first prove that the analytical solution preserves the regularity of the initial data in the random space, and then use it to establish the spectral accuracy of the proposed stochastic Galerkin method. Several numerical examples are presented to illustrate the validity of the proposed scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The PAC-MAN model: Benchmark case for linear acoustics in computational physics

Benchmark cases in the field of computational physics, on the one hand, have to contain a certain complexity to test numerical edge cases and, on the other hand, require the existence of an analytical solution, because an analytical solution allows the exact quantification of the accuracy of a numerical simulation method. This dilemma causes a need for analytical sound field formulations of complex acoustic problems. A well known example for such a benchmark case for harmonic linear acoustics is the “Cat's Eye model”, which describes the three-dimensional sound field radiated from a sphere with a missing octant analytically. In this paper, a benchmark case for two-dimensional (2D) harmonic linear acoustic problems, viz., the “PAC-MAN model”, is proposed. The PAC-MAN model describes the radiated and scattered sound field around an infinitely long cylinder with a cut out sector of variable angular width. While the analytical calculation of the 2D sound field allows different angular cut-out widths and arbitrarily positioned line sources, the computational cost associated with the solution of this problem is similar to a 1D problem because of a modal formulation of the sound field in the PAC-MAN model. - Highlights: • A benchmark case for 2D harmonic linear acoustic problems is proposed. • The PAC-MAN geometry is used as emitter/scatterer in this benchmark. • An analytic formulation of radiated and scattered sound is derived. • Line and disk sources, plane waves, and surface vibrations are considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enforcing the Courant–Friedrichs–Lewy condition in explicitly conservative local time stepping schemes

In this study, an optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves . We use an experimental fluid dynamics code that allows cubic "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a condition on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty.

79 ASTRONOMY AND ASTROPHYSICS↗

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On numerical energy conservation for an implicit particle-in-cell method coupled with a binary Monte-Carlo algorithm for Coulomb collisions

We report conventional particle-in-cell (PIC) methods suffer from enhanced numerical heating (explicit PIC) or cooling (semi-implicit PIC) when coupled with a binary Monte-Carlo algorithm for Coulomb collisions. In this work, a fully-implicit θ-PIC scheme (with adjustable time-biasing parameter 1/2 ≤ θ ≤) is considered. The discrete change in energy of a closed system after a time step for this scheme scales with (1/2 - θ)C θ , where C θ is a positive definite quantity that depends on the frequency spectrum of the energy in the fields. Collisions lead to additional energy in the field fluctuations associated with high-frequency light waves produced by a numerical Bremsstrahlung process, which can result in a large increase in the numerical cooling rate for θ > 1/2. However, for θ = 1/2, energy is exactly conserved. The energy in the field fluctuations on long time scales agrees with that calculated using the equipartition theorem for a classical system in thermodynamic equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Process intensification in hydrothermal liquefaction of biomass: A review

Hydrothermal liquefaction (HTL) presents a promising pathway for converting wet biomass resources into biofuels, offering significant advantages over conventional methods. However, numerous challenges must be addressed for HTL scale-up, including energy provision for the endothermic process, heat and mass transfer limitations, slurry concentration and pumpability challenges, char and coke formation, and continuous phase separation. This review explores key strategies such as autothermal HTL, which improves process efficiency and reduces external energy requirements by coupling exothermic and endothermic reactions within the same reactor, thereby simplifying reactor design and reducing operational costs. Additionally, multistage HTL processes are highlighted for their ability to optimize biocrude quality and yield by fractionating biomass conversion stages, resulting in higher energy returns on investment and better-quality biocrude. Solvothermal HTL and integration techniques for aqueous phase are also discussed. Furthermore, the HTL patent landscape is discussed to provide insights into current technological advancements. This review aims to offer a comprehensive understanding of process intensification in HTL, highlighting innovative solutions to enhance the efficiency and scalability of the process for sustainable biofuel production.

Biocrude oil↗