An Unusual Application of NASTRAN Contour Plotting Capability

A procedure is presented for obtaining contour plots of any physical quantity defined on a number of points of the surface of a structure. Rigid Format 1 of HEAT approach in Cosmic NASTRAN is ALTERED to enable use of contour plotting capability for scalar quantities. The ALTERED DMAP sequence is given. Examples include temperature distribution on the face of a cooled laser mirror and the angle of incidence or a radome surface

Constraining star formation rates in cool-core brightest cluster galaxies

We used broad-band imaging data for 10 cool-core brightest cluster galaxies (BCGs) and conducted a Bayesian analysis using stellar population synthesis to determine the likely properties of the constituent stellar populations. Determination of ongoing star formation rates (SFRs), in particular, has a direct impact on our understanding of the cooling of the intracluster medium (ICM), star formation and AGN-regulated feedback. Our model consists of an old stellar population and a series of young stellar components. We calculated marginalized posterior probability distributions for various model parameters and obtained 68% plausible intervals from them. The 68% plausible interval on the SFRs is broad, owing to a wide range of models that are capable of fitting the data, which also explains the wide dispersion in the star formation rates available in the literature. The ranges of possible SFRs are robust and highlight the strength in such a Bayesian analysis. The SFRs are correlated with the X-ray mass deposition rates (the former are factors of 4 to 50 lower than the latter), implying a picture where the cooling of the ICM is a contributing factor to star formation in cool-core BCGs. We find that 9 out of 10 BCGs have been experiencing starbursts since 6 Gyr ago. While four out of 9 BCGs seem to require continuous SFRs, 5 out of 9 seem to require periodic star formation on intervals ranging from 20 Myr to 200 Myr. This time scale is similar to the cooling-time of the ICM in the central (< 5 kpc) regions.Comment: 33 pages, 14 Figures, 14 Tables. Accepted for publication in MNRA

The challenging task of determining star formation rates: the case of a massive stellar burst in the brightest cluster galaxy of Phoenix galaxy cluster

Star formation in galaxies at the center of cooling-flow galaxy clusters is an important phenomenon in the context of formation and evolution of massive galaxies in the Universe. Yet, star formation rates (SFRs) in such systems continue to be elusive. We use our Bayesian-motivated spectral energy distribution (SED)-fitting code, BAYESCOOL, to estimate the plausible SFR values in the brightest cluster galaxy of a massive, X-ray luminous galaxy cluster, Phoenix. Previous studies of Phoenix have resulted in the highest measurement of SFR for any galaxy, with the estimates reaching up to 1000 solar masses/yr. However, a very small number of models have been considered in those studies. BAYESCOOL allows us to probe a large parameter space. We consider two models for star formation history, instantaneous bursts and continuous star formation, a wide range of ages for the old and the young stellar population, along with other discrete parameters, such as the initial mass function, metallicities, internal extinction and extinction law. We find that in the absence of any prior except that the maximum cooling rate < 3000 solar masses/yr, the SFR lies in the range (2230-2890) solar masses/yr. If we impose an observational prior on the internal extinction, E(B-V) < 0.6, the best-fit SFR lies in (454-494) solar masses/yr, and we consider this as the most probable range of SFR values for Phoenix. The SFR dependence on the extinction is a reflection of the standard age-extinction degeneracy, which can be overcome by using a prior on one of the two quantities in question.Comment: 12 pages, 4 figures, 1 Table, accepted for publication in MNRA

Parallel finite element simulation of 3d incompressible flows: fluid-structure interactions

Massively parallel finite element computations of 3D, unsteady incompressible flows, including those involving fluid-structure interactions, are presented. The computation with time-varying spatial domains are based on the deforming spatial domain/stabilized space-time (DSD/SST) finite element formulation. The capability to solve 3D problems involving fluid-structure interactions is demonstrated by investigating the dynamics of a flexible cantilevered pipe conveying fluid. Computations of flow past a stationary rectangular wing at Reynolds number 1000, 2500 and 107 reveal interesting flow patterns. In these computations, at each time step approximately 3 &#215; 106 non-linear equations are solved to update the flow field. Also, preliminary results are presented for flow past a wing in flapping motion. In this case a specially designed mesh moving scheme is employed to eliminate the need for remeshing. All these computations are carried out on the Army High Performance Computing Research Center supercomputers CM-200 and CM-5, with major speed-ups compared with traditional supercomputers. The coupled equation systems arising from the finite element discretizations of these large-scale problems are solved iteratively with diagonal preconditioners. In some cases, to reduce the memory requirements even further, these iterations are carried out with a matrix-free strategy. The finite element formulations and their parallel implementations assume unstructured meshes

A finite element study of incompressible flows past oscillating cylinders and aerofoils

We present our numerical results for certain unsteady flows past oscillating cylinders and aerofoils. The computations are based on the stabilized space-time finite element formulation. The implicit equation systems resulting from the space-time finite element discretizations are solved using iterative solution techniques. One of the problems studied is flow past a cylinder which is forced to oscillate in the horizontal direction. In this case we observe a change from an unsymmetric mode of vortex shedding to a symmetric one. An extensive study was carried out for the case in which a cylinder is mounted on lightly damped springs and allowed to oscillate in the vertical direction. In this case the motion of the cylinder needs to be determined as part of the solution, and under certain conditions this motion changes the vortex-shedding pattern of the flow field significantly. This non-linear fluid-structure interaction exhibits certain interesting behaviour such as 'lock-in' and 'hysteresis', which are in good agreement with the laboratory experiments carried out by other researchers in the past. Preliminary results for flow past a pitching aerofoil are also presented

Magnetic order in CaFe1-xCoxAsF (x = 0, 0.06, 0.12) superconductor compounds

A Neutron Powder Diffraction (NPD) experiment has been performed to investigate the structural phase transition and magnetic order in CaFe1-xCoxAsF superconductor compounds (x = 0, 0.06, 0.12). The parent compound CaFeAsF undergoes a tetragonal to orthorhombic phase transition at 134(3) K, while the magnetic order in form of a spin-density wave (SDW) sets in at 114(3) K. The antiferromagnetic structure of the parent compound has been determined with a unique propagation vector k = (1,0,1) and the Fe saturation moment of 0.49(5)uB aligned along the long a-axis. With increasing Co doping, the long range antiferromagnetic order has been observed to coexist with superconductivity in the orthorhombic phase of the underdoped CaFe0.94Co0.06AsF with a reduced Fe moment (0.15(5)uB). Magnetic order is completely suppressed in optimally doped CaFe0.88Co0.12AsF. We argue that the coexistence of SDW and superconductivity might be related to mesoscopic phase separation.Comment: 4pages, 4figure

Thermodynamics predicts how confinement modifies hard-sphere dynamics

We study how confining the equilibrium hard-sphere fluid to restrictive one- and two-dimensional channels with smooth interacting walls modifies its structure, dynamics, and entropy using molecular dynamics and transition-matrix Monte Carlo simulations. Although confinement strongly affects local structuring, the relationships between self-diffusivity, excess entropy, and average fluid density are, to an excellent approximation, independent of channel width or particle-wall interactions. Thus, thermodynamics can be used to predict how confinement impacts dynamics.Comment: 4 pages, 4 figure

Parallel computation of unsteady compressible flows with the EDICT

Recently, the Enhanced-Discretization Interface-Capturing Technique (EDICT) was introduced for simulation of unsteady flow problems with interfaces such as two-fluid and free-surface flows. The EDICT yields increased accuracy in representing the interface. Here we extend the EDICT to simulation of unsteady viscous compressible flows with boundary/shear layers and shock/expansion waves. The purpose is to increase the accuracy in selected regions of the computational domain. An error indicator is used to identify these regions that need enhanced discretization. Stabilized finite-element formulations are employed to solve the Navier-Stokes equations in their conservation law form. The finite element functions corresponding to enhanced discretization are designed to have two components, with each component coming from a different level of mesh refinement over the same computational domain. The primary component comes from a base mesh. A subset of the elements in this base mesh are identified for enhanced discretization by utilizing the error indicator. A secondary, more refined, mesh is constructed by patching together the second-level meshes generated over this subset of elements, and the second component of the functions comes from this mesh. The subset of elements in the base mesh that form the secondary mesh may change from one time level to other depending on the distribution of the error in the computations. Using a parallel implementation of this EDICT-based method, we apply it to test problems with shocks and boundary layers, and demonstrate that this method can be used very effectively to increase the accuracy of the base finite element formulation

Enhancement of bulk second-harmonic generation from silicon nitride films by material composition

We present a comprehensive tensorial characterization of second-harmonic generation from silicon nitride films with varying composition. The samples were fabricated using plasma-enhanced chemical vapor deposition, and the material composition was varied by the reactive gas mixture in the process. We found a six-fold enhancement between the lowest and highest second-order susceptibility, with the highest value of approximately 5 pm/V from the most silicon-rich sample. Moreover, the optical losses were found to be sufficiently small (below 6 dB/cm) for applications. The tensorial results show that all samples retain in-plane isotropy independent of silicon content, highlighting the controllability of the fabrication process.Comment: 4 pages, 3 figures, 2 tables; Re-submitted to Optics Letter