Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “arbitrary eulerian lagrangian”

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 19 records

Finite Element Simulation of a Space Shuttle Solid Rocket Booster Aft Skirt Splashdown Using an Arbitrary Lagrangian-eulerian Approach

Explicit finite element techniques employing an Arbitrary Lagrangian-Eulerian (ALE) methodology, within the transient dynamic code LS-DYNA, are used to predict splashdown loads on a proposed replacement/upgrade of the hydrazine tanks on the thrust vector control system housed within the aft skirt of a Space Shuttle Solid Rocket Booster. Two preliminary studies are performed prior to the full aft skirt analysis: An analysis of the proposed tank impacting water without supporting aft skirt structure, and an analysis of space capsule water drop tests conducted at NASA's Langley Research Center. Results from the preliminary studies provide confidence that useful predictions can be made by applying the ALE methodology to a detailed analysis of a 26-degree section of the skirt with proposed tank attached. Results for all three studies are presented and compared to limited experimental data. The challenges of using the LS-DYNA ALE capability for this type of analysis are discussed.

Melis, Matthew E.↗

Numerical simulation of the reactive flow in advanced (HSR) combustors using KIVA-2

Recent work has been done with the goal of establishing ultralow emission aircraft gas turbine combustors. A significant portion of the effort is the development of three dimensional computational combustor models. The KIVA-II computer code which is based on the Implicit Continuous Eulerian Difference mesh Arbitrary Lagrangian Eulerian (ICED-ALE) numerical scheme is one of the codes selected by NASA to achieve these goals. This report involves a simulation of jet injection through slanted slots within the Rich burn/Quick quench/Lean burn (RQL) baseline experimental rig. The RQL combustor distinguishes three regions of combustion. This work specifically focuses on modeling the quick quench mixer region in which secondary injection air is introduced radially through 12 equally spaced slots around the mixer circumference. Steady state solutions are achieved with modifications to the KIVA-II program. Work currently underway will evaluate thermal mixing as a function of injection air velocity and angle of inclination of the slots.

Winowich, Nicholas S.↗

Modified ICED-ALE method for astrogeophysical plasma flows

The Implicit-Continuous-Eulerian-Difference Mesh-Arbitrary-Lagrangian-Eulerian (ICED-ALE) algorithm of Brackbill and Pracht (1973) is modified for the study of astrophysical plasma flows in which dynamical effects are important. In the present study the general-energy-conservation law is directly applied to the iteration process, with the total (kinetic-, specific-internal-, and magnetic-) energy density being obtained implicitly at the end of the process. An example is computed in which the convergence speed of the latter method is substantially faster than that of the former. The initializing calculation, or explicit phase, in which the electric current density, magnetic diffusion of field, energy augmentation, and zero-order approximation of the flow velocity are given, is described. Consideration is given to the iteration process, or implicit phase, from which the exact Lagrangian solution for energy density, velocity, and a magnetic field is obtained.

Wu, S. T.↗

The formation of terrains antipodal to major impacts

The SALE, or Simplified Arbitrary Lagrangian Eulerian code is presently used to calculate the antipodal pressures generated by large impacts on icy satellites and planets of various compositions and core sizes. While in a high density core/lower density mantle planet there is a refraction of the incident pressure wave away from the antipode by the core, resulting in a low antipodal pressure, a planet with a low velocity core with a higher-wave speed mantle will focus the waveform to maximize antipodal pressure. Specific impact basins on the moon, Mercury, and several icy satellites are modeled to ascertain whether the impacts were capable of producing the observed antipodal features.

Watts, Alison W.↗

Solution of mixed convection heat transfer from isothermal in-line fins

Transient and steady state combined natural and forced convective flows over two in-line finite thickness fins (louvers) in a vertical channel are numerically solved using two methods. The first method of solution is based on the 'Simple Arbitrary Lagrangian Eulerian' (SALE) technique which incorporates mainly two computational phases: (1) a Lagrangian phase in which the velocity field is updated by the effects of all forces, and (2) an Eulerian phase that executes all advective fluxes of mass, momentum and energy. The second method of solution uses the finite element code entitled FIDAP. In the first part, comparison of the results by FIDAP, SALE, and available experimental work were done and discussed for steady state forced convection over louvered fins. Good agreements were deduced between the three sets of results especially for the flow over a single fin. In the second part and in the absence of experimental literature, the numerical predictions were extended to the transient transports and to the opposing flow where pressure drop is reversed. Results are presented and discussed for heat transfer and pressure drop in assisting and opposing mixed convection flows.

Khalilollahi, Amir↗

Modeling of Non-Spherical Droplet Dynamics

A two-dimensional time-dependent computer code based on the modified Arbitrary Lagrangian Eulerian (ALE) technique, has been developed to simulate non-spherical droplet dynamics and evaporation under convective flows at real rocket combustion chamber conditions. The equations of mass, momentum, energy and species are simultaneously solved for both liquid and gas phases with an accurate dynamic interface tracking. The jump boundary conditions across the deforming droplet surface are obtained by applying the integral forms of conservation of mass, momentum, and energy. At each time step, the interface geometry and flow properties at the droplet surface are implicitly solved by satisfying the interface boundary conditions. A Lagrangian technique was developed to track the arbitrarily moving interface between the liquid droplet and the external gas. An elliptic grid generator is adopted to dynamically reconstruct grids both inside and outside the droplet surface. This code has been used to study droplet oscillation, droplet deformation/breakup, nonspherical droplet evaporation in both low and high pressure convective flows. This presentation briefly describes the numerical algorithm for modeling of the nonspherical droplet dynamics and demonstrates the representative simulation results of nonspherical droplet evaporation at low and high pressure convective flows. Potential applications of this code to rocket combustor design and performance predictions are discussed.

Deng, Zheng-Tao↗

Assessment of antipodal-impact terrains on Mars

The regions anitpodal to Mars' three largest impact basins, Hellas, Isidis, and Argyre, were assessed for evidence of impact-induced disrupted terrains. Photogeology and computer modeling using the Simplified Arbitrary Lagrangian Eulerian (SALE) finite element code suggest that such terrains could have been found by the Hellas impact. Maximum antipodal pressures are 1100 MPa for Hellas, 520 MPa for Isidis, and 150 MPa for Argyre. The results suggest that if antipodal fracturing were associated with later volcanism, then Alba Patera may be related to the Hellas event, as proposed by Peterson (1978). Alba Patera is a unique volcano in the solar system, being a shield volcano which emitted large volume lava flows. This volcanism could be the result of the focusing of seismic energy which created a fractured region that served as a volcanic conduit for the future release of large volumes of magma. No disrupted terrain features are observed antipodal to the Isidis or Argyre basins, although some of the old fractures in Noctis Labyrinthus could have originated in response to the Isidis impact, and later have been reactivated by the Tharsis tectonics assumed to have produced Noctis. If the lower calculated antipodal pressures for Argyre were capable of producing disrupted terrains, then the terrains have been covered subsequently by volcanic or aeolian material, or modified beyond recognition.

Williams, David A.↗

Full-Scale Direct Numerical Simulation of Two- and Three-Dimensional Instabilities and Rivulet Formulation in Heated Falling Films

A thin film draining on an inclined plate has been studied numerically using finite element method. Three-dimensional governing equations of continuity, momentum and energy with a moving boundary are integrated in an arbitrary Lagrangian Eulerian frame of reference. Kinematic equation is solved to precisely update interface location. Rivulet formation based on instability mechanism has been simulated using full-scale computation. Comparisons with long-wave theory are made to validate the numerical scheme. Detailed analysis of two- and three-dimensional nonlinear wave formation and spontaneous rupture forming rivulets under the influence of combined thermocapillary and surface-wave instabilities is performed.

Krishnamoorthy, S.↗

The Role of Multiphysics Simulation in Multidisciplinary Analysis

This article describes the applications of the Spectrum(Tm) Solver in Multidisciplinary Analysis (MDA). Spectrum, a multiphysics simulation software based on the finite element method, addresses compressible and incompressible fluid flow, structural, and thermal modeling as well as the interaction between these disciplines. Multiphysics simulation is based on a single computational framework for the modeling of multiple interacting physical phenomena. Interaction constraints are enforced in a fully-coupled manner using the augmented-Lagrangian method. Within the multiphysics framework, the finite element treatment of fluids is based on Galerkin-Least-Squares (GLS) method with discontinuity capturing operators. The arbitrary-Lagrangian-Eulerian method is utilized to account for deformable fluid domains. The finite element treatment of solids and structures is based on the Hu-Washizu variational principle. The multiphysics architecture lends itself naturally to high-performance parallel computing. Aeroelastic, propulsion, thermal management and manufacturing applications are presented.

Rifai, Steven M.↗

Deployment Simulation Methods for Ultra-Lightweight Inflatable Structures

Two dynamic inflation simulation methods are employed for modeling the deployment of folded thin-membrane tubes. The simulations are necessary because ground tests include gravity effects and may poorly represent deployment in space. The two simulation methods are referred to as the Control Volume (CV) method and the Arbitrary Lagrangian Eulerian (ALE) method. They are available in the LS-DYNA nonlinear dynamic finite element code. Both methods are suitable for modeling the interactions between the inflation gas and the thin-membrane tube structures. The CV method only considers the pressure induced by the inflation gas in the simulation, while the ALE method models the actual flow of the inflation gas. Thus, the transient fluid properties at any location within the tube can be predicted by the ALE method. Deployment simulations of three packaged tube models; namely coiled, Z-folded, and telescopically-folded configurations, are performed. Results predicted by both methods for the telescopically-folded configuration are correlated and computational efficiency issues are discussed.

Wang, John T.↗

Dynamic Deployment Simulations of Inflatable Space Structures

The feasibility of using Control Volume (CV) method and the Arbitrary Lagrangian Eulerian (ALE) method in LSDYNA to simulate the dynamic deployment of inflatable space structures is investigated. The CV and ALE methods were used to predict the inflation deployments of three folded tube configurations. The CV method was found to be a simple and computationally efficient method that may be adequate for modeling slow inflation deployment sine the inertia of the inflation gas can be neglected. The ALE method was found to be very computationally intensive since it involves the solving of three conservative equations of fluid as well as dealing with complex fluid structure interactions.

Wang, John T.↗

Comparison of ALE and SPH Simulations of Vertical Drop Tests of a Composite Fuselage Section into Water

Simulation of multi-terrain impact has been identified as an important research area for improved prediction of rotorcraft crashworthiness within the NASA Subsonic Rotary Wing Aeronautics Program on Rotorcraft Crashworthiness. As part of this effort, two vertical drop tests were conducted of a 5-ft-diameter composite fuselage section into water. For the first test, the fuselage section was impacted in a baseline configuration without energy absorbers. For the second test, the fuselage section was retrofitted with a composite honeycomb energy absorber. Both tests were conducted at a nominal velocity of 25-ft/s. A detailed finite element model was developed to represent each test article and water impact was simulated using both Arbitrary Lagrangian Eulerian (ALE) and Smooth Particle Hydrodynamics (SPH) approaches in LS-DYNA, a nonlinear, explicit transient dynamic finite element code. Analytical predictions were correlated with experimental data for both test configurations. In addition, studies were performed to evaluate the influence of mesh density on test-analysis correlation.

Jackson, Karen E.↗

The formation of antipodal-impact terrains on Mars

The regions antipodal to Mars' three largest impact basins, Hellas, Isidis, and Argyre, were assessed for evidence of impact induced seismically disrupted terrains. Photogeology and preliminary computer modeling using the Simplified Arbitrary Lagrangian-Eulerian computer program suggest such terrain could have been formed by the Hellas and Isidis impacts. Previously, the only assessment of potential surface disruptions on Mars antipodal to large impacts was by Peterson (1978). Herein, the results give credence to Peterson's assertions. For Hellas, the computer model clearly suggests that antipodal pressures were probably strong enough to have fractured the crust and disrupted the surface, and may account for aspects of volcanism at Alba Patera. It is suggested that the Hellas impact produced sufficient focused energy at its antipode to produce deep fracture in the Martian crust centered below the current caldera for Alba Petera. With the generation and evolution of magma in the Tharsis area, the fracture system provided a ready conduit for early stage eruptions of Martians lavas derived from deep in the crust and/or upper mantle, perhaps of komatiitic composition.

Williams, D. A.↗

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods for the analysis of complete jet engines. We have entitled this program by applying the two dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition, and solution capabilities were successfully tested. We then focused attention on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion that results from these structural displacements. This is treated by a new arbitrary Lagrangian-Eulerian (ALE) technique that models the fluid mesh motion as that of a fictitious mass-spring network. New partitioned analysis procedures to treat this coupled three-component problem are developed. These procedures involved delayed corrections and subcycling. Preliminary results on the stability, accuracy, and MPP computational efficiency are reported.

Felippa, C. A.↗

Unsteady Cascade Aerodynamic Response Using a Multiphysics Simulation Code

The multiphysics code Spectrum(TM) is applied to calculate the unsteady aerodynamic pressures of oscillating cascade of airfoils representing a blade row of a turbomachinery component. Multiphysics simulation is based on a single computational framework for the modeling of multiple interacting physical phenomena, in the present case being between fluids and structures. Interaction constraints are enforced in a fully coupled manner using the augmented-Lagrangian method. The arbitrary Lagrangian-Eulerian method is utilized to account for deformable fluid domains resulting from blade motions. Unsteady pressures are calculated for a cascade designated as the tenth standard, and undergoing plunging and pitching oscillations. The predicted unsteady pressures are compared with those obtained from an unsteady Euler co-de refer-red in the literature. The Spectrum(TM) code predictions showed good correlation for the cases considered.

Lawrence, C.↗

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗

Ignition-and-Growth Modeling of NASA Standard Detonator and a Linear Shaped Charge

The main objective of this study is to quantitatively investigate the ignition and shock sensitivity of NASA Standard Detonator (NSD) and the shock wave propagation of a linear shaped charge (LSC) after being shocked by NSD flyer plate. This combined explosive train was modeled as a coupled Arbitrary Lagrangian-Eulerian (ALE) model with LS-DYNA hydro code. An ignition-and-growth (I&G) reactive model based on unreacted and reacted Jones-Wilkins-Lee (JWL) equations of state was used to simulate the shock initiation. Various NSD-to-LSC stand-off distances were analyzed to calculate the shock initiation (or failure to initiate) and detonation wave propagation along the shaped charge. Simulation results were verified by experimental data which included VISAR tests for NSD flyer plate velocity measurement and an aluminum target severance test for LSC performance verification. Parameters used for the analysis were obtained from various published data or by using CHEETAH thermo-chemical code.

Oguz, Sirri↗

TPSAS-NF1676L-18481-DND

Background: - Many of problems of interest involve involve moving or deforming geometries - Governing equations written in Arbitrary Lagrangian-Eulerian (ALE) form to account for grid speed - Nondimensionalization often more involved/confusing/critical Compatibility: - Fully compatible for compressible/incompressible flows; mixed elements; 2D/3D - Not compatible with generic gas model Status: - Compressible path with moving grids is exercised routinely; incompressible path much less so - 6-DOF option has had very limited testing / usage

Bob Biedron↗