Numerical analysis of nonlinear pneumatic structures

Numerical analysis of nonlinear behavior of inflatable structure

A hybrid-stress finite element for linear anisotropic elasticity

Standard assumed displacement finite elements with anisotropic material properties perform poorly in complex stress fields such as combined bending and shear and combined bending and torsion. A set of three dimensional hybrid-stress brick elements were developed with fully anisotropic material properties. Both eight-node and twenty-node bricks were developed based on the symmetry group theory of Punch and Atluri. An eight-node brick was also developed using complete polynomials and stress basis functions and reducing the order of the resulting stress parameter matrix by applying equilibrium constraints and stress compatibility constraints. Here the stress compatibility constraints must be formulated assuming anisotropic material properties. The performance of these elements was examined in numerical examples covering a broad range of stress distributions. The stress predictions show significant improvement over the assumed displacement elements but the calculation time is increased

Thermo-viscoplastic analysis of hypersonic structures subjected to severe aerodynamic heating

A thermoviscoplastic computational method for hypersonic structures is presented. The method employs unified viscoplastic constitutive model implemented in a finite element approach for quasi-static thermal-structural analysis. Applications of the approach to convectively cooled hypersonic structures illustrate the effectiveness of the approach and provide insight into the transient inelastic structural behavior at elevated temperatures

Solid rocket booster internal flow analysis by highly accurate adaptive computational methods

The primary objective of this project was to develop an adaptive finite element flow solver for simulating internal flows in the solid rocket booster. Described here is a unique flow simulator code for analyzing highly complex flow phenomena in the solid rocket booster. New methodologies and features incorporated into this analysis tool are described

H-P adaptive methods for finite element analysis of aerothermal loads in high-speed flows

The commitment to develop the National Aerospace Plane and Maneuvering Reentry Vehicles has generated resurgent interest in the technology required to design structures for hypersonic flight. The principal objective of this research and development effort has been to formulate and implement a new class of computational methodologies for accurately predicting fine scale phenomena associated with this class of problems. The initial focus of this effort was to develop optimal h-refinement and p-enrichment adaptive finite element methods which utilize a-posteriori estimates of the local errors to drive the adaptive methodology. Over the past year this work has specifically focused on two issues which are related to overall performance of a flow solver. These issues include the formulation and implementation (in two dimensions) of an implicit/explicit flow solver compatible with the hp-adaptive methodology, and the design and implementation of computational algorithm for automatically selecting optimal directions in which to enrich the mesh. These concepts and algorithms have been implemented in a two-dimensional finite element code and used to solve three hypersonic flow benchmark problems (Holden Mach 14.1, Edney shock on shock interaction Mach 8.03, and the viscous backstep Mach 4.08)

A dual weighted residual method applied to complex periodic gratings

An extension of the dual weighted residual (DWR) method to the analysis of electromagnetic waves in a periodic diffraction grating is presented. Using the α,0-quasi-periodic transformation, an upper bound for the a posteriori error estimate is derived. This is then used to solve adaptively the associated Helmholtz problem. The goal is to achieve an acceptable accuracy in the computed diffraction efficiency while keeping the computational mesh relatively coarse. Numerical results are presented to illustrate the advantage of using DWR over the global a posteriori error estimate approach. The application of the method in biomimetic, to address the complex diffraction geometry of the Morpho butterfly wing is also discussed

A Cahn-Hilliard-Darcy model for tumour growth with chemotaxis and active transport

Using basic thermodynamical principles we derive a Cahn--Hilliard--Darcy model for tumour growth including nutrient diffusion, chemotaxis, active transport, adhesion, apoptosis and proliferation. The model generalise earlier models and in particular include active transport mechanisms which ensures thermodynamical consistency. We perform a formally matched asymptotic expansion and develop several sharp interface models. Some of them are classical and some new ones which for example include a jump in the nutrient density at the interface. A linear stability analysis for a growing nucleus is performed and in particular the role of the new active transport term is analysed. Numerical computations are performed to study the influence of the active transport term for specific growth scenarios

Simulation of actively controlled spacecraft with flexible appendages

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76688/1/AIAA-25388-716.pd

Gyroscopic motion of superfluid trapped atomic condensates

The gyroscopic motion of a trapped Bose gas containing a vortex is studied. We model the system as a classical top, as a superposition of coherent hydrodynamic states, by solution of the Bogoliubov equations, and by integration of the time-dependent Gross-Pitaevskii equation. The frequency spectrum of Bogoliubov excitations, including quantum frequency shifts, is calculated and the quantal precession frequency is found to be consistent with experimental results, though a small discrepancy exists. The superfluid precession is found to be well described by the classical and hydrodynamic models. However the frequency shifts and helical oscillations associated with vortex bending and twisting require a quantal treatment. In gyroscopic precession, the vortex excitation modes $m=\pm 1$ are the dominant features giving a vortex kink or bend, while the $m=+2$ is found to be the dominant Kelvin wave associated with vortex twisting.Comment: 18 pages, 7 figures, 1 tabl

Reapportionment of Communicative Burden in Aphasia: A Study of Narrative Interactions

TB