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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Multi-scale calculation based on dual domain material point method combined with molecular dynamics

This dissertation combines the dual domain material point method (DDMP) with molecular dynamics (MD) in an attempt to create a multi-scale numerical method to simulate materials undergoing large deformations with high strain rates. In these types of problems, the material is often in a thermodynamically non-equilibrium state, and conventional constitutive relations are often not available. In this method, the closure quantities, such as stress, at each material point are calculated from a MD simulation of a group of atoms surrounding the material point. Rather than restricting the multi-scale simulation in a small spatial region, such as phase interfaces, or crack tips, this multi-scale method can be used to consider non-equilibrium thermodynamic e ects in a macroscopic domain. This method takes advantage that the material points only communicate with mesh nodes, not among themselves; therefore MD simulations for material points can be performed independently in parallel. First, using a one-dimensional shock problem as an example, the numerical properties of the original material point method (MPM), the generalized interpolation material point (GIMP) method, the convected particle domain interpolation (CPDI) method, and the DDMP method are investigated. Among these methods, only the DDMP method converges as the number of particles increases, but the large number of particles needed for convergence makes the method very expensive especially in our multi-scale method where we calculate stress in each material point using MD simulation. To improve DDMP, the sub-point method is introduced in this dissertation, which provides high quality numerical solutions with a very small number of particles. The multi-scale method based on DDMP with sub-points is successfully implemented for a one dimensional problem of shock wave propagation in a cerium crystal. The MD simulation to calculate stress in each material point is performed in GPU using CUDA to accelerate the computation. The numerical properties of the multiscale method are investigated as well as the results from this multi-scale calculation are compared of particles needed for convergence makes the method very expensive especially in our multi-scale method where we calculate stress in each material point using MD simulation. To improve DDMP, the sub-point method is introduced in this dissertation, which provides high quality numerical solutions with a very small number of particles. The multi-scale method based on DDMP with sub-points is successfully implemented for a one dimensional problem of shock wave propagation in a cerium crystal. The MD simulation to calculate stress in each material point is performed in GPU using CUDA to accelerate the computation. The numerical properties of the multiscale method are investigated as well as the results from this multi-scale calculation are compared with direct MD simulation results to demonstrate the feasibility of the method. Also, the multi-scale method is applied for a two dimensional problem of jet formation around copper notch under a strong impact.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2015

Pulsed laser generation of ultrasound in a metal plate between the melting and ablation thresholds

The generation of ultrasound in a metal plate exposed to nanosecond pulsed laser heating, sufficient to cause melting but not ablation, is treated. Consideration is given to the spatial and temporal profiles of the laser pulse, penetration of the laser beam into the sample, the evolution of the melt pool, and thermal conduction in the melt and surrounding solid. The excitation of the ultrasound takes place over a few nanoseconds, and occurs predominantly within the thermal diffusion length of a micron or so beneath the surface. Because of this, the output of the thermal simulations can be represented as axially symmetric transient radial and normal surface force distributions. The epicentral displacement response at the opposite surface to these forces is obtained by two methods, the one based on the elastodynamic Green’s functions for plate geometry determined by the Cagniard generalized ray method, and the other using a finite element numerical method. The two approaches are in very close agreement. Numerical simulations are reported of the epicentral displacement response of a 3.12mm thick tungsten plate irradiated with a 4 ns pulsed laser beam with Gaussian spatial profile, at intensities below and above the melt threshold. Comparison is made between results obtained using available temperature dependent thermophysical data, and room temperature materials constants except near the melting point.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography↗

DOE OSTI2014

Dynamics of multiple viscoelastic carbon nanotube based nanocomposites with axial magnetic field

Nanocomposites and magnetic field effects on nanostructures have received great attention in recent years. A large amount of research work was focused on developing the proper theoretical framework for describing many physical effects appearing in structures on nanoscale level. Great step in this direction was successful application of nonlocal continuum field theory of Eringen. In the present paper, the free transverse vibration analysis is carried out for the system composed of multiple single walled carbon nanotubes (MSWCNT) embedded in a polymer matrix and under the influence of an axial magnetic field. Equivalent nonlocal model of MSWCNT is adopted as viscoelastically coupled multi-nanobeam system (MNBS) under the influence of longitudinal magnetic field. Governing equations of motion are derived using the Newton second low and nonlocal Rayleigh beam theory, which take into account small-scale effects, the effect of nanobeam angular acceleration, internal damping and Maxwell relation. Explicit expressions for complex natural frequency are derived based on the method of separation of variables and trigonometric method for the “Clamped-Chain” system. In addition, an analytical method is proposed in order to obtain asymptotic damped natural frequency and the critical damping ratio, which are independent of boundary conditions and a number of nanobeams in MNBS. The validity of obtained results is confirmed by comparing the results obtained for complex frequencies via trigonometric method with the results obtained by using numerical methods. The influence of the longitudinal magnetic field on the free vibration response of viscoelastically coupled MNBS is discussed in detail. In addition, numerical results are presented to point out the effects of the nonlocal parameter, internal damping, and parameters of viscoelastic medium on complex natural frequencies of the system. The results demonstrate the efficiency of the suggested methodology to find the closed form solutions for the free vibration response of multiple nanostructure systems under the influence of magnetic field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1992

Calculation of fully three-dimensional separated flow with an unsteady viscous-inviscid interaction method

Three previous papers have shown that the viscous inviscid numerical methods were capable of calculating separated turbulent flows. The 'Semi-inverse method' and the models of the first author allow the computation of massive separations, stalled flows, and shock wave-boundary layer interactions, in two-dimensional or quasi-three-dimensional conditions, from low speeds to supersonic speeds. The 'Semi-implicit' method for time-consistent coupling allows for the computation of time-accurate transonic separated flow and buffer computations in two-dimensions. The present paper shows that the viscous-inviscid interaction approach is also able to compute the fully three-dimensional flow separation phenomena. The method is based on a thin-layer approximation of the theory of 'Defect-Formulation' that provides the viscous-inviscid splitting of the Navier-Stokes equations. A parametric analytical modelling of the 3D-turbulent velocity profiles is involved. Numerically, the 3D-velocity profiles are discretized in the normal z-direction and driven by parametric integral 3D-equations in direct or inverse modes in the x-direction. The viscous-inviscid coupling is fully 3D and solved the time-consistency problem with an extension of the 'Semi-implicit' method previously suggested in two-dimensions. A 3D inviscid subroutine with TSP approximation is used. Results are obtained for transonic steady flows over wings with shock-induced transonic separation. The method provides results for highly three-dimensional flow separations, such as induced by a three-dimensional through at the wall of a flat plate. The 3D viscous-inviscid coupling and the 3D model of the velocity field provide three-dimensional instantaneous skin-friction lines whose patterns exhibit the same complex topology as Navier-Stokes solvers, with foci, nodes, and saddle-points.

Leballeur, J. C.↗

DOE OSTI2020

Parallel three-dimensional simulations of quasi-static elastoplastic solids

Hypo-elastoplasticity is a flexible framework for modeling the mechanics of many hard materials under small elastic deformation and large plastic deformation. Under typical loading rates, most laboratory tests of these materials happen in the quasi-static limit, but there are few existing numerical methods tailor-made for this physical regime. Here, we extend to three dimensions a recent projection method for simulating quasi-static hypo-elastoplastic materials. The method is based on a mathematical correspondence to the incompressible Navier–Stokes equations, where the projection method of Chorin (1968) is an established numerical technique. We develop and utilize a three-dimensional parallel geometric multigrid solver employed to solve a linear system for the quasi-static projection. Our method is tested through simulation of three-dimensional shear band nucleation and growth, a precursor to failure in many materials. As an example system, we employ a physical model of a bulk metallic glass based on the shear transformation zone theory, but the method can be applied to any elastoplasticity model. We consider several examples of three-dimensional shear banding, and examine shear band formation in physically realistic materials with heterogeneous initial conditions under both simple shear deformation and boundary conditions inspired by friction welding.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

Exponential time differencing methods with Chebyshev collocation for polymers confined by interacting surfaces

We present a fast and accurate numerical method for the self-consistent field theory calculations of confined polymer systems. It introduces an exponential time differencing method (ETDRK4) based on Chebyshev collocation, which exhibits fourth-order accuracy in temporal domain and spectral accuracy in spatial domain, to solve the modified diffusion equations. Similar to the approach proposed by Hur et al. [Macromolecules 45, 2905 (2012)], non-periodic boundary conditions are adopted to model the confining walls with or without preferential interactions with polymer species, avoiding the use of surface field terms and the mask technique in a conventional approach. The performance of ETDRK4 is examined in comparison with the operator splitting methods with either Fourier collocation or Chebyshev collocation. Numerical experiments show that our exponential time differencing method is more efficient than the operator splitting methods in high accuracy calculations. This method has been applied to diblock copolymers confined by two parallel flat surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

An Immersed Boundary method with divergence-free velocity interpolation and force spreading

The Immersed Boundary (IB) method is a mathematical framework for constructing robust numerical methods to study fluid–structure interaction in problems involving an elastic structure immersed in a viscous fluid. The IB formulation uses an Eulerian representation of the fluid and a Lagrangian representation of the structure. The Lagrangian and Eulerian frames are coupled by integral transforms with delta function kernels. The discretized IB equations use approximations to these transforms with regularized delta function kernels to interpolate the fluid velocity to the structure, and to spread structural forces to the fluid. It is well-known that the conventional IB method can suffer from poor volume conservation since the interpolated Lagrangian velocity field is not generally divergence-free, and so this can cause spurious volume changes. In practice, the lack of volume conservation is especially pronounced for cases where there are large pressure differences across thin structural boundaries. The aim of this paper is to greatly reduce the volume error of the IB method by introducing velocity-interpolation and force-spreading schemes with the properties that the interpolated velocity field in which the structure moves is at least $\mathscr{C}$ 1 and satisfies a continuous divergence-free condition, and that the force-spreading operator is the adjoint of the velocity-interpolation operator. We confirm through numerical experiments in two and three spatial dimensions that this new IB method is able to achieve substantial improvement in volume conservation compared to other existing IB methods, at the expense of a modest increase in the computational cost. Further, the new method provides smoother Lagrangian forces (tractions) than traditional IB methods. Finally, the method presented here is restricted to periodic computational domains. Its generalization to non-periodic domains is important future work.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2004

Sensitivity of Simulated Warm Rain Formation to Collision and Coalescence Efficiencies, Breakup, and Turbulence: Comparison of Two Bin-Resolved Numerical Models

Numerical models that resolve cloud particles into discrete mass size distributions on an Eulerian grid provide a uniquely powerful means of studying the closely coupled interaction of aerosols, cloud microphysics, and transport that determine cloud properties and evolution. However, such models require many experimentally derived paramaterizations in order to properly represent the complex interactions of droplets within turbulent flow. Many of these parameterizations remain poorly quantified, and the numerical methods of solving the equations for temporal evolution of the mass size distribution can also vary considerably in terms of efficiency and accuracy. In this work, we compare results from two size-resolved microphysics models that employ various widely-used parameterizations and numerical solution methods for several aspects of stochastic collection.

Fridlind, Ann↗

DOE OSTI2020

Meshfree methods on manifolds for hydrodynamic flows on curved surfaces: A Generalized Moving Least-Squares (GMLS) approach

We utilize generalized moving least squares (GMLS) to develop meshfree techniques for discretizing hydrodynamic flow problems on manifolds. We use exterior calculus to formulate incompressible hydrodynamic equations in the Stokesian regime and handle the divergence-free constraints via a generalized vector potential. This provides less coordinate-centric descriptions and enables the development of efficient numerical methods and splitting schemes for the fourth-order governing equations in terms of a system of second-order elliptic operators. Using a Hodge decomposition, we develop methods for manifolds having spherical topology. We show the methods exhibit high-order convergence rates for solving hydrodynamic flows on curved surfaces. The methods also provide general high-order approximations for the metric, curvature, and other geometric quantities of the manifold and associated exterior calculus operators. The approaches also can be utilized to develop high-order solvers for other scalar-valued and vector-valued problems on manifolds.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1974

Analysis of twisting stiffness for a multifiber composite layer

The twisting stiffness of a rectangular cross section consisting of a single row of solid circular cross-section fibers embedded in a matrix is analyzed. The problem is formulated as a Dirichlet torsion problem of a multielement region and solved by the boundary-point least-squares method. Numerical results for a single-fiber square cross section compare favorably with previous relaxation-method results. New numerical results for three and five-fiber composites suggest that the torsional rigidity of a multifiber composite can be approximated from the torsional rigidities of single and three-fiber models.

Bert, C. W.↗

NASA NTRS1989

LU-SGS implicit algorithm for three-dimensional incompressible Navier-Stokes equations with source term

A numerical method is developed for solving the incompressible Navier-Stokes equations using the concept of pseudocompressibility. A lower-upper symmetric-Gauss-Seidel implicit scheme is developed for three-dimensional incompressible viscous flow computations. The present algorithm offers additional advantages when solving the flow equations with source terms. Complete vectorizability of the algorithm on oblique planes of sweep in three-dimensions is accomplished in a new flow solver, INS3D-LU code. Spatial differencing is a second-order accurate semi-discrete finite-volume method augmented by a third-order accurate numerical dissipation model which is based on spectral-radii. Comparison of numerical solutions for a curved duct with experimental data shows good agreement. The method is applied to calculate the inducer flow of the Space Shuttle Main Engine turbopump.

Yoon, Seokkwan↗

DOE OSTI2013

The coordinate transformation method for design of polarizers on HL-2A electron cyclotron resonance heating and current drive systems

Polarizers are widely used to change the polarization of millimeter waves on the electron cyclotron resonance heating and current drive (ECRH and CD) systems. A new method based on the coordinate transformation and the Fourier expansion (the so-called C-method) has been developed for design of polarizers on the HL-2A ECRH and CD systems. This method transforms the grating problem to an eigenvalue problem, making it easy and clear to understand and solve. The comparison between the C-method, the integral method, and the low power test results is presented. It indicates that the C-method can be considered as a rigorous numerical method for the design of polarizers. Finally, two polarizers were designed based on the C-method which can be used together to achieve almost arbitrary polarization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Mutual information-assisted adaptive variational quantum eigensolver

Adaptive construction of ansatz circuits offers a promising route towards applicable variational quantum eigensolvers on near-term quantum hardware. Those algorithms aim to build up optimal circuits for a certain problem and ansatz circuits are adaptively constructed by selecting and adding entanglers from a predefined pool. In this work, we propose a way to construct entangler pools with reduced size by leveraging classical algorithms. Our method uses mutual information between the qubits in classically approximated ground state to rank and screen the entanglers. The density matrix renormalization group method is employed for classical precomputation in this work. We corroborate our method numerically on small molecules. Our numerical experiments show that a reduced entangler pool with a small portion of the original entangler pool can achieve same numerical accuracy. Here, we believe that our method paves a new way for adaptive construction of ansatz circuits for variational quantum algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Self-learning quantum Monte Carlo method in interacting fermion systems

We present the self-learning Monte Carlo method is a powerful general-purpose numerical method recently introduced to simulate many-body systems. In this work, we extend it to an interacting fermion quantum system in the framework of the widely used determinant quantum Monte Carlo. This method can generally reduce the computational complexity and moreover can greatly suppress the autocorrelation time near a critical point. This enables us to simulate an interacting fermion system on a $100 × 100$ lattice even at the critical point and obtain critical exponents with high precision.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Equation of State and Thermometry of the 2D SU ( N ) Fermi-Hubbard Model

We characterize the equation of state (EoS) of the SU ( N > 2 ) Fermi-Hubbard Model (FHM) in a two-dimensional single-layer square optical lattice. We probe the density and the site occupation probabilities as functions of interaction strength and temperature for N = 3 , 4, and 6. Our measurements are used as a benchmark for state-of-the-art numerical methods including determinantal quantum Monte Carlo and numerical linked cluster expansion. By probing the density fluctuations, we compare temperatures determined in a model-independent way by fitting measurements to numerically calculated EoS results, making this a particularly interesting new step in the exploration and characterization of the SU ( N ) FHM. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Method of Moving Asymptotes Optimization Algorithm in Python

Python implementation of the Method of Moving Asymptotes optimization algorithm described in[Svanberg, K., The method of moving asymptotes- a new method for structural optimization. International journal for numerical methods in engineering, 1987. 24(2): p.359 (https://onlinelibrary.wiley.com/doi/abs/10.1002/nme.1620240207). Originally implemented in[GetDP](https://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=717799) [project](https://gitlab.onelab.info/getdp/getdp), although it has been deleted in their repository.

Salazar De Troya, Miguel↗