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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Fermions in Two Dimensions: Scattering and Many-Body Properties

Ultracold atomic Fermi gases in two dimensions (2D) are an increasingly popular topic of research. The interaction strength between spin-up and spin-down particles in two-component Fermi gases can be tuned in experiments, allowing for a strongly interacting regime where the gas properties are yet to be fully understood. We have probed this regime for 2D Fermi gases by performing T = 0 ab initio diffusion Monte Carlo calculations. The many-body dynamics are largely dependent on the two-body interactions; therefore, we start with an in-depth look at scattering theory in 2D. We show the partial-wave expansion and its relation to the scattering length and effective range. Then, we discuss our numerical methods for determining these scattering parameters. Here, we close out this discussion by illustrating the details of bound states in 2D. Transitioning to the many-body system, we also use variationally optimized wave functions to calculate ground-state properties of the gas over a range of interaction strengths. We show results for the energy per particle and parametrize an equation of state. We then proceed to determine the chemical potential for the strongly interacting gas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Strain Field Analysis in Electronic Components by ESPI: Bad Thermal Contact and Damage Evaluation

Designing electronic components and devices is a complex task that must include also considerations about thermo-mechanical behavior of the component itself. This aspect, in fact, can affect the reliability of the device as a consequence of the thermal stresses that are introduced in the component while it is working. Stress concentration can arise, in fact, due to different thermal expansion coefficients of the materials involved. Determining the thermo-mechanical response by numerically methods is not simple and not always possible so that experimental methods are preferable. In this paper, specifically, the electronic speckle pattern interferometry technique was adopted to measure full-field strain of the component during exercise. By using an optical technique it is possible to get information about the behavior of the component without being in contact with the component itself so that the dissipation coefficient is not altered. The system was able to detect gradients of deformation in the component itself connected to the different distribution of internal current in the component. Moreover, it was observed, that the presence of bad thermal contact between the package and the heat sink can be revealed because it results in a different thermos-mechanical behavior of the component. Finally, a critical defect was introduced in the component and it was observed the way the presence of damage affects the in-plane displacement of the component. Very huge differences were observed leading to the consideration that this kind of approach could be conveniently adopted also as a damage detection tool.

42 ENGINEERING↗
DOE OSTI2018

Thickness-Dependent Current Density and Flux Distribution in Thin Current-Carrying Superconducting Films Exposed to a Magnetic Field

With the rapid development of high-temperature superconductors, it is of great significance to get the precise current density and flux distribution within thin high-temperature superconducting films subjected to a transport current and an applied magnetic field. The transport current distribution and flux density in thin high-temperature superconducting films are calculated by a numerical method based on the Kim model and exponential model in this paper. The influences of transport current, applied magnetic field, width, and thickness of a superconducting film on the current distribution are discussed. The results reveal that the thickness has a significant effect on the critical current density of superconducting films.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Interaction of Harmonic Longitudinal Shear Waves with V-Shaped Inclusions

We solve the problem of determination of the stressed state in the vicinity of a tunnel rigid inclusion whose cross section consists of two segments originating out from the same point. The inclusion is placed in an infinite elastic medium with propagating plane harmonic longitudinal shear waves. The problem is reduced to the solution of a system of two singular integral equations with fixed singularities. For the approximate solution of this system, we apply a numerical method that takes into account the true asymptotics of the unknown functions and is based on the use of special quadrature formulas for the evaluation of singular integrals.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Design of all-optical XOR and XNOR logic gates based on Fano resonance in plasmonic ring resonators

In this paper, compact all-optical XOR and XNOR gates based on different configurations of metal–insulator–metal plasmonic ring resonators (PRRs) have been proposed. Square and octagon-shaped rings have been used in this case. The logic gates have been simulated using the two-dimensional finite-difference time-domain numerical method, with a conventional perfectly matched layer as the absorbing boundary condition of the area under simulation. The phenomenon of Fano resonance is employed by the proposed gates to excite ON/OFF states. As a result, a large value of contrast ratio (C.R.) is obtained. The results show that the values of C.R. of XNOR and XOR gates for square-shaped PRR are 22.66 and 22.9 dB, respectively and the values of C.R. of XNOR and XOR gates for octagon-shaped PRR are 23.01 and 23.52 dB, respectively. Also, it is shown that the octagon-shaped gates have higher transmission ratios than other proposed configurations. The proposed optical logic gates can be used as key components in optical communications and for designing compact plasmonic devices.

42 ENGINEERING↗
DOE OSTI2019

Analyzing honeycomb photonic crystal waveguides by Dirichlet-to-Neumann maps

As an analog to graphene, honeycomb photonic crystals (PhCs) have attracted a great deal of interest in recent years. The additional degrees of freedom in a honeycomb PhC are useful in designing and optimizing waveguides, can be used to control the symmetry of the system and to realize novel optical devices. In this paper, we present an efficient numerical method for analyzing two-dimensional honeycomb PhC waveguides. Our method is based on a special supercell that covers one period of the honeycomb PhC waveguide and the so-called Dirichlet-to-Neumann map of the supercell. The method gives rise to a linear eigenvalue value problem with relatively small matrices, and it is validated by comparing with previous works. As an application of the method, we calculate edge modes for a few different PhCs. For a honeycomb PhC with dielectric rods in air, we found an edge mode at the zigzag edge, without using special modified rods at the edge. For a honeycomb PhC with air holes, we found a new edge mode when the PhC is bounded by a perfectly electric conductor.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Nanoscale plasmonic detector of wave intensity difference and uni-directional waveguide

In this paper, a new nanoscale plasmonic differential detector of wave intensities is proposed. This structure can be useful for calibration applications and for finding out if a wave intensity is high enough to pass a threshold level or not. Next, a nanoscale plasmonic uni-directional waveguide is presented. The distinction of this waveguide with similar cases is that the one-way wave-guiding is provided by adding a properly located controlling source, without employing any non-reciprocal material. This feature can be used in many applications that only uni-directional wave propagation is desirable. Both structures have been designed in metal–insulator–metal structure including cavity resonators. The operations are based on the coupling of the wave from input waveguides to the cavity resonators and decoupling to output waveguides in the cavity resonance wavelengths. The simulation results of performance of the structures, as elements applicable to photonic integrated circuits, have been obtained using the numerical method of the finite difference time domain.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Generation of spiraling Bessel beams from dark/antidark Gaussian beams diffracted by a curved fork-shaped hologram

In this paper, we present a theoretical study on the generation of spiraling Bessel beams by illuminating a curved fork-shaped hologram (CFH) by spatially coherent dark/antidark Gaussian beams. The amplitude expression of the generated spiraling Bessel beam is obtained by using the framework of the Fresnel diffraction approach and by means of the stationary phase method. Numerical calculations show that upon propagating, the longitudinal intensity of the generated field reaches a threshold maximum value at a distance z{sub max}. It is demonstrated that z{sub max} increases with the increasing of the topological charge p and the order m. The result corresponding to the incident Gaussian beam diffracted by CFH is obtained as a particular case of this study.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A Multi-scale Mathematical Model of Growth and Coalescence of Bubbles Beneath the Anode in an Aluminum Reduction Cell

Modeling of bubble shapes is challenging because of the wide range of length scales in aluminum reduction cells. A 3D multi-scale mathematical model was developed to understand the nucleation, growth and coalescence of bubbles beneath the anode and to investigate the transition of bubbles from the micro- to macro level. The motion of micro-bubbles is examined using the discrete bubble model (DBM) within a Lagrangian reference. An algorithm for the transition from discrete micro-bubbles to large bubbles, which are fully resolved by the volume of fluid (VOF) approach, is achieved using a user-defined function. The two-way coupling between discrete bubbles and continuous fluids is achieved by inter-phase momentum exchange. The model involves three kinds of bubble coalescence: micro-bubble coalescence is taken into account by the DBM; large bubbles swallowing up micro-bubbles are solved by the discrete-continuum transition model; the coalescence between large bubbles is handled by the VOF method. Numerical results show that the coverage and thickness of bubbles agree well with the experimental data in the literature. Bubble thickness stops increasing when the bubble elongates with the anode bottom, which is 4.0 to 4.5 mm. Meanwhile, two asymmetrical thick heads can be observed under the anode. Bubble release frequency increases with increasing current density.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine Learning Aided Modeling of Granular Materials: A Review

Artificial intelligence (AI) has become a buzzy word since Google’s AlphaGo beat a world champion in 2017. In the past five years, machine learning as a subset of the broader category of AI has obtained considerable attention in the research community of granular materials. This work offers a detailed review of the recent advances in machine learning-aided studies of granular materials from the particle-particle interaction at the grain level to the macroscopic simulations of granular flow. This work will start with the application of machine learning in the microscopic particle-particle interaction and associated contact models. Then, different neural networks for learning the constitutive behaviour of granular materials will be reviewed and compared. Finally, the macroscopic simulations of practical engineering or boundary value problems based on the combination of neural networks and numerical methods are discussed. We hope readers will have a clear idea of the development of machine learning-aided modelling of granular materials via this comprehensive review work.

42 ENGINEERING↗
DOE OSTI2016

Spatial Ion Peak Compression and its Utility in Ion Mobility Spectrometry

A novel concept for ion spatial peak compression is described, and discussed primarily in the context of ion mobility spectrometry (IMS). Using theoretical and numerical methods, the effects of using non-constant (e.g., linearly varying) electric fields on ion distributions (e.g., an ion mobility peak) is evaluated both in the physical and temporal domains. The application of linearly decreasing electric field in conjunction with conventional drift field arrangements is shown to lead to a reduction in IMS physical peak width. When multiple ion packets in a selected mobility window are simultaneously subjected to such fields, there is ion packet compression, i.e., a reduction in peak widths of all species. This peak compression occurs with a modest reduction of resolution, but which can be quickly recovered as ions drift in a constant field after the compression event. Compression also yields a significant increase in peak intensities. In addition, approaches for peak compression in traveling wave IMS are also discussed. Ion mobility peak compression can be particularly useful for mitigating diffusion driven peak spreading over very long path length separations (e.g., in cyclic multi-pass arrangements), and for achieving higher S/N and IMS resolution over a selected mobility range.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2022

Well-Balanced Second-Order Convex Limiting Technique for Solving the Serre–Green–Naghdi Equations

In this article, we introduce a numerical method for approximating the dispersive Serre–Green–Naghdi equations with topography using continuous finite elements. The method is an extension of the hyperbolic relaxation technique introduced in Guermond et al. (J Comput Phys 450:110809, 2022). It is explicit, second-order accurate in space, third-order accurate in time, and is invariant-domain preserving. It is also well balanced and parameter free. Special attention is given to the convex limiting technique when physical source terms are added in the equations. The method is verified with academic benchmarks and validated by comparison with laboratory experimental data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Shear stress- and line length-dependent screw dislocation cross-slip in FCC Ni

Screw dislocation cross-slip is important in dynamic recovery of deformed metals. A mobile screw dislocation segment can cross slip to annihilate an immobile screw dislocation segment with opposite Burgers vector, leaving excess dislocations of one kind in a crystal. Previous studies have found that the cross-slip process depends on both the local stress state and dislocation line length, yet a quantitative study of the combined effects of these two factors has not been conducted. In this work, we employ both dynamic concurrent atomistic-continuum (CAC) [L. Xiong, G. Tucker, D.L. McDowell, Y. Chen, J. Mech. Phys. Solids 59 (2011) 160–177] and molecular dynamics simulations to explore the shear stress- and line length-dependent screw dislocation cross-slip in face-centered cubic Ni. It is demonstrated that the CAC approach can accurately describe the 3-D cross-slip process at a significantly reduced computational cost, as a complement to other numerical methods. In particular, we show that the Fleischer (FL) [R.L. Fleischer, Acta Metall. 7 (1959) 134–135] type cross-slip, in which a stair-rod dislocation is involved, can be simulated in the coarse-grained domain. Our simulations show that as the applied shear stress increases, the cross-slip mechanism changes from the Friedel-Escaig (FE) [B. Escaig, J. Phys. 29 (1968) 225–239] type to the FL type. In addition, the critical shear stress for both cross-slip mechanisms depends on the dislocation line length. Moreover, the cross-slip of a screw dislocation with a length of 6.47 nm analyzed using periodic boundary conditions occurs via only the FL mechanism, whereas a longer dislocation with length of 12.94 nm can cross-slip via either the FE or FL process in Ni subject to different shear stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Towards improved speed and accuracy of laser powder bed fusion simulations via representation of multiple time scales

Here, we report the result of investigations into alternative representations of the temporal evolution of the laser powder bed fusion (LPBF) process. In order to provide the reader with sufficient context for the following discussion, the modeling challenge and its underlying cause, as well as current attempts to provide a satisfactory solution, will be discussed. Next, the author will describe two numerical methods (multirate time integration and parallel-in-time) to better represent the temporal scales of the problem. Finally, the results of applying these methods to a two-dimensional version of the heat transfer portion of the LPBF problem will be presented, reporting on both convergence and performance behavior, which indicate the possibility of significant speedup of solution calculations (40–100 times), along with well-characterized effects on solution accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

A smoothed particle hydrodynamics model for droplet and film flow on smooth and rough fracture surfaces

Flow on fracture surfaces has been identified by many authors as an important flow process in unsaturated fractured rock formations. Given the complexity of flow dynamics on such small scales, robust numerical methods have to be employed in order to capture the highly dynamic interfaces and flow intermittency. In this work we present microscale free-surface flow simulations using a three-dimensional multiphase Smoothed Particle Hydrodynamics (SPH) code. Pairwise solid-fluid and fluid-fluid interaction forces are used to control the wetting behavior and cover a wide range of static and transient contact angles as well as Reynolds numbers encountered in droplet flow on rock surfaces. We validate our model via comparison with existing empirical and semi-analyical solutions for droplet flow. We use the model to investigate the occurence of adsorbed trailing films of droplets under various flow conditions and its importance for the flow dynamics when films and droplets coexist. We show that flow velocities are higher on prewetted surfaces covered by a thin film which is qualitatively attributed to the enhanced dynamic wetting and dewetting at the trailing and advancing contact line.

unsaturated flow, fractured rocks, smoothed partic↗
DOE OSTI2018

Pore scale study of multiphase multicomponent reactive transport during CO 2 dissolution trapping

Solubility trapping is crucial for permanent CO 2 sequestration in deep saline aquifers. For the first time, a pore-scale numerical method is developed to investigate coupled scCO 2 -water two-phase flow, multicomponent (CO 2 (aq), H + , HCO 3 – , CO 3 2 – and OH – ) mass transport, heterogeneous interfacial dissolution reaction, and homogeneous dissociation reactions. Pore-scale details of evolutions of multiphase distributions and concentration fields are presented and discussed. Time evolutions of several variables including averaged CO 2 (aq) concentration, scCO 2 saturation, and pH value are analyzed. Specific interfacial length, an important variable which cannot be determined but is required by continuum models, is investigated in detail. Mass transport coefficient or efficient dissolution rate is also evaluated. The pore-scale results show strong non-equilibrium characteristics during solubility trapping due to non-uniform distributions of multiphase as well as slow mass transport process. Complicated coupling mechanisms between multiphase flow, mass transport and chemical reactions are also revealed. Lastly, effects of wettability are also studied. The pore-scale studies provide deep understanding of non-linear non-equilibrium multiple physicochemical processes during CO 2 solubility trapping processes, and also allow to quantitatively predict some important empirical relationships, such as saturation-interfacial surface area, for continuum models.

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

From Bessel beam to complex-source-point cylindrical wave-function

This investigation shows that a scalar Bessel beam can be transformed into the non-paraxial complex-source-point cylindrical wave (CSPCW). High-order CSPCW solutions, termed here high-order quasi-Gaussian cylindrical beams, which exactly satisfy the Helmholtz equation, are derived analytically. Moreover, partial-derivatives of the high-order CSPCW solutions satisfy the Helmholtz equation. In addition, the CSPCW solutions satisfy the nonrelativistic Schrödinger equation within standard quantum mechanics, thus, the results can be used in the description of elementary particle/matter motion and related applications in quantum scattering theory. Furthermore, the analysis is extended to the case of vector beams in which the components of the electromagnetic (EM) field are obtained based on different polarizations of the magnetic and electric vector potentials, which exactly satisfy Maxwell’s vectorial equations and Lorenz’ gauge condition. An attractive feature of the high-order solutions is the rigorous description of strongly focused (or strongly divergent) cylindrical wave-fields without any approximations, nor the need for numerical methods. Possible applications are in beam-forming design using high-aperture or collimated cylindrical laser/electron quasi-Gaussian beams in imaging microscopy, particle manipulation, optical tweezers, and the study of the scattering, and radiation forces on objects. - Highlights: • Bessel beam is transformed into the non-paraxial cylindrical complex-source-point. • Exact high-order tightly focused solutions are derived without any approximations. • The exact solutions also satisfy the nonrelativistic Schrödinger equation. • Electromagnetic beams are obtained as solutions of Maxwell’s vectorial equations. • Applications are in laser/electron beam imaging, tweezers, and radiation force.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗