Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dimensionless Numbers”

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.

79 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dimensionless parameters to categorize the failure modes of ductile plate perforation

In this work, two dimensionless parameters were derived organically from impact mechanics equations to remove the need for empirical constants. These two dimensionless parameters effectively collapse the ballistic performance of ductile metal plates along two distinct loci, depending on the projectile geometry and the corresponding failure mode. Projectiles that initiate ductile hole enlargement via cavity expansion fall along a linear curve, while shear-plugging failure modes tend to fall along a cubic curve. Results were demonstrated for various projectiles of different calibers impacting aluminum and Ti-6Al-4V alloys, and further explored for more complex target combinations. The main advantage of the approach in this work is to drastically reduce the number of features for input into machine learning algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Drop Interactions with the Conical Shock Structure Generated by a Mach 4.5 Projectile

This work presents measurements of liquid drop deformation and breakup time behind approximately conical shock waves and evaluates the predictive capabilities of low-order models and correlations developed using planar shock experiments. A conical shock was approximated by firing a bullet at Mach 4.5 past a vertical column of water drops with a mean initial diameter of $192$ $\mu$m. The time-resolved drop position and maximum transverse dimension were characterized using backlit stereo images taken at 500 kHz. The gas density and velocity fields experienced by the drops were estimated using a Reynolds-averaged Navier–Stokes simulation of the bullet. Classical correlations predict drop breakup times and deformation in error by a factor of 3 or more. The Taylor analogy breakup (TAB) model predicts deformed drop diameters that agree within the confidence bounds of the ensemble-averaged experimental values using a dimensionless constant $C$ 2 = $2$ compared to the accepted value $C$ 2 = $2/3$. In conclusion, results demonstrate existing correlations are inadequate for predicting the drop response to the three-dimensional relaxation of the flowfield downstream of a conical-like shock and suggest the TAB model results represent a path toward improved predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Error estimates of finite difference methods for the Dirac equation in the massless and nonrelativistic regime

We present four frequently used finite difference methods and establish the error bounds for the discretization of the Dirac equation in the massless and nonrelativistic regime, involving a small dimensionless parameter 0 < ε &NestedLessLess; 1 inversely proportional to the speed of light. In the massless and nonrelativistic regime, the solution exhibits rapid motion in space and is highly oscillatory in time. Specifically, the wavelength of the propagating waves in time is at O(ε), while in space, it is at O(1) with the wave speed at O(ε -1 ). We adopt one leap-frog, two semi-implicit, and one conservative Crank-Nicolson finite difference methods to numerically discretize the Dirac equation in one dimension and establish rigorously the error estimates which depend explicitly on the time step τ, mesh size h, and the small parameter ε. The error bounds indicate that, to obtain the “correct” numerical solution in the massless and nonrelativistic regime, i.e., 0 < ε &NestedLessLess; 1, all these finite difference methods share the same ε-scalability as time step τ = O(ε 3/2 ) and mesh size h = O(ε 1/2 ). A large number of numerical results are reported to verify the error estimates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Detecting cosmic strings with lensed fast radio bursts

Correlated red noise recently reported from pulsar timing observations may be an indication of stochastic gravitational waves emitted by cosmic strings that formed during a primordial phase transition near the grand unification energy scale. Unfortunately, known probes of cosmic strings, namely the cosmic microwave background anisotropies and string lensing of extragalactic galaxies, are not sensitive enough for low dimensionless string tensions of Gμc -2 = 10 -10 – 10 -7 (where the tension μ is the string energy per unit length) that are needed to explain this putative signal. We show that strong gravitational lensing of fast radio bursts (FRBs) by cosmic strings is a potentially unambiguous avenue to probe that range of string tension values. The image pair of string lensing are expected to have identical magnification factor and parity, and have a typical time delay of ~10 2 (Gμc -2 /10 -8 ) 2 s. Here, the unique spectral fingerprint of each FRB, as well as the possibility to detect correlations in the time series of the electric field of the radio waves, will enable verification of the string lensing interpretation. Very-long-baseline interferometry observations can spatially resolve the image pair and provide a lower bound on the string tension based on the image separation. We calculate the FRB lensing rate as a function of the FRB detection number for several different models of the FRB redshift distribution. We find that a survey detecting ~10 5 FRBs, in line with estimates for the detection rate of the forthcoming survey CHORD, can uncover a strong lensing event for a string tension of Gμc -2 ≃ 10 -7 . Larger FRB surveys, such as Phase 2 of the Square Kilometre Array, have the potential to significantly improve the sensitivity on the string tension to Gμc -2 ≃ 10 -9 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Mesoscopic structure of mixed type domain walls in multiaxial ferroelectrics

The structure of a 180° uncharged rotational domain wall in a multiaxial ferroelectric film is studied in the framework of an analytical Landau-Ginzburg-Devonshire (LGD) approach. Finite element modeling (FEM) is used to solve numerically the system of the coupled nonlinear Euler-Lagrange (EL) second-order differential equations for two components of polarization. We show that the structure of the domain wall and corresponding metastable or stable phase of the film are controlled by a single parameter—the dimensionless ferroelectric anisotropy factor μ. We fitted the static profile of a solitary domain wall, calculated by FEM, with kinklike functions for polarization components, and extracted the five μ-dependent parameters from the fitting to FEM curves. The surprisingly high accuracy of the fitting results for two polarization components in the entire μ range allows us to conclude that the analytical functions, which are trial functions in the direct variational method, can be treated as a high-accuracy variational solution of the static EL equations. We further derive exact two-component analytical solutions of the static EL equations for a polydomain 180° domain structure in a multiaxial ferroelectric film. Using these, we derive analytical expressions for the system free energy and analyze its dependence on the film thickness and boundary conditions at the film surfaces. The single-domain state is ground for zero polarization derivative at the surfaces, while the polydomain states minimize the system energy for zero polarization at the surfaces. Counterintuitively, the energy of the polydomain states split into two levels, “0” and “1,” for zero polarization at the surfaces, and each of the levels contains a large number of close-energy sublevels, whose structure is characterized by a different number and type of domain walls. The analytical solutions can become a useful tool for Bayesian analysis of high-resolution scanning transmission electron microscopy images in ferroelectric films.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phonon-informed Neural Thermal Scattering (NeTS) Optimization for Crystalline Graphite and Beryllium Metal

Fast neutrons born from fission lose energy through scattering interactions in the process of slowing-down. As neutrons thermalize to the order of $k$ $b$ $T$ (where $k$ $b$ is the Boltzmann constant, and $T$ is the temperature of the medium), their de Broglie wavelength and energy approaches the order of inter-atomic spacing and quantized lattice vibrations, i.e., phonons. At thermal energies, the thermal scattering law (TSL), i.e., $S$($α, β$), captures crystal binding contributions to the total reaction rate, or cross section. This dimensionless material property describes the energy ($β$) and momentum ($α$) exchanges available in a medium. Currently, $S$($α, β$) is evaluated in the Full Law Analysis Scattering System Hub (FLASSH) code for discrete inputs and stored as ENDF/B File 7 for 0-phonon elastic (MT 2) and n-phonon inelastic (MT 4) processes. Further processing recasts $S$($α, β$) into cumulative distribution functions for sampling post-collision scattering kinematics. In practice, interpolation schemes are employed to access data between tabulated values. An improvement to this juncture of the nuclear data pipeline is supplying cross sections on-the-fly (OTF), as has been developed for the un-resolved resonance region to minimize non-physical interpolation errors. This capability may improve simulation accuracy for accident and transient analyses, where rapidly varying changes in temperature and pressure are difficult to predict beforehand. To do so, deep artificial neural networks (ANNs) can be employed which collapse non-linear, complex data into a lightweight dictionary of neural weights and biases. This has been successfully demonstrated for the hydrogen in light water $S$($α, β$) dataset in the form of a Neural Thermal Scattering (NeTS) module. In this work, the NeTS framework is extended to consider the impact of material-dependent dynamical features on optimal neural pre-processing and architecture design decisions, such as number of neurons per hidden layer, residual skip connections and neural depth. New NeTS modules for crystalline graphite and beryllium metal illuminate a novel correlation between dynamical nonlinearity and optimal neural parametrization when deploying $S$($α, β$) on-the-fly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗