Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximated interface conditions”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Ultrasonic waves from radial mode excitation of a piezoelectric disc on the surface of an elastic solid

In this work, we discuss the ultrasonic wave propagation characteristics in an isotropic elastic solid due to the radial mode excitation of a piezoelectric disc actuator bonded to its surface. Finite element simulations using coupled electromechanical modeling is employed to investigate the wave propagation behavior. We find that the radial mode vibrations on the surface of the elastic solid generates all the three types of ultrasonic waves: longitudinal, shear, and surface waves. The waves in the solid are comprised of a central lobe and multiple side lobes based on the frequency of the radial mode excitation. The central lobe is predominantly composed of longitudinal waves and the side lobes are composed of shear waves. While the longitudinal waves have a strong central lobe and weak side lobes that are not fully developed, shear waves consists of fully developed side lobes. In addition, we observe that the longitudinal waves have fewer side lobes inside the solid compared to the shear waves. A semi-analytical approach is presented to explain the above observations. The interface traction boundary condition between the piezoelectric disc and the elastic solid is approximated as a truncated 'Bessel' excitation over the area of the piezoelectric disc in contact with the solid. A Fourier–Bessel series expansion technique is used to investigate the wave propagation behavior from the above traction boundary condition. The results obtained explain the observations from the finite element simulations with regard to the number and distribution of side lobes pertaining to the longitudinal and shear waves generated from the radial modes of the piezoelectric disc. The proposed semi-analytical approach is general and can be applied to any arbitrary axisymmetric excitations on the surface of an elastic solid.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Integration Study of Converter-Interfaced Combined Heat and Power Plants

Combined Heat and Power (CHP) systems, as a proven technology, provide numerous benefits to the plant owner by saving electricity bill, reducing carbon dioxide emissions, and improving electricity service reliability. However, their integration into the distribution grid faces many barriers including the complexity of grid code requirements and the lack of technical expertise in commercial and small industrial plants to interface with the utility and accelerate the interconnection process. To solve these issues, this study proposes a converter-interfaced solution for CHP integration. Simulations including various fault scenarios and load flow under different grid conditions are carried out to compare technical performances of the proposed solution with directly-coupled CHP. Results show that the generator size of converter-interfaced CHP can be reduced by approximately 25% and it also provides advantageous features such as reduced stress level under fault conditions and a better dynamic stability margin. In addition, a power factor controller is designed and reactive power from CHP is dispatched in such a way that the CHP plant is compliant with power factor requirements from utilities under all operating conditions.

Converter, Combined heat and power, Distributed en↗
DOE OSTI2020

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗
DOE OSTI2019

An Agile Design-to-Simulation Workflow Using a New Conforming Moving Least Squares Method

This report summarizes the accomplishments and challenges of a two year LDRD effort focused on improving design-to-simulation agility. The central bottleneck in most solid mechanics simulations is the process of taking CAD geometry and creating a discretization of suitable quality, i.e., the "meshing" effort. This report revisits meshfree methods and documents some key advancements that allow their use on problems with complex geometries, low quality meshes, nearly incompressible materials or that involve fracture. The resulting capability was demonstrated to be an effective part of an agile simulation process by enabling rapid discretization techniques without increasing the time to obtain a solution of a given accuracy. The first enhancement addressed boundary-related challenges associated with meshfree methods. When using point clouds and Euclidean metrics to construct approximation spaces, the boundary information is lost, which results in low accuracy solutions for non-convex geometries and mate rial interfaces. This also complicates the application of essential boundary conditions. The solution involved the development of conforming window functions which use graph and boundary information to directly incorporate boundaries into the approximation space.

42 ENGINEERING↗
DOE OSTI2024

Investigating stainless steel/aluminum bimetallic structures fabricated by cold metal transfer (CMT)-based wire-arc directed energy deposition

Here, this study investigated the process of fabricating a bimetallic structure of 316L stainless steel (SS) and 4043 aluminum (Al) using wire-arc directed energy deposition (DED) based on the cold metal transfer (CMT) process. The impact of heat input on the fabricated structure’s geometry, porosity, and microstructures at the interface, with a specific emphasis on the intermetallic compound (IMC) formation, and the subsequent impact on the joint strength of the structure, were studied. The IMC layer at the interface was predominantly comprised of FeAl 2 Si. For the various heat input conditions studied, the IMC layer's thickness varied from 5 µm to 18 µm. The tensile strength reached up to approximately 130 MPa, which is among the highest reported in the literature for steel/Al bimetallic structures. The specimens fabricated with high heat input conditions had a thicker IMC layer at the steel/Al interface, resulting in a more brittle interface and degradation of the mechanical properties.

36 MATERIALS SCIENCE↗
NASA NTRS2018

Radiative Transfer in a Discrete Random Medium Adjacent to a Half-Space with a Rough Interface

For a macroscopically plane-parallel discrete random medium, the boundary conditions for the specific coherency dyadic at a rough interface are derived. The derivation is based on a modification of the Twersky approximation for a scattering system consisting of a group of particles and the rough surface, and reduces to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave propagating in a discrete random medium with non-scattering boundaries. In a matrix-form setting, the boundary conditions for the specific coherency dyadic imply the boundary conditions for specific intensity column vectors which in turn, yield the expressions for the reflection and transmission matrices. The derived expressions are shown to be identical to those obtained by applying a phenomenological approach based on a facet model to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave.

Adrian Doicu↗
DOE OSTI2020

A passive model for the evolution of subgrid-scale instabilities in turbulent flow regimes

In this work, a zero-dimensional model for the development of turbulent mixing layers driven by Rayleigh–Taylor (RT) or Richtmyer–Meshkov (RM) instabilities is presented, for application in Reynolds-Averaged Navier–Stokes (RANS) simulations. The model approximates the behavior of the turbulent mixing layer at a passively advected tracer particle that sits on the material interface of interest in the simulation, and is used to provide initial conditions for the RANS model once the mixing layer width has developed to a size that may be resolved on the computational mesh. The model relies on approximations to the Besnard–Harlow–Rauenzahn (BHR) RANS closure equations and is thus consistent the turbulent development of the mixing layer, unlike previous models of this type which have focused on the early-time laminar development of RT or RM instabilities. These instabilities will often grow through an initially laminar state before developing into turbulence, and criteria for switching from an existing laminar model to the proposed turbulent model are considered. The proposed model reproduces the mixing layer development of RT and RM instabilities seen in mixed width resolving BHR RANS simulations to a reasonable degree of accuracy, and relaxes the need for the RANS simulation to resolve the mixed width as soon as the mixing layer becomes turbulent, which may be impractical in cases where the mixing layer grows through a wide range of scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Propagation of Spin Waves Through an Interface Between Ferromagnetic and Antiferromagnetic Materials

Boundary conditions for order parameters at an interface between ferromagnetic (FM) and two-sublattice antiferromagnetic (AFM) materials are obtained in the continuous medium approximation similarly to the approach which allows one to take into account the finite thickness of the FM/FM interface, which is much less than spin wavelength. Three order parameters are considered inside an interface of finite thickness with the magnetization M of FM, magnetizations of both sublattices M{sub 1} and M{sub 2} of AFM. The uniform and non-uniform exchange between all order parameters are taken into account to the interface energy. Using these boundary conditions, the excitation of a surface evanescent spin wave is considered in AFM when the spin wave in FM falls onto this interface. The coefficients and the phases of transmission and reflection of spin wave through the FM/AFM interface are derived.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Can We Describe Graphene Confined Water Structures as Overlapping of Approaching Graphene–Water Interfacial Structures?

We investigate the microscopic mechanisms of the overlap of interfacial structures in confined fluids and attempt to answer the question whether the confined structures can be predicted from the original density profiles of individual solid–fluid interfaces. For that purpose, we perform (globally) isobaric–isothermal (locally, grand canonical) molecular dynamics simulations to extract not only the axial distribution functions of the water sites for the uncoupled graphene–water interfaces, but also those corresponding to the confined aqueous environments over the interplate range $8 ≤ h (Å) ≤ 28$ typically at ambient conditions. We have tested two (i.e., an arithmetic and a geometric) superposition approximations for the singlet density of confined water between flat graphene plates, as well as for a combination of flat and corrugated graphene plates. The outcome of this study suggests that the answer to the title’s question is a “yes”, provided that the interplate distance h is large enough to avoid fluid geometric packing frustration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

Application of the MHD energy principle to magnetostatic atmospheres

The MHD energy principle is applied to the stability of a magnetized atmosphere which is bounded below by much denser fluid, as is the solar corona. The two fluids are treated as ideal; the approximation is consistent with the energy principle, and the dynamical conditions that must hold at a fluid-fluid interface are used to show that if vertical displacements of the lower boundary are premitted, then the lower atmosphere must be perturbed as well. However, displacements which do not perturb the coronal boundary can be properly treated as isolated perturbations of the corona alone.

Zweibel, E. G.↗
NASA NTRS1985

Application of the MHD energy principle to magnetostatic atmospheres

The MHD energy principle is applied to the stability of a magnetized atmosphere which is bounded below by much denser fluid, as is the solar corona. The two fluids are treated as ideal; the approximation is consistent with the energy principle, and the dynamical conditions that must hold at a fluid-fluid interface are used to show that if vertical displacements of the lower boundary are permitted, then the lower atmosphere must be perturbed as well. However, displacements which do not perturb the coronal boundary can be properly treated as isolated perturbations of the corona alone.

Zweibel, E. G.↗
DOE OSTI2014

Numerical analysis of flows of rarefied gases in long channels with octagonal cross section shapes

Isothermal, pressure driven rarefied gas flows through long channels with octagonal cross section shapes are analyzed computationally. The capillary is between inlet and outlet reservoirs. The cross section is constant along the axial direction. The boundary condition at the solid-gas interface is assumed to be diffuse reflection. Since the channel is long, the gaseous velocity is small compared to the average molecular speed. Consequently, a linearized description can be used. The flow is described by the linearized Bhatnagar-Gross-Krook kinetic model. The solution of the problem is divided into two stages. First, the local flow field is determined by assuming the local pressure gradient. Secondly, the global flow behavior is deduced by the consideration of the conservation of the mass along the axis of the capillary. The kinetic equation is solved by the discrete velocity method on the cross section. Both spatial and velocity spaces are discretized. A body fitted rectangular grid is used for the spatial space. Near the boundary, first-order, while in the interior part of the flow domain, second-order finite-differences are applied to approximate the spatial derivatives. This combination results into an efficient and straightforward numerical treatment. The velocity space is represented by a Gauss-Legendre quadrature. The kinetic equation is solved in an iterative manner. The local dimensionless flow rate is calculated and tabulated for a wide range of the gaseous rarefaction for octagonal cross sections with various geometrical parameters. It exhibits the Knudsen minimum phenomenon. The flow rates in the octagonal channel are compared to those through capillaries with circular and square cross sections. Typical velocity profiles are also shown. The mass flow rate and the distribution of the pressure are determined and presented for global pressure driven flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Experimental study of vibration isolation in thin-walled structural assemblies with embedded total-internal-reflection metasurfaces

The concept of total-internal-reflection elastic metasurface (TIR-MS) was recently proposed [1] and employed within flexible planar waveguides in order to create highly subwavelength sound-hard barriers impenetrable to low frequency elastic waves. The underlying physical mechanism relies on the design of engineered interfaces exhibiting extreme phase gradients such that any incoming wave at, approximately, any incidence will experience total-internal-reflection conditions. At the design frequency, the metasurface exhibits a large phase gradient such that, in accordance with the generalized Snell's law, the first critical angle is virtually always exceeded. It is worth noting that in practical realizations, the actual total reflection performance might vary depending on the angle of incidence. This dependence is due to the discrete implementation of the metasurface which results in diffraction effects. This paper presents the results of an experimental study that explores the vibration isolation performance of TIR-MS when applied to structures made of complex combinations of different elastic waveguides (e.g. bolted assemblies of beams, plates, and shells). Such system can be seen as a prototypical structure emulating mechanical assemblies of practical interest for many engineering applications. Experimental results confirm that, when the TIR-MS is embedded in the host waveguide, significant vibration isolation capabilities are achieved under quasi-omnidirectional incidence and highly subwavelength excitation conditions (i.e. the ratio of the operating wavelength to the width of the TIR-MS is approximately 5.25). Furthermore, these experimental results suggest new interesting directions to achieve vibration isolation and mechanical energy filtering for practical engineering systems.

42 ENGINEERING↗
NASA NTRS1997

Analysis of Hydrodynamic (Landau) Instability in Liquid-Propellant Combustion at Normal and Reduced Gravity

The burning of liquid propellants is a fundamental combustion problem that is applicable to various types of propulsion and energetic systems. The deflagration process is often rather complex, with vaporization and pyrolysis occurring at the liquid/gas interface and distributed combustion occurring either in the gas phase or in a spray. Nonetheless, there are realistic limiting cases in which combustion may be approximated by an overall reaction at the liquid/gas interface. In one such limit, the gas flame occurs under near-breakaway conditions, exerting little thermal or hydrodynamic influence on the burning propellant. In another such limit, distributed combustion occurs in an intrusive regime, the reaction zone lying closer to the liquid/gas interface than the length scale of any disturbance of interest. Finally, the liquid propellant may simply undergo exothermic decomposition at the surface without any significant distributed combustion, such as appears to occur in some types of HydroxylAmmonium Nitrate (HAN)-based liquid propellants at low pressures. Such limiting models have recently been formulated,thereby significantly generalizing earlier classical models that were originally introduced to study the hydrodynamic stability of a reactive liquid/gas interface. In all of these investigations, gravity appears explicitly and plays a significant role, along with surface tension, viscosity, and, in the more recent models, certain reaction-rate parameters associated with the pressure and temperature sensitivities of the reaction itself. In particular, these parameters determine the stability of the deflagration with respect to not only classical hydrodynamic disturbances, but also with respect to reactive/diffusive influences as well. Indeed, the inverse Froude number, representing the ratio of buoyant to inertial forces, appears explicitly in all of these models, and consequently, in the dispersion relation that determines the neutral stability boundaries beyond which steady, planar burning is unstable to nonsteady, and/or nonplanar (cellular) modes of burning. These instabilities thus lead to a number of interesting phenomena, such as the sloshing type of waves that have been observed in mixtures of HAN and TriEthanolAmmonium Nitrate (TEAN) with water. Although the Froude number was treated as an O(1) quantity in these studies, the limit of small inverse Froude number corresponding to the microgravity regime is increasingly of interest and can be treated explicitly, leading to various limiting forms of the models, the neutral stability boundaries, and, ultimately, the evolution equations that govern the nonlinear dynamics of the propagating reaction front. In the present work, we formally exploit this limiting parameter regime to compare some of the features of hydrodynamic instability of liquid-propellant combustion at reduced gravity with the same phenomenon at normal gravity.

Margolis, Stephen B.↗
DOE OSTI2021

A Poisson–Nernst–Planck Model of Ion Transport and Interface Segregation in Metal–Insulator–Semiconductor Structures and Solar Cells

A numerical model that describes the transport of mobile ionic species in metal–insulator–semiconductor (MIS) and photovoltaic (PV) devices subject to temperature and voltage stress is presented. The finite element method (FEM) is used to solve the Nernst–Planck equation while imposing Poisson's equation self-consistently as a restriction for the electrostatic potential. This allows the contribution of the ionic species to the potential to be taken into account. Using a variational formulation eases the implementation of diverse boundary conditions, including the incorporation of segregation kinetics at the device interfaces. Segregation across the dielectric–semiconductor interface is relevant to modeling the electronic device degradation in systems where contamination reaches the semiconductor. The model in closed systems with no-flux boundary conditions is validated first. In the limiting case of low contamination levels with respect to the gate bias, the FEM solution matches analytically derived approximations. Then, the implementation is broadened to include an open boundary at the dielectric–semiconductor interface to account for leakage of ions. The predicted time dependence of the flatband voltage in Na-contaminated MIS test structures agrees well with measurements. This model successfully captures the role of long-range ion transport at concentrations of relevance to electronic and PV device instability and neuromorphic computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

GUI Control System for the Mu2e Electrostatic Septum High Voltage at Fermilab

The Mu2e Experiment has stringent beam structure requirements; namely, its proton bunches with a time structure of 1.7 $\mu$s in the Fermilab Delivery Ring. This beam structure will be delivered using the Fermilab 8-GeV Booster, the 8-GeV Recycler Ring, and the Delivery Ring. The 1.7-$\mu$s period of the Delivery Ring will generate the required beam structure by means of a third order resonant extraction system operating on a single circulating bunch. The electrostatic septum (ESS) for this system is particularly challenging, requiring mechanical precision in a ultra high vacuum of 1 x 10$^-8$ Torr to generate 100 kV across 15 mm. This paper describes a graphical user interface that has been developed to automate the conditioning and commissioning process for the electrostatic septa. It is based on an interface to the Fermilab ACNET system using the ACSys Python Data Pool Manager (DPM) Client produced and maintained by Fermilab Accelerator Controls. Network interfacing between data pool managers made by the application and ACNET devices introduce an inherent (approximately 1 s) latency in throughput of the readouts. This delay is utilized to process and graph incoming data events of devices crucial to conditioning of a electrostatic septum (ESS). 'Ramping' and 'Monitoring' modes adjust settings of the power supply based on internal logic to efficaciously increase and maintain the high voltage (HV) in the ESS, easing the voltage setting on incidence of sparking or other possibly damaging events. A timestamped log file is produced as the application runs.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

Accurate gradient approximation for complex interface problems in 3D by an improved coupling interface method

Most elliptic interface solvers become complicated for complex interface problems at those “exceptional points” where there are not enough neighboring interior points for high order interpolation. Such complication increases especially in three dimensions. Usually, the solvers are thus reduced to low order accuracy. In this paper, we classify these exceptional points and propose two recipes to maintain order of accuracy there, aiming at improving the previous coupling interface method [26]. Yet the idea is also applicable to other interface solvers. The main idea is to have at least first order approximations for second order derivatives at those exceptional points. Recipe 1 is to use the finite difference approximation for the second order derivatives at a nearby interior grid point, whenever this is possible. Recipe 2 is to flip domain signatures and introduce a ghost state so that a second-order method can be applied. This ghost state is a smooth extension of the solution at the exceptional point from the other side of the interface. The original state is recovered by a post-processing using nearby states and jump conditions. The choice of recipes is determined by a classification scheme of the exceptional points. The method renders the solution and its gradient uniformly second-order accurate in the entire computed domain. Numerical examples are provided to illustrate the second order accuracy of the presently proposed method in approximating the gradients of the original states for some complex interfaces which we had tested previous in two and three dimensions, and a real molecule ( (1D63)) which is double-helix shape and composed of hundreds of atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A mixture-energy-consistent six-equation two-phase numerical model for fluids with interfaces, cavitation and evaporation waves

We model liquid–gas flows with cavitation by a variant of the six-equation single-velocity two-phase model with stiff mechanical relaxation of Saurel–Petitpas–Berry (Saurel et al., 2009) [9]. In our approach we employ phasic total energy equations instead of the phasic internal energy equations of the classical six-equation system. This alternative formulation allows us to easily design a simple numerical method that ensures consistency with mixture total energy conservation at the discrete level and agreement of the relaxed pressure at equilibrium with the correct mixture equation of state. Temperature and Gibbs free energy exchange terms are included in the equations as relaxation terms to model heat and mass transfer and hence liquid–vapor transition. The algorithm uses a high-resolution wave propagation method for the numerical approximation of the homogeneous hyperbolic portion of the model. In two dimensions a fully-discretized scheme based on a hybrid HLLC/Roe Riemann solver is employed. Thermo-chemical terms are handled numerically via a stiff relaxation solver that forces thermodynamic equilibrium at liquid–vapor interfaces under metastable conditions. We present numerical results of sample tests in one and two space dimensions that show the ability of the proposed model to describe cavitation mechanisms and evaporation wave dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗