Interface solutions of partial differential equations with point singularity
Not Available
97 MATHEMATICS AND COMPUTING↗
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.
Pseudospectral collocation methods for fourth order differential equations
Collocation schemes are presented for solving linear fourth order differential equations in one and two dimensions. The variational formulation of the model fourth order problem is discretized by approximating the integrals by a Gaussian quadrature rule generalized to include the values of the derivative of the integrand at the boundary points. Collocation schemes are derived which are equivalent to this discrete variational problem. An efficient preconditioner based on a low-order finite difference approximation to the same differential operator is presented. The corresponding multidomain problem is also considered and interface conditions are derived. Pseudospectral approximations which are C1 continuous at the interfaces are used in each subdomain to approximate the solution. The approximations are also shown to be C3 continuous at the interfaces asymptotically. A complete analysis of the collocation scheme for the multidomain problem is provided. The extension of the method to the biharmonic equation in two dimensions is discussed and results are presented for a problem defined in a nonrectangular domain.
Malek, Alaeddin↗
A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction
This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.
97 MATHEMATICS AND COMPUTING↗
Flight Test Results from the Rake Airflow Gage Experiment on the F-15B
The results are described of the Rake Airflow Gage Experiment (RAGE), which was designed and fabricated to support the flight test of a new supersonic inlet design using Dryden's Propulsion Flight Test Fixture (PFTF) and F-15B testbed airplane (see figure). The PFTF is a unique pylon that was developed for flight-testing propulsion-related experiments such as inlets, nozzles, and combustors over a range of subsonic and supersonic flight conditions. The objective of the RAGE program was to quantify the local flowfield at the aerodynamic interface plane of the Channeled Centerbody Inlet Experiment (CCIE). The CCIE is a fixed representation of a conceptual mixed-compression supersonic inlet with a translating biconic centerbody. The primary goal of RAGE was to identify the relationship between free-stream and local Mach number in the low supersonic regime, with emphasis on the identification of the particular free-stream Mach number that produced a local Mach number of 1.5. Measurements of the local flow angularity, total pressure distortion, and dynamic pressure over the interface plane were also desired. The experimental data for the RAGE program were obtained during two separate research flights. During both flights, local flowfield data were obtained during straight and level acceleration segments out to steady-state test points. The data obtained from the two flights showed small variations in Mach number, flow angularity, and dynamic pressure across the interface plane at all flight conditions. The data show that a free-stream Mach number of 1.65 will produce the desired local Mach number of 1.5 for CCIE. The local total pressure distortion over the interface plane at this condition was approximately 1.5%. At this condition, there was an average of nearly 2 of downwash over the interface plane. This small amount of downwash is not expected to adversely affect the performance of the CCIE inlet.
Frederick, Michael↗
On the coupling of hyperbolic and parabolic systems: Analytical and numerical approach
The coupling of hyperbolic and parabolic systems is discussed in a domain Omega divided into two distinct subdomains omega(+) and omega(-). The main concern is to find the proper interface conditions to be fulfilled at the surface separating the two domains. Next, they are used in the numerical approximation of the problem. The justification of the interface conditions is based on a singular perturbation analysis, i.e., the hyperbolic system is rendered parabolic by adding a small artifical viscosity. As this goes to zero, the coupled parabolic-parabolic problem degenerates into the original one, yielding some conditions at the interface. These are taken as interface conditions for the hyperbolic-parabolic problem. Actually, two alternative sets of interface conditions are discussed according to whether the regularization procedure is variational or nonvariational. It is shown how these conditions can be used in the frame of a numerical approximation to the given problem. Furthermore, a method of resolution is discussed which alternates the resolution of the hyperbolic problem within omega(-) and of the parabolic one within omega(+). The spectral collocation method is proposed, as an example of space discretization (different methods could be used as well); both explicit and implicit time-advancing schemes are considered. The present study is a preliminary step toward the analysis of the coupling between Euler and Navier-Stokes equations for compressible flows.
Gastaldi, Fabio↗
Interface flux recovery framework for constructing partitioned heterogeneous time-integration methods
A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.
97 MATHEMATICS AND COMPUTING↗
Stability of interfaces with mesh refinement
A study is conducted of the stability of mesh refinement in space and time for several different interface equations and finite-difference approximations. First, a root condition which implies stability for the initial-boundary value problem for this type of interface is derived. From the root condition, the stability of several interface equations is proved, using the maximum principle. In some cases, the final verification steps can be done analytically; in other cases, a simple computer program has been written to check the condition for values of a parameter along the boundary of the unit circle. Using this method, stability for Lax-Wendroff with all the interface conditions considered, and for Leapfrog with interpolation interface conditions when the fine and coarse grids overlap is proved.
Berger, M. J.↗
Bubble departure radii at solidification interfaces
A model has been developed for the prediction of bubble departure radii from flat solidification interfaces, including the effect of thermocapillary forces. Under normal gravity conditions, the necessary gas bubble radius required for departure from a CBr4 solidification interface is predicted to be approximately 1/2 mm in agreement with measured values. Under microgravity conditions, however, where surface forces predominate, the model predicts a seemingly prohibitive value of 40 mm. This result is at least in agreement with the microgravity tests conducted on the NASA SPAR I and SPAR III sounding rockets (1978) where the bubbles were not larger than 2 mm in radius and no bubble detachment was observed.
Cole, R.↗
A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction
Herein this paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but computational tests indicate that the methodology introduced here remains stable for selected test cases across a broad range of solid-fluid mass density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.
97 MATHEMATICS AND COMPUTING↗
Stresses and deformations in cross-ply composite tubes subjected to a uniform temperature change: Elasticity and Approximate Solutions
The effects of a uniform temperature change on the stresses and deformations of composite tubes are investigated. The accuracy of an approximate solution based on the principle of complementary virtual work is determined. Interest centers on tube response away from the ends and so a planar elasticity approach is used. For the approximate solution a piecewise linear variation of stresses with the radial coordinate is assumed. The results from the approximate solution are compared with the elasticity solution. The stress predictions agree well, particularly peak interlaminar stresses. Surprisingly, the axial deformations also agree well. This, despite the fact that the deformations predicted by the approximate solution do not satisfy the interface displacement continuity conditions required by the elasticity solution. The study shows that the axial thermal expansion coefficient of tubes with a specific number of axial and circumferential layers depends on the stacking sequence. This is in contrast to classical lamination theory which predicts the expansion to be independent of the stacking arrangement. As expected, the sign and magnitude of the peak interlaminar stresses depends on stacking sequence.
Hyer, M. W.↗
Stresses and deformations in cross-ply composite tubes subjected to a uniform temperature change
This study investigates the effects of a uniform temperature change on the stresses and deformations of composite tubes and determines the accuracy of an approximate solution based on the principle of complementary virtual work. Interest centers on tube response away from the ends and so a planar elasticity approach is used. For the approximate solution a piecewise linear variation of stresses with the radial coordinate is assumed. The results from the approximate solution are compared with the elasticity solution. The stress predictions agree well, particularly peak interlaminar stresses. Surprisingly, the axial deformations also agree well, despite the fact that the deformations predicted by the approximate solution do not satisfy the interface displacement continuity conditions required by the elasticity solution. The study shows that the axial thermal expansion coefficient of tubes with a specific number of axial and circumferential layers depends on the stacking sequence. This is in contrast to classical lamination theory, which predicts that the expansion will be independent of the stacking arrangement. As expected, the sign and magnitude of the peak interlaminar stresses depend on stacking sequence. For tubes with a specific number of axial and circumferential layers, thermally induced interlaminar stresses can be controlled by altering stacking arrangement.
Hyer, M. W.↗
An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing
The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Development and Validation of a Slosh-Based Ullage Collapse Model
Ullage collapse modeling is an essential tool for the success of spacecraft powered by cryogenic propellants. The degree of ullage collapse can have a major impact on the design of cryogenic propellant tanks as well as the mission’s concept of operations. Additional commodities (e.g., helium gas for pressure control) and their quantities are dependent on the amount of heat transfer and phase change of propellants expected during flight. Currently, high-fidelity computational fluid dynamics (CFD) with heat transfer and phase change is the most accurate method for predicting ullage collapse. CFD simulations of this kind, however, can become very costly for problems requiring large domains, long simulation times, or significant liquid motion. Consequently, there is a need for a reduced order model to provide fast responses which can be used in engineering design decisions. The proposed numerical model is a nodal approach, augmented by outputs from CFD analyses that are generated rapidly and reliably. These outputs include transient values for liquid-vapor interface area, interface velocity, and liquid position using a volume of fluid (VoF) method. There are several benefits of combining a nodal approach with slosh predictions from CFD simulations: leveraging of the Loci/STREAM-VoF solver that has been validated in high-g and low-g environments, rapid modeling of VoF with isothermal propellant simplifications and relatively coarse meshes, and improved model predictions from the nodal tool. The nodal approach, informed by the slosh predictions from CFD, assumes uniform conditions in the liquid and gas domains. Modifications to account for temperature stratification in the liquid and gas may be implemented as a future improvement. Thermodynamic changes within the tank are solved at user-defined timesteps. First, impacts due to mass entering or leaving the control volume are captured. Second, heat transfer impacts are accounted for by tracking the volume of the gas which has been contacted by liquid in that timestep and rebalancing the internal energy. Finally, phase change is approximated using the area and velocity of the interface. The tank conditions are re-computed by conserving mass and energy. Calibration of the model to several test flights and experiments are required to baseline the predictive power and uncertainty. The final product will enable designers and engineers to create quick-turnaround predictions for ullage collapse and inform design feasibility.
Computational Fluid Dynamics↗
Theoretical and Experimental Investigations on Droplet Evaporation and Droplet Ignition at High Pressures
The evaporation of fuel droplets under high ambient pressure and temperature in normal gravity and microgravity has been investigated experimentally. For subcritical ambient conditions, droplet evaporation after a heat-up period follows the d(exp 2)-law. For all data the evaporation constant increases as the ambient temperature increases. At identical ambient conditions the evaporation constant under microgravity is smaller compared to normal gravity. This effect can first be observed at 1 bar and increases with ambient pressure. Preliminary experiments on ignition delay for self-igniting fuel droplets have been performed. Above a 1 s delay time, at identical ambient conditions, significant differences in the results of the normal and microgravity data are observed. Self-ignition occurs within different temperature ranges due to the influence of gravity. The time dependent behavior of the droplet is examined theoretically. In the calculations two different approaches for the gas phase are applied. In the first approach the conditions at the interface are given using a quasi steady theory approximation. The second approach uses a set of time dependent governing equations for the gas phase which are then evaluated. In comparison, the second model shows a better agreement with the drop tower experiments. In both cases a time dependent gasification rate is observed.
Ristau, R.↗
The optical depth of the 158 micrometer (C-12 II) line: Detection of the F=1 yields 0 (C-13 III) hyperfine-structure component
The first detection of the F = 1 yields 0 hyperfine component of the 158 micrometer (C-13 II) fine structure line in the interstellar medium is reported. A twelve point intensity map was obtained of the (C-13 II) distribution over the inner 190 inch (right ascension) by 190 inch (declination) regions of the Orion nebula using an imaging Fabry-Perot interferometer. The (C-12 II)/(C-13 II) line intensity ratio varied significantly over the region mapped. It is highest (86 plus or minus 9) in the core of the Orion H II region and significantly lower (62 plus or minus 7) in the outer regions of the map, reflecting higher optical depth in the (C-12 II) line here. It is suggested that this enhanced optical depth is the result of limb brightening of the optically thin (C-13 II) line at the edges of the bowl-shaped H II region blister. If the C-12/C-13 abundance ratio is 43, the (C-12 II) line in the inner regions of the Orion nebula, has a low optical depth: tau sub 12 approximately = 0.75 plus or minus 0.25. The optical depth together with the large brightness temperature of the (C-12 II) line (approximately 160 K) requires that the excitation temperature of the P-2 sub 3/2 level be approximately 310 K, in very good agreement with the previous analysis of the physical conditions of the Orion interface region based on fine structure line intensity ratios and photodissociation region models. If the C-12/C-13 abundance ratio is 67, the line optical depth is somewhat larger (tau sub 12 approximately = 1.85), and the transition excitation temperature is somewhat smaller (approximately 190 K) than that predicted by these models. The present results therefore support values approximately = 43 for the C-12/C-13 abundance ratio in the Orion nebula.
Stacey, G. J.↗
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↗
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↗
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↗