Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerics”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Numerical methods for fractional Fokker–Planck equation with multiplicative Marcus Lévy noises

The Fokker–Planck equation (FPE) is an important deterministic tool for investigating stochastic dynamical systems. In this paper, we consider the space-time fractional FPE driven by multiplicative Marcus Lévy noises. Efficient numerical schemes are presented to solve the equations. Stability and convergence of the methods are also discussed. We give some numerical experiments to validate our schemes, and examine the effects of parameters on solutions. Additionally, we analyze the maximal likely trajectories and the critical time for the change of the most probability location.

Mathematics↗
DOE OSTI2022

A Comprehensive Numerical and Experimental Study for the Passive Thermal Management in Battery Modules and Packs

Cooling plates in battery packs of electric vehicles play critical roles in passive thermal management systems to reduce risks of catastrophic thermal runaway. In this work, a series of numerical simulations and experiments are carried out to unveil the role of cooling plates (both between cells and a bottom plate parallel to the cell stack) on the thermal behavior of battery modules and packs under nail penetrations. First, we investigated the role of side cooling plates on the thermal runaway propagation mitigation in battery modules (1S3P) and packs (3S3P) by varying the key parameters of the side cooling plates, such as plate thicknesses, thermal contact resistances, and materials. Then, three important factors for passive thermal management systems are identified: (i) thermal mass of side cooling plates, (ii) interfacial thermal contact resistances, and (iii) the effective heat transfer coefficients at exterior surfaces. The roles of bottom cooling plates on thermal runaway propagation mitigation in 1S3P and 1S5P battery modules are numerically investigated by comparing the thermal behavior of the modules with only side cooling plates and with both side and bottom cooling plates.

25 ENERGY STORAGE↗
DOE OSTI2021

Reduced Order Modeling for Accelerating Numerical Simulations (ROMANS)

SAND2021-15064 O ROMANS, or Reduced Order Modeling for Accelerating Numerical Simulations, is an open source software for accelerating numerical simulations. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Martin, Shawn↗
DOE OSTI2024

FTTN: Feature-Targeted Testing for Numerical Properties of NVIDIA & AMD Matrix Accelerators

FTTN is a test suite to evaluate the numerical behaviors of matrix accelerators of GPUs (NVIDIA Tensor Cores and AMD Matrix Cores) in a quick and simple setting. Matrix accelerators are heavily used in today's computationally intense applications to speed up matrix multiplications. This test suite provides a comprehensive study on the numerical behaviors of these accelerators, including support for subnormals, rounding modes, extra precision bits and FMA features. Is there

Laguna Peralta, Ignacio↗
DOE OSTI2024

Efficient general method for numerically modeling laser pulse propagation, overlap, and lifetime effects in amplifiers

An efficient numerical time-dependent general method is developed to address incoherent pulse overlap and lifetime effects in laser amplifiers. The alternating propagation-population laser energetics method (APPLE) has been validated against a semi-discrete coupled rate equation numerical method (SDRE) and analytic formalisms in bounding cases. APPLE is based on decoupled rates applied to a time-dependent framework where both space-time-dependent populations and pulse energetics are consistently updated in each time step. A significant advantage of APPLE lies in its conceptual simplicity, ease of implementation, and relatively small computational cost. SDRE tracks the populations through coupled rates and uses the method of lines to discretize the hyperbolic partial differential transport equations allowing for use of ordinary differential equation solvers. With reasonably sized mesh, we report both energetic and power pulse shape relative differences on the order of one percent between the models over a large range of initial conditions.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Experimental and numerical evaluation of laser shinethrough in polystyrene at 3 ω and 4 ω laser frequencies

In laser-direct-drive inertial confinement fusion (ICF), laser shinethrough is the process associated with laser propagation through the target plastic ablator, polystyrene is often used in ICF experiments, and leads to energy deposition in the target until the critical density is established at the ablator front surface. This process causes material modifications that may change the shock propagation and lead to density and pressure perturbations at the interface between the polystyrene and cryogenic deuterium-tritium (DT) layer. This work reports on experimental observations and numerical investigations of laser shinethrough at the third and the fourth harmonic of ICF class laser systems (1 ω corresponds to a wavelength of ∼1.06 μm ). The experiments evidence that, at 4 ω , the laser is absorbed over a distance of about a few microns and laser shinethrough is completely eliminated in contrast with the behavior at 3 ω . These experimental observations are reproduced by numerical simulations. The laser fluence required to reach the critical density at the ablator front surface is also found to be more than one order of magnitude higher at the third harmonic. Finally, modeling suggests that the back surface of the ablator exhibits ion temperatures greater than the melting temperature of cryogenic DT and electron pressure larger than the bulk and shear moduli of cryogenic DT for the case of using a 3 ω laser.

Pineau, A. [Univ. of Rochester, NY (United States)↗
DOE OSTI2022

NUMERICAL INVESTIGATION OF AIR DEHUMIDIFICATION THROUGH WATER DROPLETS DIELECTROPHORESIS

This paper numerically investigated the dehumidification of moist airflow in a converging rectangular duct with electrically enhanced condensation. The aim was to develop a condensation model that predicted water vapor condensation by employing the dielectrophoresis principle. The charged droplets were injected into the computational domain to produce a gradient electric field. The dielectrophoretic interaction between the droplets and vapor molecules of the humid air resulted in local water vapor condensation at the vapor and liquid droplet interface. This phenomenon is described in detail in this paper. A hybrid Eulerian-Lagrangian solver SprayFoam was developed based on OpenFOAM® to simulate condensation in two-phase gas-liquid mixtures. The following developments were made based on the standard compressible multiphase flow solver SprayFoam in OpenFOAM®: (1) Eulerian solver for gas phase (2) multi-component species transport, and (3) Lagrangian solver for gas-droplet two-phase flows and sub-models for liquid droplets. The newly developed numerical model was experimentally validated with data from a series of tests conducted in the authors' laboratory. The results showed that the simulations followed the same trends as the data, and the model predicted the condensation due to the electro-spray injection process in the air stream. The injected water droplets, which are electrically charged to their Rayleigh limits, increased in size while wiping out the humidity from the air. The simulations indicated that the dehumidification was about 1% for 0.5 cubic feet per minute airflow rate. Scaling up to larger flows is a future follow-up work.

Yel Mahi, *Maliha↗
DOE OSTI2023

Numerical Modeling of Air-Blast Suppression as a Function of Explosive-Charge Burial Depth

As a chemical explosion is buried, the mechanism for acoustic wave generation transitions from fully gas-generated at the surface to completely spall-induced at full containment depth. The fully gas-generated and completely spall-induced signals in the acoustic waveform are well described; however, the transition between these two end-members eludes numerical modeling because of the complex phenomena that are involved. The phenomena of crater formation and explosive cloud evolution are simulated using an Eulerian hydrocode that incorporates geomaterials with strength and porosity. Having accurately modeled these phenomena, we can confidently predict the propagation and relative strength of the gas-generated and spall-induced pulses in the recorded acoustic waveform. The numerical predictions agree with observations from the historical Stagecoach experiment as well as modern recordings from the Source Physics Experiment. In particular, the peak pressure p generated by an explosion is initially due to the gas-generated mechanism and decays with scaled depth of burial d s (depth d scaled by the cube-root of explosive yield w 1/3 ) as exp(-d s ) but then transitions near a scaled depth of 6 m/ton 1/3 to the spall-generated mechanism in which the decay is d$_{s}^{-7/4}$. This decay form is related to the strong ground-motion attenuation relationship that affects spall strength. So these results can improve seismoacoustic inverse models for the explosive source that need to account for the gas-generated and spall-induced signals and their effect on peak pressures and other acoustic signal features.

58 GEOSCIENCES↗
DOE OSTI2022

Assessment of Lost Circulation Material Particle-Size Distribution on Fracture Sealing: A Numerical Study

Summary Lost circulation materials (LCMs) are essential to combat fluid loss while drilling and may put the whole operation at risk if a proper LCM design is not used. The focus of this research is understanding the function of LCMs in sealing fractures to reduce fluid loss. One important consideration in the success of fracture sealing is the particle-size distribution (PSD) of LCMs. Various studies have suggested different guidelines for obtaining the best size distribution of LCMs for effective fracture sealing based on limited laboratory experiments or field observations. Hence, there is a need for sophisticated numerical methods to improve the LCM design by providing some predictive capabilities. In this study, computational fluid dynamics (CFD) and discrete element methods (DEM) numerical simulations are coupled to investigate the influence of PSD of granular LCMs on fracture sealing. Dimensionless variables were introduced to compare cases with different PSDs. We validated the CFD-DEM model in reproducing specific laboratory observations of fracture-sealing experiments within the model boundary parameters. Our simulations suggested that a bimodally distributed blend would be the most effective design in comparison to other PSDs tested here.

Engineering↗
DOE OSTI2020

Analytic and numerical bootstrap of CFTs with $O(m)\times O(n)$ global symmetry in 3D

Motivated by applications to critical phenomena and open theoretical questions, we study conformal field theories with O(m)\times O(n) O ( m ) × O ( n ) global symmetry in d=3 d = 3 spacetime dimensions. We use both analytic and numerical bootstrap techniques. Using the analytic bootstrap, we calculate anomalous dimensions and OPE coefficients as power series in \varepsilon=4-d ε = 4 − d and in 1/n 1 / n , with a method that generalizes to arbitrary global symmetry. Whenever comparison is possible, our results agree with earlier results obtained with diagrammatic methods in the literature. Using the numerical bootstrap, we obtain a wide variety of operator dimension bounds, and we find several islands (isolated allowed regions) in parameter space for O(2)\times O(n) O ( 2 ) × O ( n ) theories for various values of n n . Some of these islands can be attributed to fixed points predicted by perturbative methods like the \varepsilon ε and large- n n expansions, while others appear to arise due to fixed points that have been claimed to exist in resummations of perturbative beta functions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗
DOE OSTI2021

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Numerical Simulations of Non-Proliferation Experiments

We have implemented a set of numerical simulations to numerically investigate the outcomes of non-proliferation experiments. The high explosive (HE) in the simulations is in the form of a cylinder. The two-dimensional simulations are carried out in the r-z cylindrical coordinate, and the three-dimensional simulation is in the Cartesian coordinate. From the simulation, the cavity formed after detonation stops its expansion after about 0.15 sec. At the time the cavity stops expanding, the cavity shape seems quite spherical although the initial HE is in the shape of cylinder, but the radius in the vertical direction, which is about 19.7m, is a little less than the one in the horizonal direction, which is about 21.5m. Although there is no layer of air between the HE and tuff in real experiments, in addition to the simulations without air, we add a layer of air in one simulation to assess any possible influence of the air on the shape of the cavity.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

Bringing randomized algorithms to mainstream numerical linear algebra

Numerical linear algebra (NLA) underpins huge swaths of computational science and engineering. For scientists and engineers to make the most of the DOE’s computing resources, it is essential that they have access to high-performance implementations of algorithms with best-in-class scalability and reliability. Despite this, prevailing NLA libraries have little to no support for breakthrough algorithms from the field of randomized numerical linear algebra (RandNLA) that have been developed over the past twenty years. The goal of this LDRD was to break a log-jam that had prevented broad adoption of RandNLA. Our work had two thrusts. The first was to develop RandBLAS: a trustworthy and high-performance C++ library for randomized dimension reduction (an operation widely known as sketching). The second was the development of a novel randomized algorithm for computing a challenging type of matrix decomposition known as Householder QR with column pivoting (Householder QRCP). In this one-year late-start LDRD we successfully delivered RandBLAS 1.0 and new CPU and GPU codes for Householder QRCP. RandBLAS has extensive documentation at https://randblas.readthedocs.io/en/stable/. Papers on RandBLAS and and our high-performance QRCP codes are forthcoming.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A numerical investigation of coal particle modeling for the inlet section of a pressurized oxy-combustion burner

Concerns over climate change have led to numerous efforts to develop low-carbon technologies, and pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in the power industry. Designing an effective burner for such technologies is vital for the development of these and other coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the pressurized oxy-combustor burner is close to or even higher than the maximum limit for commercial CFD codes (e.g., ANSYS FLUENT). At these high particle volume fractions, the interactions among particles, fluid flow and wall need to be re-evaluated. The present computational work is a first step in a systematic analysis of the influence of various parameters, like method of particle release, release location, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In the POC process, pulverized coal is burned under elevated pressure in an O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing Ansys FLUENT commercial sodtware using the Reynolds-Averaged Navier-Stokes (RANS) approach. It is revealed that for this pilot-scale, pressurized burner, velocity profiles in the near-wall region exhibit some anomalies. In order to investigate the influence of particle loading in the near-wall region, particle release location was investigated. The numerical simulations incorporate the coupling between the turbulent flow and the particles. A sensitivity investigation of particle release location found that by tuning the release location, the velocity profile can be consistent with the pure gas flow velocity profile. More importantly, the particle releasing location also affects flame stability. Particle size was also found to have a significant impact on particle trajectory, flame stability and temperature. Finally, Large Eddy Simulations (LES) were performed and compared with the results from two-dimensional RANS.

Li, Lei↗
DOE OSTI2023

Numerical Investigation of Wall-Cooling Effect on Aero-Optical Distortions for Hypersonic Boundary Layer

Compressible wall-modeled large-eddy simulations of Mach 8 turbulent boundary-layer flows over a flat plate were carried out for the conditions of the hypersonic wind tunnel at Sandia National Laboratories. The simulations provide new insight into the effect of wall cooling on the aero-optical path distortions for hypersonic turbulent boundary-layer flows. Four different wall-to-recovery temperature ratios, 0.3, 0.48, 0.71, and 0.89, are considered. Despite the much lower grid resolution, the mean velocity, temperature, and resolved Reynolds stress profiles from the simulation for a temperature ratio of 0.48 are in good agreement with those from a reference direct numerical simulation. The normalized root-mean-square optical path difference obtained from the present simulations is compared with that from reference direct numerical simulations, Sandia experiments, as well as predictions obtained with a semi-analytical model by Notre Dame University. Here the present analysis focuses on the effect of wall cooling on the wall-normal density correlations, on key underlying assumptions of the aforementioned model such as the strong Reynolds analogy, and on the elevation angle effect on the optical path difference. Wall cooling is found to increase the velocity fluctuations and decrease the density fluctuations, resulting in an overall reduction of the normalized optical path distortion. Compared to the simulations, the basic strong Reynolds analogy overpredicts the temperature fluctuations for cooled walls. Also different from the strong Reynolds analogy, the velocity and temperature fluctuations are not perfectly anticorrelated. Finally, as the wall temperature is raised, the density correlation length, away from the wall but inside the boundary layer, increases significantly for beam paths tilted in the downstream direction.

42 ENGINEERING↗
DOE OSTI2020

Numerical Considerations in the Modeling of a High Explosive Cylinder Experiment Using an ALE Continuum Mechanics Code

A series of experiments involving the detonation of PBX 9501 encased in a copper cylinder are modeled with the objective of evaluating a proposed set of phenomenological parameters for the Wescott–Stewart–Davis reactive burn model. The numerical analysis is conducted using the Los Alamos continuum mechanics code FLAG. Numerical considerations pertaining to various aspects of modeling the experiments using FLAG are discussed. It is shown that use of the proposed set of phenomenological parameters results in predictions of free-surface velocity that match empirically measured velocities reasonably well.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Gravitational Fluctuations as an Alternative to Inflation III. Numerical Results

Power spectra play an important role in the theory of inflation, and their ability to reproduce current observational data to high accuracy is often considered a triumph of inflation, largely because of a lack of credible alternatives. In previous work we introduced an alternative picture for the cosmological power spectra based on the nonperturbative features of the quantum version of Einstein’s gravity, instead of currently popular inflation models based on scalar fields. The key ingredients in this new picture are the appearance of a nontrivial gravitational vacuum condensate (directly related to the observed cosmological constant), and a calculable renormalization group running of Newton’s G on cosmological scales. More importantly, one notes the absence of any fundamental scalar fields in this approach. Results obtained previously were largely based on a semi-analytical treatment, and thus, while generally transparent in their implementation, often suffered from the limitations of various approximations and simplifying assumptions. In this work, we extend and refine our previous calculations by laying out an updated and extended analysis, which now utilizes a set of suitably modified state-of-the-art numerical programs (ISiTGR, MGCAMB and MGCLASS) developed for observational cosmology. As a result, we are able to remove some of the approximations employed in our previous studies, leading to a number of novel and detailed physical predictions. These should help in potentially distinguishing the vacuum condensate picture of quantum gravity from that of other models such as scalar field inflation. Here, besides the matter power spectrum P m (k) , we work out, in detail, predictions for what are referred to as the TT, TE, EE, BB angular spectra, as well as their closely related lensing spectra. However, the current limited precision of observational data today (especially on large angular scales) does not allow us yet to clearly prove or disprove either set of ideas. Nevertheless, by exploring in more details the relationship between gravity and cosmological matter and radiation both analytically and numerically, together with an expected future influx of increasingly accurate observational data, one can hope that the new quantum gravitational picture can be subjected to further stringent tests in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗