Simulations of Normal Spiral Galaxies

Results are presented of numerical simulations of normal isolated late type spiral galaxies. Specifically the galaxy NGC 628 is used as a template. The method employs a TREESPH code including stellar particles, gas particles, cooling and heating of the gas, star formation according to a Jeans criterion, and Supernova feedback. A regular spiral disc can be generated as an equilibrium situation of two opposing actions. On the one hand cooling and dissipation of the gas, on the other hand gas heating by the FUV field of young stars and SN mechanical forcing. The disc exhibits small and medium scale spiral structure of which the multiplicity increases as a function of radius. The theory of swing amplification can explain, both qualitatively and quantitatively, the emerging spiral structure. In addition, swing amplification predicts that the existence of a grand design m=2 spiral is only possible if the disc is massive. The simulations show that the galaxy is then unstable to bar formation. A general criterion is derived for the transition between bar stable and unstable, depending on disc mass contribution and on disc thickness. It seems that bar stability hardly depends on the presence of gas. A detailed quantitative analysis is made of the emerging spiral structure and a comparison is made with observations. That demonstrates that the structure of the numerical isolated galaxies is not as strong and has a larger multiplicity compared to the structure of some exemplary real galaxies. It is argued that a grand design can only be generated by a central bar or by tidal forces resulting from an encounter with another galaxy.Comment: Accepted for publication in: Monthly Notices of the Royal Astron. Soc. The astro-ph copy has a few figures with degraded resolution. A copy with high quality graphics (4281 kb) can be downloaded from the kapteyn institute weg page at http://www.astro.rug.nl (Goto preprints and preprints 2003

The opaque nascent starburst in NGC 1377: Spitzer SINGS observations

We analyze extensive data on NGC 1377 from the Spitzer Infrared Nearby Galaxies Survey (SINGS). Within the category of nascent starbursts that we previously selected as having infrared-to-radio continuum ratios in large excess of the average and containing hot dust, NGC 1377 has the largest infrared excess yet measured. Optical imaging reveals a morphological distortion suggestive of a recent accretion event. Infrared spectroscopy reveals a compact and opaque source dominated by a hot, self-absorbed continuum (τ~ 20 in the 10μm silicate band). We provide physical evidence against nonstellar activity being the heating source. H II regions are detected through the single [Ne II] line, probing 85% of ionizing photons are suppressed by dust. The only other detected emission features are molecular hydrogen lines, arguably excited mainly by shocks, besides photodissociation regions, and weak aromatic bands. The new observations support our interpretation in terms of an extremely young starburst (<1 Myr). More generally, galaxies deficient in radio synchrotron emission are likely observed within a few Myr of the onset of a starburst and after a long quiescence, prior to the replenishment of the interstellar medium with cosmic rays. The similar infrared-radio properties of NGC 1377 and some infrared-luminous galaxies suggest that NGC 1377 constitutes an archetype that will be useful to better understand starburst evolution. Although rare locally because observed in a brief evolutionary stage, nascent starbursts may represent a nonnegligible fraction of merger-induced starbursts that dominate deep infrared counts. Since they differ dramatically from usual starburst templates, they have important consequences for the interpretation of deep surveys

Modeling, optimization, and sensitivity analysis of a continuous multi-segment crystallizer for production of active pharmaceutical ingredients

We have investigated the simulation-based, steady-state optimization of a new type of crystallizer for the production of pharmaceuticals. The multi-segment, multi-addition plug-flow crystallizer (MSMA-PFC) offers better control over supersaturation in one dimension compared to a batch or stirred-tank crystallizer. Through use of a population balance framework, we have written the governing model equations of population balance and mass balance on the crystallizer segments. The solution of these equations was accomplished through either the method of moments or the finite volume method. The goal was to optimize the performance of the crystallizer with respect to certain quantities, such as maximizing the mean crystal size, minimizing the coefficient of variation, or minimizing the sum of the squared errors when attempting to hit a target distribution. Such optimizations are all highly nonconvex, necessitating the use of the genetic algorithm. Our results for the optimization of a process for crystallizing flufenamic acid showed improvement in crystal size over prior literature results. Through the use of a novel simultaneous design and control (SDC) methodology, we have further optimized the flowrates and crystallizer geometry in tandem.^ We have further investigated the robustness of this process and observe significant sensitivity to error in antisolvent flowrate, as well as the kinetic parameters of crystallization. We have lastly performed a parametric study on the use of the MSMA-PFC for in-situ dissolution of fine crystals back into solution. Fine crystals are a known processing difficulty in drug manufacture, thus motivating the development of a process that can eliminate them efficiently. Prior results for cooling crystallization indicated this to be possible. However, our results show little to no dissolution is used after optimizing the crystallizer, indicating the negative impact of adding pure solvent to the process (reduced concentration via dilution, and decreased residence time) outweighs the positive benefits of dissolving fines. The prior results for cooling crystallization did not possess this coupling between flowrate, residence time, and concentration, thus making fines dissolution significantly more beneficial for that process. We conclude that the success observed in hitting the target distribution has more to do with using multiple segments and having finer control over supersaturation than with the ability to go below solubility. Our results showed that excessive nucleation still overwhelms the MSMA-PFC for in-situ fines dissolution when nucleation is too high