Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical method”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Thermal Modelling of the in-Situ Consolidation of Automated Fiber Placement of Thermoplastic Composites

NASA is developing carbon fiber reinforced thermoplastic composites processing methods under the Hi-rate Composites Aircraft Manufacturing (HiCAM) project. The in-situ consolidation automated fiber placement (AFP) of thermoplastics (ICAT) process has high potential to increase manufacturing throughput and lower costs as it combines the repeatability and fast lay down rates of AFP with the thermoformability of thermoplastic materials to achieve part consolidation out of the autoclave. Physics-based process models are employed in the ICAT process development to understand, fundamentally, the thermal response of the carbon fiber (CF)/polyaryletherketone (PAEK) material during the rapid heating and cooling associated with laser-assisted AFP. A one-dimensional (1-D) closed-form model was developed capable of predicting the temperature profile through-the-thickness of the tape material. In addition, a two-dimensional (2-D) model was developed to predict the heat transfer through-the-thickness and in the laydown x-direction of the AFP head motion. The solution of the Fourier heat-transfer equations in the 2-D model is approximated using the explicit finite difference numerical method. The temperature profile during placement of PAEK slit-tape materials at various laser target temperatures and placement speeds were measured during ICAT process development trials at Electroimpact®, Inc. †† and the resulting experimental data is compared with the predictions of the 1-D and 2-D models.

Thermoplastics↗
NASA NTRSDate not supplied

Thermal Modelling of the In-Situ Consolidation of Automated Fiber Placement Of Thermoplastic Composites

NASA is developing carbon fiber reinforced thermoplastic composites processing methods under the Hi-rate Composites Aircraft Manufacturing (HiCAM) project. The in-situ consolidation automated fiber placement (AFP) of thermoplastics (ICAT) process has high potential to increase manufacturing throughput and lower costs as it combines the repeatability and fast lay down rates of AFP with the thermoformability of thermoplastic materials to achieve part consolidation out of the autoclave. Physics-based process models are employed in the ICAT process development to understand, fundamentally, the thermal response of the carbon fiber (CF)/polyaryletherketone (PAEK) material during the rapid heating and cooling associated with laser-assisted AFP. A one-dimensional (1-D) closed-form model was developed capable of predicting the temperature profile through-the-thickness of the tape material. In addition, a two-dimensional (2-D) model was developed to predict the heat transfer through-the-thickness and in the laydown x-direction of the AFP head motion. The solution of the Fourier heat-transfer equations in the 2-D model is approximated using the explicit finite difference numerical method. The temperature profile during placement of PAEK slit-tape materials at various laser target temperatures and placement speeds were measured during ICAT process development trials at Electroimpact®, Inc. †† and the resulting experimental data is compared with the predictions of the 1-D and 2-D models.

Christopher J. Stelter↗
DOE OSTI2023

Insights from high-fidelity modeling of industrial rotary bell atomization

The global automotive industry sprayed over 2.6 billion liters of paint in 2018, much of which through electrostatic rotary bell atomization, a highly complex process involving the fluid mechanics of rapidly rotating thin films tearing apart into micrometer-thin filaments and droplets. Coating operations account for 65% of the energy usage in a typical automotive assembly plant, representing 10,000s of gigawatt-hours each year in the United States alone. Optimization of these processes would allow for improved robustness, reduced material waste, increased throughput, and significantly reduced energy usage. Here, we introduce a high-fidelity mathematical and algorithmic framework to analyze rotary bell atomization dynamics at industrially relevant conditions. Our approach couples laboratory experiment with the development of robust non-Newtonian fluid models; devises high-order accurate numerical methods to compute the coupled bell, paint, and gas dynamics; and efficiently exploits high-performance supercomputing architectures. These advances have yielded insight into key dynamics, including i) parametric trends in film, sheeting, and filament characteristics as a function of fluid rheology, delivery rates, and bell speed; ii) the impact of nonuniform film thicknesses on atomization performance; and iii) an understanding of spray composition via primary and secondary atomization. These findings result in coating design principles that are poised to improve energy- and cost-efficiency in a wide array of industrial and manufacturing settings.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

Structure and Stability of One-Dimensional Detonations in Ethylene-Air Mixtures

The propagation of one-dimensional detonations in ethylene-air mixtures is investigated numerically by solving the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing scheme and a point implicit, first-order-accurate, time marching algorithm. The ethylene-air combustion is modeled with a 20-species, 36-step reaction mechanism. A multi-level, dynamically adaptive grid is utilized, in order to resolve the structure of the detonation. Parametric studies over an equivalence ratio range of 0.5 less than phi less than 3 for different initial pressures and degrees of detonation overdrive demonstrate that the detonation is unstable for low degrees of overdrive, but the dynamics of wave propagation varies with fuel-air equivalence ratio. For equivalence ratios less than approximately 1.2 the detonation exhibits a short-period oscillatory mode, characterized by high-frequency, low-amplitude waves. Richer mixtures (phi greater than 1.2) exhibit a low-frequency mode that includes large fluctuations in the detonation wave speed; that is, a galloping propagation mode is established. At high degrees of overdrive, stable detonation wave propagation is obtained. A modified McVey-Toong short-period wave-interaction theory is in excellent agreement with the numerical simulations.

Yungster, S.↗
DOE OSTI2022

Application of the Continuum Damage Mechanics Wilshire-Cano-Stewart (WCS) Model

In this study, the applications of the continuum damage mechanics-based Wilshire-Cano-Stewart (WCS) model are explored to predict rupture time, minimum-creep-strain-rate (MCSR), damage, damage evolution, and creep deformation. Increase knowledge in manufacturing methods has pushed the limit of material science and the development of new materials. Conventional testing is required to qualify materials against creep which according to the ASME B&PV III code, 10,000+ hours of experiments are necessary for each heat before materials are put into service. This process is costly and not feasible for new materials. As an alternative, models have been employed to predict creep behaviors and reduce the amount of time necessary for material qualification. Many models have been developed to predict distinct creep behaviors and the question of which model is best remains. Amongst current models, the WCS model has emerge with the ability to predict multiple behaviors using an explicit analytical approach with the ability to predict long-term creep. In this study the novel continuum damage mechanics WCS is employed in multiple applications. The goals of the study are (a) to discuss and determine the framework of the WCS model and validated it mathematically and using parametric simulations, (b) applied the model to accelerated creep data to show the capabilities of the model with non-conventional data, and (c) applied a novel numerical method, the datum temperature method (DTM) to show the model extrapolations and interpolations capabilities with limited and reduce data sets. To accomplish these goals, data is gathered for alloy P91 and Inconel 718 to develop and post-audit validate the model. The benefits of using the WCS model is that it provides an explicit stress and temperature dependency ideal for extrapolations, the coupled equations are suitable for finite element analysis (FEA) implementation, and it follows an explicit calibration approach. The model also proves that it can be applied to accelerated testing data and using the DTM.

Cano, Jaime A↗
NASA NTRS2017

A New Formulation of Time Domain Boundary Integral Equation for Acoustic Wave Scattering in the Presence of a Uniform Mean Flow

It has been well-known that under the assumption of a constant uniform mean flow, the acoustic wave propagation equation can be formulated as a boundary integral equation, in both the time domain and the frequency domain. Compared with solving partial differential equations, numerical methods based on the boundary integral equation have the advantage of a reduced spatial dimension and, hence, requiring only a surface mesh. However, the constant uniform mean flow assumption, while convenient for formulating the integral equation, does not satisfy the solid wall boundary condition wherever the body surface is not aligned with the uniform mean flow. In this paper, we argue that the proper boundary condition for the acoustic wave should not have its normal velocity be zero everywhere on the solid surfaces, as has been applied in the literature. A careful study of the acoustic energy conservation equation is presented that shows such a boundary condition in fact leads to erroneous source or sink points on solid surfaces not aligned with the mean flow. A new solid wall boundary condition is proposed that conserves the acoustic energy and a new time domain boundary integral equation is derived. In addition to conserving the acoustic energy, another significant advantage of the new equation is that it is considerably simpler than previous formulations. In particular, tangential derivatives of the solution on the solid surfaces are no longer needed in the new formulation, which greatly simplifies numerical implementation. Furthermore, stabilization of the new integral equation by Burton-Miller type reformulation is presented. The stability of the new formulation is studied theoretically as well as numerically by an eigenvalue analysis. Numerical solutions are also presented that demonstrate the stability of the new formulation.

Hu, Fang↗
DOE OSTI2024

Quantifying local and global mass balance errors in physics-informed neural networks

Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Hamiltonian Truncation Effective Theory

Hamiltonian truncation is a non-perturbative numerical method for calculating observables of a quantum field theory. The starting point for this method is to truncate the interacting Hamiltonian to a finite-dimensional space of states spanned by the eigenvectors of the free Hamiltonian H_0 H 0 with eigenvalues below some energy cutoff E_\text{max} E max . In this work, we show how to treat Hamiltonian truncation systematically using effective field theory methodology. We define the finite-dimensional effective Hamiltonian by integrating out the states above E_\text{max} E max . The effective Hamiltonian can be computed by matching a transition amplitude to the full theory, and gives corrections order by order as an expansion in powers of 1/E_\text{max} 1 / E max . The effective Hamiltonian is non-local, with the non-locality controlled in an expansion in powers of H_0/E_\text{max} H 0 / E max . The effective Hamiltonian is also non-Hermitian, and we discuss whether this is a necessary feature or an artifact of our definition. We apply our formalism to 2D \lambda\phi^4 λ ϕ 4 theory, and compute the the leading 1/E_\text{max}^2 1 / E max 2 corrections to the effective Hamiltonian. We show that these corrections nontrivially satisfy the crucial property of separation of scales. Numerical diagonalization of the effective Hamiltonian gives residual errors of order 1/E_\text{max}^3 1 / E max 3 , as expected by our power counting. We also present the power counting for 3D \lambda \phi^4 λ ϕ 4 theory and perform calculations that demonstrate the separation of scales in this theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2016

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.

Majumdar, A. K.↗
DOE OSTI2024

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimization of flow in additively manufactured porous columns with graded permeability

Chemical engineering systems often involve a functional porous medium, such as in catalyzed reactive flows, fluid purifiers, and chromatographic separations. Ideally, the flow rates throughout the porous medium are uniform, and all portions of the medium contribute efficiently to its function. The permeability is a property of a porous medium that depends on pore geometry and relates flow rate to pressure drop. Additive manufacturing techniques raise the possibilities that permeability can be arbitrarily specified in three dimensions, and that a broader range of permeabilities can be achieved than by traditional manufacturing methods. Using numerical optimization methods, we show that designs with spatially varying permeability can achieve greater flow uniformity than designs with uniform permeability. We consider geometries involving hemispherical regions that distribute flow, as in many glass chromatography columns. By several measures, significant improvements in flow uniformity can be obtained by modifying permeability only near the inlet and outlet.

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

Lagrangian–Eulerian multidensity topology optimization with the material point method

Abstract In this paper, a hybrid Lagrangian–Eulerian topology optimization (LETO) method is proposed to solve the elastic force equilibrium with the Material Point Method (MPM). LETO transfers density information from freely movable Lagrangian carrier particles to a fixed set of Eulerian quadrature points. This transfer is based on a smooth radial kernel involved in the compliance objective to avoid the artificial checkerboard pattern. The quadrature points act as MPM particles embedded in a lower‐resolution grid and enable a subcell multidensity resolution of intricate structures with a reduced computational cost. A quadrature‐level connectivity graph‐based method is adopted to avoid the artificial checkerboard issues commonly existing in multiresolution topology optimization methods. Numerical experiments are provided to demonstrate the efficacy of the proposed approach.

Li, Yue↗
DOE OSTI2025

AMR-Wind: A Performance-Portable, High-Fidelity Flow Solver for Wind Farm Simulations

We present AMR-Wind, a verified and validated high-fidelity computational-fluid-dynamics code for wind farm flows. AMR-Wind is a block-structured, adaptive-mesh, incompressible-flow solver that enables predictive simulations of the atmospheric boundary layer and wind plants. It is a highly scalable code designed for parallel high-performance computing with a specific focus on performance portability for current and future computing architectures, including graphical processing units (GPUs). In this paper, we detail the governing equations, the numerical methods, and the turbine models. Establishing a foundation for the correctness of the code, we present the results of formal verification and validation. The verification studies, which include a novel actuator line test case, indicate that AMR-Wind is spatially and temporally second-order accurate. The validation studies demonstrate that the key physics capabilities implemented in the code, including actuator disk models, actuator line models, turbulence models, and large eddy simulation (LES) models for atmospheric boundary layers, perform well in comparison to reference data from established computational tools and theory. We conclude with a demonstration simulation of a 12-turbine wind farm operating in a turbulent atmospheric boundary layer, detailing computational performance and realistic wake interactions.

17 WIND ENERGY↗
DOE OSTI2020

Notes on the complex Sachdev-Ye-Kitaev model

We describe numerous properties of the Sachdev-Ye-Kitaev model for complex fermions with N $\gg$ 1 flavors and a global U(1) charge. We provide a general definition of the charge in the ( G, Σ) formalism, and compute its universal relation to the infrared asymmetry of the Green function. The same relation is obtained by a renormalization theory. The conserved charge contributes a compact scalar field to the effective action, from which we derive the many-body density of states and extract the charge compressibility. We compute the latter via three distinct numerical methods and obtain consistent results. Finally, we present a two dimensional bulk picture with free Dirac fermions for the zero temperature entropy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Calabi-Yau CFTs and random matrices

Using numerical methods for finding Ricci-flat metrics, we explore the spectrum of local operators in two-dimensional conformal field theories defined by sigma models on Calabi-Yau targets at large volume. Focusing on the examples of K3 and the quintic, we show that the spectrum, averaged over a region in complex structure moduli space, possesses the same statistical properties as the Gaussian orthogonal ensemble of random matrix theory.

superstring vacua↗
DOE OSTI2024

BPS fivebrane stars. Part II. Fluctuations

We investigate quantum fluctuations of metric components in coherent 1/2-BPS bound states of n1 fundamental strings and n 5 NS5-branes. The leading order contribution in an expansion in 1/(n 1 n 5 ) is calculated via a combination of analytical and numerical methods. We find that the fluctuations are small away from a tiny distance from the source, comparable to the 6d Planck scale. Comparing this result with an analysis in the literature of fluctuations in the maximally mixed state, we conclude that the large fluctuations previously found for the latter are statistical rather than quantum in nature, and that perturbative string theory provides an accurate description of these backgrounds.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Continuum-mediated self-interacting dark matter

Dark matter may self-interact through a continuum of low-mass states. This happens if dark matter couples to a strongly-coupled nearly-conformal hidden sector. This type of theory is holographically described by brane-localized dark matter interacting with bulk fields in a slice of 5D anti-de Sitter space. The long-range potential in this scenario depends on a non-integer power of the spatial separation, in contrast to the Yukawa potential generated by the exchange of a single 4D mediator. The resulting self-interaction cross section scales like a non-integer power of velocity. We identify the Born, classical and resonant regimes and investigate them using state-of-the-art numerical methods. We demonstrate the viability of our continuum-mediated framework to address the astrophysical small-scale structure anomalies. Investigating the continuum-mediated Sommerfeld enhancement, we demonstrate that a pattern of resonances can occur depending on the non-integer power. We conclude that continuum mediators introduce novel power-law scalings which open new possibilities for dark matter self-interaction phenomenology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗