Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Structure Interaction”

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 631 records · Page 35

Efficient simulation of press hardening process through integrated structural and CFD analyses

Press hardened steel parts are being increasingly used in automotive structures for their higher strength to meet safety standards while reducing vehicle weight to improve fuel consumption. However, manufacturing of sheet metal parts by press hardening process to achieve desired properties is extremely challenging as it involves complex interaction of plastic deformation, metallurgical change, thermal distribution, and fluid flow. Numerical simulation is critical for successful design of the process and to understand the interaction among the numerous process parameters to control the press hardening process in order to consistently achieve desired part properties. Until now there has been no integrated commercial software solution that can efficiently model the complete process from forming of the blank, heat transfer between the blank and tool, microstructure evolution in the blank, heat loss from tool to the fluid that flows through water channels in the tools. In this study, a numerical solution based on Altair HyperWorks® product suite involving RADIOSS®, a non-linear finite element based structural analysis solver and AcuSolve®, an incompressible fluid flow solver based on Galerkin Least Square Finite Element Method have been utilized to develop an efficient solution for complete press hardening process design and analysis. RADIOSS is used to handle the plastic deformation, heat transfer between the blank and tool, and microstructure evolution in the blank during cooling. While AcuSolve is used to efficiently model heat loss from tool to the fluid that flows through water channels in the tools. The approach is demonstrated through some case studies.

36 MATERIALS SCIENCE↗

The numerical simulation of subsonic flutter

The present paper describes a numerical simulation of unsteady, subsonic aeroelastic responses. The technique accounts for aerodynamic nonlinearities associated with angles of attack, vortex-dominated flow, static deformations, and unsteady behavior. The fluid and the wing together are treated as a single dynamic system, and the equations of motion for the structure and flowfield are integrated simultaneously and interactively in the time domain. The method employs an iterative scheme based on a predictor-corrector technique. The aerodynamic loads are computed by the general unsteady vortex-lattice method and are determined simultaneously with the motion of the wing. Two models are used to demonstrate the technique: a rigid wing on an elastic support experiencing plunge and pitch about the elastic axis, and a continuous wing rigidly supported at the root chord experiencing spanwise bending and twisting. The time domain solution coupled with the unsteady vortex-lattice method provides the capability of graphically depicting wing and wake motion. Several graphs that illustrate the time domain behavior of the wing and wake are presented.

Strganac, Thomas W.↗

An Assessment of CFD/CSD Prediction State-of-the-Art by Using the HART II International Workshop Data

Over the past decade, there have been significant advancements in the accuracy of rotor aeroelastic simulations with the application of computational fluid dynamics methods coupled with computational structural dynamics codes (CFD/CSD). The HART II International Workshop database, which includes descent operating conditions with strong blade-vortex interactions (BVI), provides a unique opportunity to assess the ability of CFD/CSD to capture these physics. In addition to a baseline case with BVI, two additional cases with 3/rev higher harmonic blade root pitch control (HHC) are available for comparison. The collaboration during the workshop permits assessment of structured, unstructured, and hybrid overset CFD/CSD methods from across the globe on the dynamics, aerodynamics, and wake structure. Evaluation of the plethora of CFD/CSD methods indicate that the most important numerical variables associated with most accurately capturing BVI are a two-equation or detached eddy simulation (DES)-based turbulence model and a sufficiently small time step. An appropriate trade-off between grid fidelity and spatial accuracy schemes also appears to be pertinent for capturing BVI on the advancing rotor disk. Overall, the CFD/CSD methods generally fall within the same accuracy; cost-effective hybrid Navier-Stokes/Lagrangian wake methods provide accuracies within 50% the full CFD/CSD methods for most parameters of interest, except for those highly influenced by torsion. The importance of modeling the fuselage is observed, and other computational requirements are discussed.

Smith, Marilyn J.↗

Direct Analysis of Early Solar System Aqueous Fluids

We lack fundamental information on early solar system aqueous fluids. Fluid inclusions are present in carbonates, sulfides and halides in astro materials, whose characterization is finally becoming feasible. The ordinary chondrite regolith breccias Monahans 1998 (H5) and Zag (H3-6)contain fluid inclusion-bearing halite (NaCl) crystals dated at~4.5 billion years old. Freezing studies to measure the eutectic temperature demonstrated that the brines likely contain divalent cations(Ca, Mg, Fe). Halite is effective at preservation (at least 250 MY) of organic phases and structures. Thus, compositional data on fluid inclusions in halite will reveal unique information on the origin and activity of aqueous fluids in the early solar system, and interactions with organics. To measure fluid compositions from individual fluid inclusions we performed TOF-SIMS depth profiling a tultra-high vacuum at ~110 Kusing a 30 kVBi+analysis ion beam and a1 kVO2+sputtering ion beam. Halite fluid inclusion FI1 contains H2O, K,Na,Cl and NaCl•2H2O.FI2 contains H2O, K, Na, Cl, Al, Fe,Ca,NaCl•2H2O and C-N-O-H species originating from organic molecules. Mass overlap with NaOH species prevented definitive identification of Cain FI1.

M. E Zolensky↗

The delineation and interpretation of the earth's gravity field

In an attempt to understand the mechanical interaction of a growing lithosphere containing fracture zones with small and large scale mantle convection, which gives rise to geoid anomalies in oceanic regions, a series of fluid dynamical experiments is in progress to investigate: (1) the influence of lithosphere structure, fluid depth and viscosity field on the onset, scale, and evolution of sublithospheric convection; (2) the role of this convection in determining the rate of growth of lithosphere, especially in light of the flattening of the lithosphere bathymetry and heat flow at late times; and (3) combining the results of both numerical and laboratory experiments to decide the dominate factors in producing geoid anomalies in oceanic regions through the thermo-mechanical interaction of the lithosphere and subjacent mantle. The clear existence of small scale convection associated with a downward propagating solidification front (i.e., the lithosphere) and a larger scale flow associated with a discontinuous upward heat flux (i.e., a fracture zone) has been shown. The flows exist simultaneously and each may have a significant role in deciding the thermal evolution of the lithosphere and in understanding the relation of shallow mantle convection to deep mantle convection. This overall process is reflected in the geoid, gravity, and topographic anomalies in the north-central Pacific. These highly correlated fields of intermediate wavelength (approx. 200 to 2000 km) show isostatic compensation by a thin lithosphere for shorter (less than or equal to approx. 500 km), but not the longer, wavelengths. The ultimate, dynamic origin of this class of anomalies is being investigated.

Marsh, Bruce D.↗

Molecular dynamics simulations of soliton-like structures in a dusty plasma medium

The existence and evolution of soliton-like structures in a dusty plasma medium are investigated in a first principles approach using molecular dynamic (MD) simulations of particles interacting via a Yukawa potential. These localized structures are found to exist in both weakly and strongly coupled regimes with their structures becoming sharper as the correlation effects between the dust particles get stronger. A surprising result, compared to fluid simulations, is the existence of rarefactive soliton-like structures in our non-dissipative system, a feature that arises from the charge conjugation symmetry property of the Yukawa fluid. Our simulation findings closely resemble many diverse experimental results reported in the past.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Absence of first-order unbinding transitions of fluid and polymerized membranes

Unbinding transitions of fluid and polymerized membranes are studied by renormalization-group (RG) methods. Two different RG schemes are used and found to give rather consistent results. The fixed-point structure of both RG's exhibits a complex behavior as a function of the decay exponent tau for the fluctuation-induced interaction of the membranes. For tau greater than tau(S2) interacting membranes can undergo first-order transitions even in the strong-fluctuation regime. These estimates for tau(S2) imply, however, that both fluid and polymerized membranes unbind in a continuous way in the absence of lateral tension.

Grotehans, Stefan↗

Three-dimensional lattice Boltzmann simulations of microdroplets including contact angle hysteresis on topologically structured surfaces

In this study, the dynamical behavior of a droplet on topologically structured surface is investigated by using a three-dimensional color-gradient lattice Boltzmann model. A wetting boundary condition is proposed to model fluid-surface interactions, which is advantageous to improve the accuracy of the simulation and suppress spurious velocities at the contact line. The model is validated by the droplet partial wetting test and reproduction of the Cassie and Wenzel states. A series of simulations are conducted to investigate the behavior of a droplet when subjected to a shear flow. It is found that in Cassie state, the droplet undergoes a transition from stationary, to slipping and finally to detachment states as the capillary number increases, while in Wenzel state, the last state changes to the breakup state. The critical capillary number, above which the droplet slipping occurs, is small for the Cassie droplet, but is significantly enhanced for the Wenzel droplet due to the increased contact angle hysteresis. In Cassie state, the receding contact angle nearly equals the prediction by the Cassie relation, and the advancing contact angle is close to 180°, leading to a small contact angle hysteresis. In Wenzel state, however, the contact angle hysteresis is extremely large (around 100°). Finally, high droplet mobility can be easily achieved for Cassie droplets, whereas in Wenzel state, extremely low droplet mobility is identified.

97 MATHEMATICS AND COMPUTING↗

Lewis Information Network (LINK): Background and overview

The NASA Lewis Research Center supports many research facilities with many isolated buildings, including wind tunnels, test cells, and research laboratories. These facilities are all located on a 350 acre campus adjacent to the Cleveland Hopkins Airport. The function of NASA-Lewis is to do basic and applied research in all areas of aeronautics, fluid mechanics, materials and structures, space propulsion, and energy systems. These functions require a great variety of remote high speed, high volume data communications for computing and interactive graphic capabilities. In addition, new requirements for local distribution of intercenter video teleconferencing and data communications via satellite have developed. To address these and future communications requirements for the next 15 yrs, a project team was organized to design and implement a new high speed communication system that would handle both data and video information in a common lab-wide Local Area Network. The project team selected cable television broadband coaxial cable technology as the communications medium and first installation of in-ground cable began in the summer of 1980. The Lewis Information Network (LINK) became operational in August 1982 and has become the backbone of all data communications and video.

Schulte, Roger R.↗

A Numerical Model of Unsteady, Subsonic Aeroelastic Behavior

A method for predicting unsteady, subsonic aeroelastic responses was developed. The technique accounts for aerodynamic nonlinearities associated with angles of attack, vortex-dominated flow, static deformations, and unsteady behavior. The fluid and the wing together are treated as a single dynamical system, and the equations of motion for the structure and flow field are integrated simultaneously and interactively in the time domain. The method employs an iterative scheme based on a predictor-corrector technique. The aerodynamic loads are computed by the general unsteady vortex-lattice method and are determined simultaneously with the motion of the wing. Because the unsteady vortex-lattice method predicts the wake as part of the solution, the history of the motion is taken into account; hysteresis is predicted. Two models are used to demonstrate the technique: a rigid wing on an elastic support experiencing plunge and pitch about the elastic axis, and an elastic wing rigidly supported at the root chord experiencing spanwise bending and twisting. The method can be readily extended to account for structural nonlinearities and/or substitute aerodynamic load models. The time domain solution coupled with the unsteady vortex-lattice method provides the capability of graphically depicting wing and wake motion.

Strganac, Thomas W.↗

Simulation-Based Analysis and Prediction of Thrust Vector Servoelastic Coupling

A method of analysis and prediction of servoelastic coupling in launch vehicles is presented, surveying the discovery and subsequent resolution of a predicted servoelastic resonance phenomenon affecting the NASA Space Launch System launch vehicle at specific flight conditions. A physics-based linearized multibody mechanization of the governing equations is combined with first principles analysis to demonstrate that antisymmetric bending of the solid rocket motors leads to a reduction of equivalent viscous modal damping through coupling with the thrust vector control actuators. The sensitivity to parameters and the effects of the resonance phenomenon on flight control performance and stability are confirmed through extensive simulation verification in the time and frequency domain. A novel enhancement in model fidelity that accounts for Coriolis effects of fluid flow on bending within the solid rocket motor case and nozzle is shown to add sufficient damping to reduce the risk of adverse control-structure interaction.

Jeb S Orr↗

Three-body interactions in complex fluids: Virial coefficients from simulation finite-size effects

A simulation technique is described for quantifying the contribution of three-body interactions to the thermodynamical properties of coarse-grained representations of complex fluids. The method is based on a new approach for determining virial coefficients from the measured volume-dependent asymptote of a certain structural function. By comparing the third virial coefficient B{sub 3} for a complex fluid with that of an approximate coarse-grained model described by a pair potential, three body effects can be quantified. The strategy is applicable to both Molecular Dynamics and Monte Carlo simulation. Its utility is illustrated via measurements of three-body effects in models of star polymers and in highly size-asymmetrical colloid-polymer mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Applying Time-Parallelization to Turbulent Flows

Parallelization of the temporal domain is explored for the solution of turbulent flows. Multigrid reduction-in-time (MGRIT) is used to advance the large-scale fluid dynamics in time sequentially on the coarsest space-time grid but propagate the information in time parallel on all other levels. The goal of this process is to accurately and efficiently resolve the coarse-scale turbulence structure and use that to drive the fine-scales of the turbulent flow. The extra forcing from nonlinear multigrid facilitates the coupling and interaction between fine and coarse scales, through which the multiscale nonlinear physics is properly captured. Adaptive mesh refinement is employed to finely resolve only the regions with strong gradients, which provides further computational efficiency. The underlying computational fluid dynamics solver is a fourth-order finite-volume scheme with the standard 4-stage Runge-Kutta method. An advanced approach is devised and implemented to enable MGRIT to solve highly turbulent flows successfully. Furthermore, the method is applied to solve a Taylor-Green vortex problem and a doubleshear-layer turbulent mixing flow. Results are promising, validating that MGRIT with the filtering approach has the potential to efficiently solve general turbulent flows.

Computational Fluid Dynamics↗

A detailed numerical model of a supersonic reacting mixing layer

A current research effort is underway at the NASA Langley Research Center to achieve a detailed understanding of important phenomena present when a supersonic flow undergoes chemical reaction. A computer program has been developed to study the details of such flows. The program has been constructed to consider the multicomponent diffusion and convection of important species, the finite-rate reaction of these species, and the resulting interaction between the fluid mechanics and chemistry. Code results from the analysis of a spatially developing and reacting mixing layer are presented, and conclusions are drawn regarding the structure of the evolving layer and its associated flame.

Drummond, J. P.↗

The coalescence of two merged interaction regions between 6.2 and 9.5 AU - September 1979 event

A simulation study, based on an unsteady, one-dimensional, one-fluid MHD model, and using the plasma and magnetic field data from the Voyager 1 at 6.2 AU, was conducted on the evolution and interaction of solar wind structures to explain the two interaction regions observed by Voyager 1 within a large-scale interplanetary compound stream that was recorded at a heliocentric distance of 6.2 AU. A strong forward shock F(D) with a speed of 960 km/s was present at the front of the second interaction region, and two reverse shocks, R1 and R2, were at the end of the first interaction region. The model shows that the forward shock passed through the two reverse shocks and into the first interaction region, becoming weaker in each of these interactions. The reverse shocks coalesced to form a stronger reverse shock R; thus, the shock signature changed from R1-R2-F(D) to F(D)-R between 6.2 and 9.5 AU. The major stream structures at 9.5 AU predicted by the simulation model agree well with those directly observed from Pioneer 11.

Whang, Y. C.↗

Cosmic string wakes

Accretion of matter onto wakes left behind by horizon-sized pieces of cosmic string is investigated, and the effects of wakes on the large-scale structure of the universe are determined. Accretion of cold matter onto wakes, the effects of a long string on fluids with finite velocity dispersion or sound speeds, the interactions between loops and wakes, and the conditions for wakes to survive disruption by loops are discussed. It is concluded that the most important wakes are those which were formed at the time of equal matter and radiation density. This leads to sheetlike overdense regions of galaxies with a mean separation in agreement with the scale of the bubbles of de Lapparent, Geller, and Huchra (1986). However, for the value of G(mu) favored from galaxy formation considerations in a universe with cold dark matter, a wake accretes matter from a distance of only about 1.5 Mpc, which is much less than the distance between the wakes.

Stebbins, Albert↗

A numerical model for supersonic reacting mixing layers

A current research effort is underway at the NASA Langley Research Center to achieve a detailed understanding of important phenomena present when a supersonic flow undergoes a chemical reaction. A computer program has been developed to study the details of such flows. The program has been constructed to consider the multicomponent diffusion and convection of important species, the finite-rate reaction of these species, and the resulting interaction between the fluid mechanics and chemistry. Code results from the analysis of a spatially developing and reacting mixing layer are presented, and conclusions are drawn regarding the structure of the evolving layer and its associated flame.

Drummond, J. Philip↗

Reflection and transmission of GMIR shock at the heliopause and their relation to the 2- and 3-kHz radio emissions

We use Voyager 2 plasma and magnetic field data together with a one-fluid MHD model to study the interactions of the 1991 Global Merged Interaction Region (GMIR) shock with the heliopause. The 1991 GMIR is an extraordinarily large global solar wind structure in radial, longitudinal and latitudinal extents. It has a strong shock at the leading edge. After its penetration through the termination shock, the GMIR shock first propagates through the subsonic solar wind, then interacts with the heliopause. The interaction produces a transmitted shock propagating outward in the interstellar medium, and a reflected shock propagating backward in the subsonic solar wind. We identify the reflected shock and the transmitted shock as the possible source of the radio noise detected at Voyagers. The plasma frequency behind the reflected and the transmitted shock can be, respectively, responsible for the 2- and 3-kHz radio emissions. The two bands of radio noise are emitted from sources on both sides of the heliopause starting at about the same time. If the emission is generated by f(sub p)-radiation then the heliopause is located at R approximately 130 AU. If the emission is generated by 2f(sub p)-radiation the n R approximately 150 AU. Because the relative speed of the interstellar plasma with respect to the sun appears to be sub-Alfvenic, it is very unlikely there is a fast-bow shock of the heliosphere.

Whang, Y. C.↗