Least-Squares Finite Element Formulation for Fluid-Structure Interaction

Fluid-structure interaction problems prove difficult due to the coupling between fluid and solid behavior. Typically, different theoretical formulations and numerical methods are used to solve fluid and structural problems separately. The least-squares finite element method is capable of accurately solving both fluid and structural problems. This capability allows for a simultaneously coupled fluid structure interaction formulation using a single variational approach to solve complex and nonlinear aeroelasticity problems. The least-squares finite element method was compared to commonly used methods for both structures and fluids individually. The fluid analysis was compared to finite differencing methods and the structural analysis type compared to traditional Weak Galerkin finite element methods. The simultaneous solution method was then applied to aeroelasticity problems with a known solution. Achieving these results required unique iterative methods to balance each domain\u27s or differential equation\u27s weighting factor within the simultaneous solution scheme. The scheme required more computational time but it did provide the first hands-off method capable of solving complex fluid-structure interaction problems using a simultaneous least-squares formulation. A sequential scheme was also examined for coupled problems

Proceedings of the ECCOMAS Thematic Conference on Multibody Dynamics 2015

This volume contains the full papers accepted for presentation at the ECCOMAS Thematic Conference on Multibody Dynamics 2015 held in the Barcelona School of Industrial Engineering, Universitat Politècnica de Catalunya, on June 29 - July 2, 2015. The ECCOMAS Thematic Conference on Multibody Dynamics is an international meeting held once every two years in a European country. Continuing the very successful series of past conferences that have been organized in Lisbon (2003), Madrid (2005), Milan (2007), Warsaw (2009), Brussels (2011) and Zagreb (2013); this edition will once again serve as a meeting point for the international researchers, scientists and experts from academia, research laboratories and industry working in the area of multibody dynamics. Applications are related to many fields of contemporary engineering, such as vehicle and railway systems, aeronautical and space vehicles, robotic manipulators, mechatronic and autonomous systems, smart structures, biomechanical systems and nanotechnologies. The topics of the conference include, but are not restricted to: ● Formulations and Numerical Methods ● Efficient Methods and Real-Time Applications ● Flexible Multibody Dynamics ● Contact Dynamics and Constraints ● Multiphysics and Coupled Problems ● Control and Optimization ● Software Development and Computer Technology ● Aerospace and Maritime Applications ● Biomechanics ● Railroad Vehicle Dynamics ● Road Vehicle Dynamics ● Robotics ● Benchmark ProblemsPostprint (published version

Ground state properties of the bond alternating spin-12\frac{1}{2} anisotropic Heisenberg chain

Ground state properties, dispersion relations and scaling behaviour of spin gap of a bond alternating spin-$\frac{1}{2}$ anisotropic Heisenberg chain have been studied where the exchange interactions on alternate bonds are ferromagnetic (FM) and antiferromagnetic (AFM) in two separate cases. The resulting models separately represent nearest neighbour (NN) AFM-AFM and AFM-FM bond alternating chains. Ground state energy has been estimated analytically by using both bond operator and Jordan-Wigner representations and numerically by using exact diagonalization. Dispersion relations, spin gap and several ground state orders have been obtained. Dimer order and string orders are found to coexist in the ground state. Spin gap is found to develop as soon as the non-uniformity in alternating bond strength is introduced in the AFM-AFM chain which further remains non-zero for the AFM-FM chain. This spin gap along with the string orders attribute to the Haldane phase. The Haldane phase is found to exist in most of the anisotropic region similar to the isotropic point.Comment: 16 pages, 6 figures, 1 tabl

Mathematical Framework and Numerical Methods for the Modeling of Mechanochemistry in Multi-Phase Materials

This dissertation presents a suite of mathematical formulations and numerical methods for modeling the interactions between solid mechanics and chemistry in multi-phase materials. In all cases, the treatments rely on the free energy of the system, which potentially includes the strain energy, the chemical free energy, and the interfacial energy. Variational methods are applied to the free energy functionals to derive equilibrium conditions for mechanics and to identify constraints on kinetic laws for chemistry. The applications of this class of variational methods include evolving material configurations associated with phase changes, both diffusive (e.g. oxidation) and non-diffusive (e.g. martensitic transformations). Motivated by the need to represent multi-well, oscillatory, free energy densities, a study is presented comparing spline and polynomial forms for these functions. An alternative approach to phase-field dynamics for finding a minimum energy state is demonstrated, with Mg alloy precipitates as an example. It involves learning the free energy surface as a function of key geometric features with machine learning techniques, which are then used to predict a minimum energy state. This collection of mathematical formulations and numerical methods is aimed at explorations of the physics underlying observed phenomena in multi-phase materials, with potential use in materials' design.PHDMechanical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/138690/1/greght_1.pd

Numerical simulation of the compressible Orszag-Tang vortex 2. Supersonic flow

The numerical investigation of the Orszag-Tang vortex system in compressible magnetofluids will consider initial conditions with embedded supersonic regions. The simulations have initial average Mach numbers 1.0 and 1.5 and beta 10/3 with Lundquist numbers 50, 100, or 200. The behavior of the system differs significantly from that found previously for the incompressible and subsonic analogs. Shocks form at the downstream boundaries of the embedded supersonic regions outside the central magnetic X-point and produce strong local current sheets which dissipate appreciable magnetic energy. Reconnection at the central X-point, which dominates the incompressible and subsonic systems, peaks later and has a smaller impact as M increases from 0.6 to 1.5. Similarly, correlation between the momentum and magnetic field begins significant growth later than in subsonic and incompressible flows. The shocks bound large compression regions, which dominate the wavenumber spectra of autocorrelations in mass density, velocity, and magnetic field

Einstein equations in the null quasi-spherical gauge III: numerical algorithms

We describe numerical techniques used in the construction of our 4th order evolution for the full Einstein equations, and assess the accuracy of representative solutions. The code is based on a null gauge with a quasi-spherical radial coordinate, and simulates the interaction of a single black hole with gravitational radiation. Techniques used include spherical harmonic representations, convolution spline interpolation and filtering, and an RK4 "method of lines" evolution. For sample initial data of "intermediate" size (gravitational field with 19% of the black hole mass), the code is accurate to 1 part in 10^5, until null time z=55 when the coordinate condition breaks down.Comment: Latex, 38 pages, 29 figures (360Kb compressed