Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Dynamic Performance Enhancement for Nonlinear Stochastic Systems Using RBF Driven Nonlinear Compensation with Extended Kalman Filter

In this paper, a novel hybrid control method is proposed to enhance the control performance of the PI based control system for a class of nonlinear and non-Gaussian stochastic dynamic processes with unmeasurable states. Firstly, to enhance the tracking performance of the PI controller where the PI parameters are fixed in many actual control systems, the compensative signal is formed using the extended Kalman filter(EKF) based state estimator and driven by the radial basis function (RBF) neural network based compensator. In addition, the weights of RBF is trained to minimize the entropy criterion of tracking error as the process is subjected to non-Gaussian disturbances. Meanwhile, since the precise statistical property of noises is hard to obtain for many industrial processes, the kernel density estimation (KDE) technique is employed in this paper to estimate the entropy. The convergence of RBF network is discussed and the stability of the resulted closed-loop hybrid control system is analyzed in mean square sense. Finally, a numerical example and a practical system testing are given to illustrate the effectiveness of the proposed control method.

stochastic systems, PI Controller, Entropy, Extend↗
DOE OSTI2020

On differentiable local bounds preserving stabilization for Euler equations

This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. Furthermore, the importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. Additionally, in the case of transient problem, we also observe a reduction in the computational cost.

42 ENGINEERING↗
DOE OSTI2021

How Equivalent Are Equivalent Porous Media?

Abstract Geoenergy and geoengineering applications usually involve fluid injection into and production from fractured media. Accounting for fractures is important because of the strong poromechanical coupling that ties pore pressure changes and deformation. A possible approach to the problem uses equivalent porous media to reduce the computational cost and model complexity instead of explicitly including fractures in the models. We investigate the validity of this simplification by comparing these two approaches. Simulation results show that pore pressure distribution significantly differs between the two approaches even when both are calibrated to predict identical values at the injection and production wells. Additionally, changes in fracture stability are not well captured with the equivalent porous medium. We conclude that explicitly accounting for fractures in numerical models may be necessary under some circumstances to perform reliable coupled thermohydromechanical simulations, which could be used in conjunction with other tools for induced seismicity forecasting.

58 GEOSCIENCES↗
DOE OSTI2024

A comprehensive review on liquid electrolyte design for low-temperature lithium/sodium metal batteries

Lithium/sodium metal batteries (LMBs/SMBs) possess immense potential for various applications due to their high energy density. Nevertheless, LMBs/SMBs are highly susceptible to the detrimental effects of an unstable solid electrolyte interphase (SEI) and dendrites during practical applications, particularly pronounced in low-temperature environments. Furthermore, sluggish ion transportation further compromises the cycling stability of LMBs/SMBs at low temperatures. To achieve stable operation of LMBs/SMBs at low temperatures, researchers have made numerous efforts including electrolyte optimization aimed at creating stable SEIs and suppressing the metal dendrites under low temperature conditions. Despite the significant advancements made recently in the liquid electrolyte design, there remains considerable hurdle in electrolyte engineering for practical low-temperature, high energy density LMBs/SMBs, calling for a profound comprehension of the intricate interplay between the electrochemical reaction kinetics and electrolyte compositions. Here, this review provides a thorough overview of various strategies in optimizing liquid electrolytes covering weakly solvating electrolytes, concentration-designed electrolytes, and solvation structure-designed electrolytes, to address the challenges faced by LMBs/SMBs at low temperatures, including slow reaction kinetics and the difficulties in Li + /Na + solvation/desolvation. Furthermore, this review discusses future prospects for the advancement of this field, intending to provide valuable insights and support for subsequent research undertakings.

25 ENERGY STORAGE↗
DOE OSTI2025

FreeMHD: Validation and verification of the open-source, multi-domain, multi-phase solver for electrically conductive flows

The extreme heat fluxes in the divertor region of tokamaks may require an alternative to solid plasma-facing components, for the extraction of heat and the protection of the surrounding walls. Flowing liquid metals are proposed as an alternative, but raise additional challenges that require investigation and numerical simulations. Free surface designs are desirable for plasma-facing components, but steady flow profiles and surface stability must be ensured to limit undesirable interactions with the plasma. Previous studies have mainly used steady-state, 2D, or simplified models for internal flows and have not been able to adequately model free-surface liquid metal (LM) experiments. Therefore, FreeMHD has been recently developed as an open-source magnetohydrodynamics (MHD) solver for free-surface electrically conductive flows subject to a strong external magnetic field. The FreeMHD solver computes incompressible free-surface flows with multi-region coupling for the investigation of MHD phenomena involving fluid and solid domains. The model utilizes the finite-volume OpenFOAM framework under the low magnetic Reynolds number approximation. FreeMHD is validated using analytical solutions for the velocity profiles of closed channel flows with various Hartmann numbers and wall conductance ratios. Next, experimental measurements are then used to verify FreeMHD, through a series of cases involving dam breaking, 3D magnetic fields, and free-surface LM flows. These results demonstrate that FreeMHD is a reliable tool for the design of LM systems under free surface conditions at the reactor scale. Furthermore, it is flexible, computationally inexpensive, and can be used to solve fully 3D transient MHD flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Uniform-density Bose-Einstein condensates of the Gross-Pitaevskii equation found by solving the inverse problem for the confining potential

Here, in this work, we consider a “reverse-engineering” approach to construct confining potentials that support exact, constant density kovaton solutions to the classical Gross-Pitaevskii equation (GPE) also known as the nonlinear Schr¨odinger equation (NLSE). In the one-dimensional case, the exact solution is the sum of stationary kink and anti-kink solutions, i.e. a kovaton, and in the overlapping region, the density is constant. In higher dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed out finite square well with minima also at the edges. When self-interactions are added, a term proportional to ±gψ*ψ gets added to the confining potential and ±gM, where M is the norm, gets added to the total energy. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of M (related to the number of particles in the BEC). Comparing the stability criteria from Derrick’s theorem and Bogoliubov-de Gennes analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass M for attractive interactions. However, the numerical analysis gives a much lower critical mass. The numerical analysis shows further that the initial instabilities violate the symmetry x → -x assumed by Derrick’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Vorticity-induced flow-focusing leads to bubble entrainment in an inkjet printhead: Synchrotron x-ray and volume-of-fluid visualizations

The oscillatory flows present in an inkjet printhead can lead to strong deformations of the air-liquid interface at the nozzle exit. Such deformations may lead to an inward directed air jet with bubble pinch-off and the subsequent entrainment of an air bubble, which is highly detrimental to the stability of inkjet printing. Understanding the mechanisms of bubble entrainment is therefore crucial to improving print stability. Here, in the present work, we use ultrafast x-ray phase-contrast imaging and direct numerical simulations based on the volume-of-fluid method to study the mechanisms underlying the bubble entrainment in a piezoacoustic printhead. We first demonstrate good agreement between experiments and numerics. We then show the different classes of bubble pinch-off obtained in experiments, and that those were also captured numerically. The numerical results are then used to show that the baroclinic torque, which is generated at the gas-liquid interface due to the misalignment of density and pressure gradients, results in a flow-focusing effect that drives the formation of the air jet from which a bubble can pinch off.

42 ENGINEERING↗
DOE OSTI2021

Data-Driven Adaptive Damping Controller for Wind Power Plants with Doubly-Fed Induction Generators

This paper presents an adaptive damping controller for wind power plants in which the turbines are equipped with doubly-fed induction generators. The controller is designed to respond to an input control signal that is triggered according to the system operating conditions. A processing unit continuously estimates the electromechanical modes of oscillation based on real-time streaming data acquired from a phasor measurement unit that is strategically positioned on the grid. The decision to trigger (or not trigger) the control signal is automatic, based on the relative damping of the dominant mode. The modes are estimated using the dynamic mode decomposition algorithm with time-delay embedding. Numerical simulations performed on the two-area system demonstrate that the proposed controller enhances the rotor angle stability for both small-signal and large disturbances, and is adaptive to changing grid conditions.

damping controller↗
DOE OSTI2021

Beam Spill in Downstream Transport due to Injector Faults [PowerPoint]

We have investigated injector faults with numerical experiments. For this investigation we used the nominal tune that is also being used to assess beam stability. Beam spill in the DST due to a serious injector fault was evaluated using the LSP slice algorithm. Only the most serious fault investigated (-350 keV) caused beam spill in the DST. DST transport focusing solenoids limit the beam spill. Beam spill in the DST is the result of large halo caused by betatron oscillations of the mismatched beam. Bremsstrahlung radiation from the DST beam spill was estimated from Tom Martin scaling normalized to DARHT-I. In conclusion, Scorpius is reasonably safe from beam spill radiation in the DST due to injector faults.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Land-based wind turbines with flexible rail-transportable blades – Part 2: 3D finite element design optimization of the rotor blades

Abstract. Increasing growth in land-based wind turbine blades to enable higher machine capacities and capacity factors is creating challenges in design, manufacturing, logistics, and operation. Enabling further blade growth will require technology innovation. An emerging solution to overcome logistics constraints is to segment the blades spanwise and chordwise, which is effective, but the additional field-assembled joints result in added mass and loads, as well as increased reliability concerns in operation. An alternative to this methodology is to design slender flexible blades that can be shipped on rail lines by flexing during transport. However, the increased flexibility is challenging to accommodate with a typical glass-fiber, upwind design. In a two-part paper series, several design options are evaluated to enable slender flexible blades: downwind machines, optimized carbon fiber, and active aerodynamic controls. Part 1 presents the system-level optimization of the rotor variants as compared to conventional and segmented baselines, with a low-fidelity representation of the blades. The present work, Part 2, supplements the system-level optimization in Part 1 with high-fidelity blade structural optimization to ensure that the designs are at feasible optima with respect to material strength and fatigue limits, as well as global stability and structural dynamics constraints. To accommodate the requirements of the design process, a new version of the Numerical Manufacturing And Design (NuMAD) code has been developed and released. The code now supports laminate-level blade optimization and an interface to the International Energy Agency Wind Task 37 blade ontology. Transporting long, flexible blades via controlled flapwise bending is found to be a viable approach for blades of up to 100 m. The results confirm that blade mass can be substantially reduced by going either to a downwind design or to a highly coned and tilted upwind design. A discussion of active and inactive constraints consisting of material rupture, fatigue damage, buckling, deflection, and resonant frequencies is presented. An analysis of driving load cases revealed that the downwind designs are dominated by loads from sudden, abrupt events like gusts rather than fatigue. Finally, an analysis of carbon fiber spar caps for downwind machines finds that, compared to typical carbon fibers, the use of a new heavy-tow carbon fiber in the spar caps is found to yield between 9 % and 13 % cost savings.

17 WIND ENERGY↗
NASA NTRS1973

Survey of computational methods for three-dimensional supersonic inviscid flows with shocks

Hyperbolic systems of partial differential equations governing supersonic inviscid flows are discussed and analyzed. Finite-difference analogues for integrating these systems in the interior of fluid domains are described from two points of view: a differential form approach and an integral form approach. The algorithms presented are analyzed for stability and accuracy. The concept of time splitting is discussed and supplied to these methods to achieve increased numerical efficiency. Techniques for treating conditions at the boundaries of the fluid domain and shock-wave discontinuities at surfaces within the domain are described.

Maccormack, R. W.↗
NASA NTRS1983

Viscous, resistive MHD stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one dimensional magnetohydrodynamic (MHD) quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds like numbers. Marginal stability curves, growth rates versus Reynolds like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result which appears general is that instability was found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗
NASA NTRS1983

Heavy-lift airship dynamics

The basic aerodynamic and dynamic properties of an example heavy-lift airship (HLA) configuration are analyzed using a nonlinear, multibody, 6-degrees-of-freedom digital simulation. The slung-payload model is described, and a preliminary analysis of the coupled vehicle-payload dynamics is presented. Trim calculations show the importance of control mixing selection and suggest performance deficiencies in crosswind stationkeeping for the unloaded example HLA. Numerically linearized dynamics of the unloaded vehicle exhibit a divergent yaw mode and an oscillatory pitch mode whose stability characteristic is sensitive to flight speed. An analysis of the vehicle-payload dynamics shows significant coupling of the payload dynamics with those of the basic HLA. It is shown that significant improvement in the vehicle's dynamic behavior can be achieved with the incorporation of a simple flight controller having proportional, rate, and integral-error feedbacks.

Tischler, M. B.↗
NASA NTRS1983

Aeroelastic stability of forward swept composite winged aircraft

This paper reviews the author's past and present aeroelastic stability and performance studies related to forward swept, composite wing aircraft. The influence of laminate elastic bend/twist coupling upon wing divergence, lateral control, and lift effectiveness will be illustrated by means of closed-form solutions, numerical analysis and simple wind-tunnel experiments. In addition, results of analyses of a freely flying flexible FSW aircraft are discussed to indicate the possible effects of the flexible forward swept wing on aircraft dynamic stability. These studies show, both theoretically and experimentally, that, if the aircraft is not carefully designed, a phenomenon referred to as body freedom flutter may appear.

Weisshaar, T. A.↗
NASA NTRS1984

Optimum design of rotor blades for vibration reduction in forward flight

Modern structural optimization techniques are applied to vibration reduction of helicopter rotor blades in forward flight. The objective function minimized consists of the oscillatory vertical hub shears or the hub rolling moments at one particular advance ratio. The behavior constraints are the frequency placements of the blade and the requirement that aeroelastic stability margins, in hover, remain unaffected by the optimization process. The aeroelastic stability and response analysis is based on a fully coupled flap-lag-torsional analysis of the blade. Numerical results are presented for some typical soft-in-plane hingeless rotor configurations indicating a 15-40 percent reduction in vibration levels, as well as a blade which is 20 percent lighter than the initial design. These results imply that structural optimization techniques can yield substantial practical benefits in the design process of rotor systems.

Friedmann, P. P.↗
NASA NTRS1985

Stability of limit cycles in frictionally damped and aerodynamically unstable rotor stages

This paper deals with the stability of limit cycles (Steady-State Oscillations) associated with the multi-degree-of-freedom model of a frictionally damped and aerodynamically unstable rotor stage. By using the first order averaging technique, a generalized criterion has been established to sort out those unstable limit cycles which govern the maximum transient amplitude beyond which the rotor stage becomes unstable.The stability of the remaining steady-state solutions is analyzed by linearizing the averaged system of differential equations. Numerical results are discussed for three-, four- and five-bladed disks.

Sinha, A.↗
NASA NTRS1986

Motion of an electron bunch through a plasma

In space research, pulse-modulated beams, or bunches, of energetic electrons have been used for tracing magnetic field lines and for remote sounding of magnetospheric electric fields. Doubts have been expressed regarding the possibility that an electron bunch would be able to travel through the plasma in space. However, experiments in the ionosphere have shown that, at least on some occasions, bunches of electrons with energies of a few keV have traveled over long distances with no apparent energy loss. The present study is concerned with the theoretical explanation for this phenomenon. The method employed in this study makes use of a direct numerical integration of the Vlasov equation. The results obtained provide the basis for a provisional conclusion concerning the stability of a moving electron bunch in the considered case. It appears that the system is mildly unstable.

Shoucri, M. M.↗
NASA NTRS1989

Numerical Study of Three-dimensional Spatial Instability of a Supersonic Flat Plate Boundary Layer

The behavior of spatially growing three-dimensional waves in a supersonic boundary layer was studied numerically by solving the complete Navier-Stokes equations. Satisfactory comparison with linear parallel and non-parallel stability theories, and experiment are obtained when a small amplitude inflow disturbance is used. The three-dimensional unsteady Navier-Stokes equations are solved by a finite difference method which is fourth-order and second-order accurate in the convection and viscous terms respectively, and second-order accurate in time. Spanwise periodicity is assumed. The inflow disturbance is composed of eigenfunctions from linear stability theory. By increasing the amplitude of the inflow disturbance, nonlinear effects in the form of a relaxation type oscillation of the time signal of rho(u) are observed.

Maestrello, Lucio↗