Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical analysis”

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 919 records · Page 51

The impacts of convex piecewise linear cost formulations on AC optimal power flow

Despite strong connections through shared application areas, research efforts on power market optimization (e.g., unit commitment) and power network optimization (e.g., optimal power flow) remain largely independent. A notable illustration of this is the treatment of power generation cost functions, where nonlinear network optimization has largely used polynomial representations and market optimization has adopted piecewise linear encodings. This work combines state-of-the-art results from both lines of research to understand the best mathematical formulations of the nonlinear AC optimal power flow problem with piecewise linear generation cost functions. An extensive numerical analysis of non-convex models, linear approximations, and convex relaxations across fifty-four realistic test cases illustrates that nonlinear optimization methods are surprisingly sensitive to the mathematical formulation of piecewise linear functions. The results indicate that a poor formulation choice can slow down algorithm performance by a factor of ten, increasing the runtime from seconds to minutes. Furthermore, these results provide valuable insights into the best formulations of nonlinear optimal power flow problems with piecewise linear cost functions, an important step towards building a new generation of energy markets that incorporate the nonlinear AC power flow model.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Performance study of the cross-flow Maisotsenko cycle in humid climate conditions

This paper presents a first numerical analysis of the cross-flow Maisotsenko cycle (M-Cycle) below dew-point evaporative cooler for humid climates with model including the water vapor condensation phenomenon. Proposed solution allows to cool ambient air below its dew-point temperature, by utilizing cooling energy in air extracted from the building, due to that it can be used as effective pre-cooling solution for traditional systems. Proposed unit has significant advantage over traditional counter and cross-flow units operating in heat recovery mode, hence it is almost insensitive on temperature of exhaust air. This is caused by effective pre-cooling of the exhaust air in its initial part. Operation in humid climates using cold extracted air results in condensation occurrence in the primary air channels, and hence partial dehumidification of the primary airflow (which allows to cool the air before its dew- point temperature). The Maisotsenko cycle exchanger used to recover heat in the air-conditioning system is compared with Maisotsenko exchanger that operates only on the ambient airflow and does not allow to recover heat from the exhaust airflow from a conditioned space. Achieved results showed that M-Cycle below dew-point air cooler can operate more effectively in humid climates than traditional cross-flow Maisotsenko cycle exchanger.

Condensation↗

Multi-dimensional modeling of mixture preparation in a direct injection engine fueled with gaseous hydrogen

With the recent advances of direct injection (DI) technology, introducing hydrogen into the combustion chamber through DI is being considered as a viable approach to circumvent backfire and pre-ignition encountered in early generations of hydrogen engines. As part of a broader vision to develop a robust numerical model to study hydrogen spark ignition (SI) combustion in internal combustion (IC) engines, the present numerical investigation focuses on mixture preparation in a hydrogen DI SI engine. This study is carried out with a single hole injector with gaseous hydrogen injected at 100 bar injection pressure. Simulations are carried out for high and low tumble configurations and validated against optical data acquired from planar laser induced fluorescence (PLIF) measurements. Varying mesh configurations are investigated for the impact on in-cylinder mixture distribution. A particular emphasis is placed on the effect of nozzle geometry and mesh orientation near the wall. Overall, the computational model is found to predict the mixture distribution in the combustion cylinder reasonably well. The results showed that the alignment of mesh with the flow direction is important to achieve good agreement between numerical analysis and optical measurement data.

08 HYDROGEN↗

The mimetic finite difference method for elliptic and parabolic problems with a staggered discretization of diffusion coefficient

Numerical schemes for nonlinear parabolic equations based on the harmonic averaging of cell-centered diffusion coefficients break down when some of these coefficients go to zero or their ratio grows. To tackle this problem, we propose new mimetic finite difference schemes that use a staggered discretization of the diffusion coefficient. The primary mimetic operator approximates div (k•); the derived (dual) mimetic operator approximates - ∇(•). The new mimetic schemes preserve symmetry and positive-definiteness of the continuum problem which allows us to use algebraic solvers with optimal complexity. We perform detailed numerical analysis of the new schemes for linear elliptic problems and a specially designed linear parabolic problem that has solution dynamics typical for nonlinear problems. We also show that the new schemes are competitive with the state-of-the-art schemes for steady-state problems but provide much more accurate solution dynamics for the transient problem.

97 MATHEMATICS AND COMPUTING↗

Design, performance, and cost savings of using GF-PC additively manufactured tooling for stamping of HSS 590 sheet metal

Conventional sheet metal forming tooling in the automotive industry is made up of hardened steel and used for mass-production. Prototype tooling made of metal is durable, but it is only used for a small number of parts despite its high cost, contributing heavily to the vehicle development cost. Additive Manufacturing (AM) offers a low-cost alternative for fabricating tooling suitable for low-volume production. Here, this study investigates the feasibility of using AM polymer composite tooling for the stamping of 1.5mm thickness HSS 590 steel sheets through a two-pronged approach – experimental and numerical analysis. Sheet metal stamping experiments were conducted with glass fiber reinforced polycarbonate (GF-PC) AM polymer tooling and their performance was evaluated based on various metrics such as tool deformation and part accuracy. Finite element simulations of the stamping process were performed, which accurately captured the tool deformations seen in experiments with the use of anisotropic material models. The effective use of simulations in optimizing process parameters to achieve the desired final part geometry is also demonstrated. The GF-PC AM polymer composite tooling was found to be well suited for low volume production applications with potential cost savings.

36 MATERIALS SCIENCE↗

The effect of fuel thermal conductivity on the behavior of LWR cores during loss-of-coolant accidents

The effect of variation in thermal conductivity of light water reactor fuel elements on core response during loss-of-coolant accident scenarios is examined. Initially, a simplified numerical analysis is utilized to determine the time scales associated with dissipation of stored energy from the fuel into the coolant once the fission reaction is stopped. The analysis is then followed by full reactor system thermal-hydraulics analysis of a typical boiling and pressurized water reactor subjected to a large break loss-of-coolant accident scenario using the TRACE code. Accordingly, sensitivity analyses to examine the effect of an increase in fuel thermal conductivity, up to 500%, on fuel temperature evolution during these transients are performed. Given the major differences in thermal-hydraulics design aspects of boiling and pressurized water reactors, different fuel and temperature responses during the simulated loss-of-coolant transients are observed.

36 MATERIALS SCIENCE↗

Micromechanical response of SiC-OPyC layers in TRISO fuel particles

Tristructural isotropic (TRISO)–coated particle fuel is a proposed fuel for multiple advanced reactor concepts. The performance of the particle depends on whether the silicon carbide (SiC) layer remains intact to prevent the release of metallic and gaseous fission products. Mechanical fracture of the SiC layer is a potential failure mode under various fuel configurations and operating environments, including the potential transmission of matrix-originating cracks through TRISO particles. Furthermore, this study uses instrumented indentation techniques on cross-sectioned surrogate particles to examine the mechanical stability of the critical interface between SiC and the outer pyrolytic carbon (OPyC) layer. The observed behavior at the interface is rationalized by examining the radially dependent fracture behavior of the SiC layer and performing a numerical analysis to quantify the residual stresses that develop during the processing and cross-sectioning of the as-fabricated particle. Characterizing the SiC-OPyC interface of surrogate TRISO particles using nanoindentation provides unique insight into the interface's room-temperature residual stress and mechanical stability. The modeling efforts were used to investigate the experimental procedure further, and the results are presented herein to validate this fuel form's potential mechanical failure modes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An event-triggered coincidence algorithm for fast-neutron multiplicity assay corrected for cross-talk and photon breakthrough

Here, a model quantifying detector cross-talk and the misidentification of events in fast neutron coincidence distributions is described. This is demonstrated for two experimental arrangements comprising rings of 8 and 15 organic liquid scintillation detectors. Correction terms developed as part of this model are tested with 252 Cf and a relationship is developed between the 235 U enrichment of U 3 O 8 and the order of correlated, fast neutron multiplets induced by an americium–lithium source. The model is also supported by Geant4 simulations. The results suggest that a typical assay, for experimental arrangements that are similar to the examples investigated in this research, will exhibit cross-talk for less than 1% of all detected fast neutrons but, if not accounted for, this can bias the numerical analysis by a margin of 10% and 35% in second- and third-order coincidences (i.e. couplet and triplet counts), respectively. Further, for the case of 252 Cf, it is shown that a relatively low proportion of 4% breakthrough by rays (that is, photons misidentified as neutrons by the pulse-shape discrimination process) can lead to an erroneous increase of 20% in total neutron counts in the assay of a mixed-field, in this case of 252 Cf. These findings will help direct the developments needed to enable organic scintillation detectors with pulse shape discriminators to be applied reliably to nuclear safeguards and non-proliferation verification tasks.

42 ENGINEERING↗

Oscillating surge wave energy converter using a novel above-water power takeoff with belt-arc speed amplification

We investigate the performance of a novel power takeoff (PTO) featuring belt-arc speed amplification for oscillating surge wave energy converters (OSWECs), aiming to address the challenges of extremely low rotary speed and large torque under the low-frequency ocean wave excitations. The belt-arc design significantly increases the rotary speed of the generator, enabling generator downsizing and decreasing the powertrain friction losses. The design also allows for placing the generator above water, eliminating the need for high Ingress Protection ratings for the generator and powertrain, and potentially leading to substantial reductions in capital and maintenance costs. Using the linear potential wave theory, the dynamics of the integrated system are analyzed, and key parameters are identified. To validate the numerical analysis, a 1:10 scale model is designed, fabricated, and tested in a wave tank. Performance evaluations are conducted under regular and irregular wave conditions, with quantified effects of parameter tuning. The results reveal an optimal wave-to-electric efficiency of 48% under regular wave excitation and 20% under irregular excitation. Furthermore, these findings underscore the effectiveness of the proposed novel PTO design in addressing the challenges of low rotation speed and large torque inherent in OSWECs, demonstrating its ability to efficiently convert wave power into electricity.

16 TIDAL AND WAVE POWER↗

Particle multiplicities in the central region of high-energy collisions from $k_T$-factorization with running coupling corrections

Horowitz and Kovchegov have derived a $k_T$-factorization formula for particle production at small x which includes running coupling corrections. We perform a first numerical analysis to confront the theory with data on the energy and centrality dependence of particle multiplicities at midrapidity in high-energy p+A (and A+A) collisions. Moreover, we point out a strikingly different dependence of the multiplicity per participant on $N$ part in p+Pb vs. Pb+Pb collisions at LHC energies, and argue that the observed behavior follows rather naturally from the convolution of the gluon distributions of an asymmetric vs. symmetric projectile and target.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Control-inspired design and power optimization of an active mechanical motion rectifier based power takeoff for wave energy converters

Ocean waves have high energy density and are persistent and predictable. Yet, converting wave energy to a useable form remains challenging. A significant hurdle is the oscillatory nature of waves resulting in the alternating loads, which necessitate the use of rectification at some stage of the energy conversion. This research effort presents a novel design of active mechanical motion rectifier (AMMR) for the power takeoff (PTO), which provides enhanced controllability and better power performance when compared to passive mechanical motion rectifiers (MMR). Inspired by transistors used in synchronous electrical rectifiers, the proposed design uses controllable electromagnetic clutches in the mechanical transmission to allow active engagement-disengagement control; thus, rectifying the oscillatory motion into a unidirectional rotation for high energy conversion efficiency and allowing the generator in unidirectional rotation to control the bidirectional wave capture structure for maximizing the power output. A semi-analytical computational approach is developed to efficiently evaluate the optimal power achieved using the proposed AMMR-based PTO and active control. It is found that the AMMR-based PTO design yields a higher optimal power than the previous passive MMR design across the wave spectrum. The influences of generator inertia and reactive power are discussed. Furthermore, the effects of control parameters on the power output and the optimal trajectories are analyzed. Wave tank tests with the AMMR prototype demonstrated the effectiveness of AMMR based PTO design and validated the numerical analysis.

16 TIDAL AND WAVE POWER↗

Recent advances in computational mathematics and applications

We are honored to bring you this special issue dedicated to recent advances in computational mathematics and applications in science and engineering. The papers in this special issue were presented in ”Conference on Computational Mathematics and Applications (CCMA)” held at the University of Nevada Las Vegas (UNLV) during October 25–27, 2019. The conference has attracted about 90 attendees from six countries. The 15 papers in the special issue comprise a diverse collection of theoretical numerical analysis as well as applications in various subjects. Topics cover optimal control of PDEs, discontinuous Galerkin methods, fluid flow in porous media, turbulence flow, static and moving interface problems, complex fluids composed by the mixture of Newtonian fluid and nematic (liquid crystal) flows, modeling of ocean–atmosphere system, image segmentation algorithm, uncertainty quantification problems for turbulence model and Maxwell’s equations, iterative solver for systems resulting from mixed methods, and novel methods for solving systems of nonlinear equations. As the Guest Editors we would like to thank Dr. Leland Jameson at NSF (National Science Foundation) for kindly supporting our conference proposal which made this conference possible. We would also like to thank Darren Sugrue from Elsevier who provided partial support for our conference, and Thennarasu Gunasekaran and his production team from Elsevier who produced this nice special issue. We are grateful to many people (Dr. Zhijian Wu, Lori Ornelas, Elsa Juarez, and many Ph.D. students) at the Department of Mathematical Sciences of UNLV who provided tremendous support for the conference. Finally, we very much appreciate all authors who contributed their precious results to this special issue.

97 MATHEMATICS AND COMPUTING↗

A supply chain network design model for biomass co-firing in coal-fired power plants

We propose a framework for designing the supply chain network for biomass co-firing in coal-fired power plants. This framework is inspired by existing practices with products with similar physical characteristics to biomass. We present a hub-and-spoke supply chain network design model for long-haul delivery of biomass. This model is a mixed integer linear program solved using benders decomposition algorithm. Numerical analysis indicates that 100 million tons of biomass are located within 75 miles from a coal plant and could be delivered at $8.53/dry-ton; 60 million tons of biomass are located beyond 75 miles and could be delivered at $36/dry-ton.

09 BIOMASS FUELS↗

Convective heat transfer and friction factor characteristics of molten salts in spirally fluted tubes

Spirally fluted tubes have been widely used for heat exchangers due to their superior heat transfer enhancement. However, most of the previous studies focused on the effects of a limited number of geometric parameters, i.e., the flute pitch and flute depth, on convective heat transfer and friction factor characteristics of low-Prandtl-number fluids, i.e., air and water. The correlations developed in these studies may not be accurate or applicable for medium-Prandtl-number fluids, such as molten salts. A numerical analysis using a Computational Fluid Dynamics (CFD) tool, STAR–CCM+, is therefore carried out in this study to systematically investigate the effects of four geometric parameters, including the flute pitch ρ, flute depth e, flute start number N s (or flute helix angle θ), and trough length L tr on convective heat transfer and friction factor characteristics of a medium-Prandtl-number fluid, FLiNaK (46.5LiF-11.5NaF-42KF mol %), in spirally fluted tubes. Additionally, the convective heat transfer and Darcy friction factor correlations are proposed and validated, with ± 20% uncertainties, for medium-Prandtl-number fluids under the following conditions: Re = 88–1600, Pr = 2.5–40, ρ/D c = 0.44–3.51, e/D c = 0.10–0.40, θ/90= 0.20–0.81, and L tr /D c = 0.71–2.16. The correlations proposed help improve the design of spirally fluted-tube heat exchangers.

42 ENGINEERING↗

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Task-oriented machine learning surrogates for tipping points of agent-based models

We present a machine learning framework bridging manifold learning, neural networks, Gaussian processes, and Equation-Free multiscale approach, for the construction of different types of effective reduced order models from detailed agent-based simulators and the systematic multiscale numerical analysis of their emergent dynamics. The specific tasks of interest here include the detection of tipping points, and the uncertainty quantification of rare events near them. Our illustrative examples are an event-driven, stochastic financial market model describing the mimetic behavior of traders, and a compartmental stochastic epidemic model on an Erdös-Rényi network. We contrast the pros and cons of the different types of surrogate models and the effort involved in learning them. Importantly, the proposed framework reveals that, around the tipping points, the emergent dynamics of both benchmark examples can be effectively described by a one-dimensional stochastic differential equation, thus revealing the intrinsic dimensionality of the normal form of the specific type of the tipping point. This allows a significant reduction in the computational cost of the tasks of interest.

97 MATHEMATICS AND COMPUTING↗

X-ray free electron laser heating of water and gold at high static pressure

Probing of reactive materials such as H 2 O ices and fluids at the high pressures and temperatures of planetary interiors is limited by unwanted chemical reactions and confinement failure. Faster experiments can mitigate such issues, but the common approach of adiabatic compression limits the conditions achieved. This study demonstrates a fast experimental strategy for the creation and probing of selected extreme states using static compression coupled with ultrafast X-ray laser heating. Indirect X-ray heating of H 2 O through the use of a gold absorber is evidenced by sample melting inferred from textural changes in the H 2 O diffraction lines and inter-dispersion of gold and H 2 O melts. Coupled with numerical analysis of femtosecond energy absorption, thermal equilibration, and heat transfer, all evidence indicates that temperatures in excess of an electron volt have been reached in the H 2 O at high pressure. Even after repeated heating, samples stayed chemically unchanged from the starting material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Planck 2013 results. XXII. Constraints on inflation

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to be n s = 0.9603 ± 0.0073, ruling out exact scale invariance at over 5σ. Planck establishes an upper bound on the tensor-to-scalar ratio of r< 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V "< 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n = 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining d n s /dlnk =-0.0134 ± 0.0090. Here, we verify these conclusions through a numerical analysis, which makes no slow-roll approximation, and carry out a Bayesian parameter estimation and model-selection analysis for a number of inflationary models including monomial, natural, and hilltop potentials. For each model, we present the Planck constraints on the parameters of the potential and explore several possibilities for the post-inflationary entropy generation epoch, thus obtaining nontrivial data-driven constraints. We also present a direct reconstruction of the observable range of the inflaton potential. Unless a quartic term is allowed in the potential, we find results consistent with second-order slow-roll predictions. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit by Δχ 2 eff &approx; 10 ; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with Planck bounds on f NL . Planck constrains with unprecedented accuracy the amplitude and possible correlation (with the adiabatic mode) of non-decaying isocurvature fluctuations. The fractional primordial contributions of cold dark matter (CDM) isocurvature modes of the types expected in the curvaton and axion scenarios have upper bounds of 0.25% and 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelated isocurvature component can improve the χ 2 eff by approximately 4 as a result of slightly lowering the theoretical prediction for the &ell; ≳ 40 multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗