Approved for public release; distribution is unlimited / A FORTRAN 77 computer code employing an adaptation of the finite differencing algorithm proposed by Brian was developed for the solution of transient heat conduction problems in cylindrical geometries. Validation of code was accomplished by comparison with an ana­lytic solution derived for a model with symmetric, linear boundary conditions. Accuracy of results for asymmetric and non-linear boundary conditions was determined by comparison with a similarly vali­dated code employing the explicit method. Code effectiveness was then demonstrated by conducting a transient temperature analysis for a simulated earth-orbiting satellite. Brian's method demonstrated unconditional stability with associated significant reductions in execu­tion time compared to the explicit method. The effects of discretization error on the accuracy of results require further investigation. / http://archive.org/details/applicationofbri00hein / Lieutenant Commander, United States Navy

The Richtmyer-Meshkov instability occurs when density-stratified interfaces are impulsively accelerated, typically by a shock wave. We present a numerical method to simulate the Richtmyer-Meshkov instability in cylindrical geometry. The ideal MHD equations are linearized about a time-dependent base state to yield linear partial differential equations governing the perturbed quantities. Convergence tests demonstrate that second order accuracy is achieved for smooth flows, and the order of accuracy is between first and second order for flows with discontinuities. Numerical results are presented for cases of interfaces with positive Atwood number and purely azimuthal perturbations. In hydrodynamics, the Richtmyer-Meshkov instability growth of perturbations is followed by a Rayleigh-Taylor growth phase. In MHD, numerical results indicate that the perturbations can be suppressed for sufficiently large perturbation wavenumbers and magnetic fields.

We investigate the Richtmyer-Meshkov instability (RMI) that occurs when an incident shock impulsively accelerates the interface between two different fluids. RMI is important in many technological applications such as Inertial Confinement Fusion (ICF) and astrophysical phenomena such as supernovae. We consider RMI in the presence of the magnetic field in converging geometry through both simulations and analytical means in the framework of ideal magnetohydrodynamics (MHD). In this thesis, we perform linear stability analyses via simulations in the cylindrical geometry, which is of relevance to ICF. In converging geometry, RMI is usually followed by the Rayleigh-Taylor instability (RTI). We show that the presence of a magnetic field suppresses the instabilities. We study the influence of the strength of the magnetic field, perturbation wavenumbers and other relevant parameters on the evolution of the RM and RT instabilities. First, we perform linear stability simulations for a single interface between two different fluids in which the magnetic field is normal to the direction of the average motion of the density interface. The suppression of the instabilities is most evident for large wavenumbers and relatively strong magnetic fields strengths. The mechanism of suppression is the transport of vorticity away from the density interface by two Alfv ́en fronts. Second, we examine the case of an azimuthal magnetic field at the density interface. The most evident suppression of the instability at the interface is for large wavenumbers and relatively strong magnetic fields strengths. After the shock interacts with the interface, the emerging vorticity breaks up into waves traveling parallel and anti-parallel to the magnetic field. The interference as these waves propagate with alternating phase causing the perturbation growth rate of the interface to oscillate in time. Finally, we propose incompressible models for MHD RMI in the presence of normal or azimuthal magnetic field. The linearized equations are solved numerically using inverse Laplace transform. The incompressible models show that the magnetic field suppresses the RMI, and the mechanism of this suppression depends on the orientation of the initially applied magnetic field. The incompressible model agrees reasonably well with compressible linear simulations.

Submitted by Johnny Rodrigues (johnnyrodrigues@ufcg.edu.br) on 2018-04-23T16:23:25Z No. of bitstreams: 1 FERNANDO JOSÉ DE ALMEIDA GAMA - TESE PPGEP 2014..pdf: 3506912 bytes, checksum: c91b94e7df5d978098c9d79aa8eb7487 (MD5) / Made available in DSpace on 2018-04-23T16:23:25Z (GMT). No. of bitstreams: 1 FERNANDO JOSÉ DE ALMEIDA GAMA - TESE PPGEP 2014..pdf: 3506912 bytes, checksum: c91b94e7df5d978098c9d79aa8eb7487 (MD5) Previous issue date: 2014-12-19 / Este trabalho tem como objetivo estudar o fenômeno da difusão transiente de transferência de calor e massa em sólidos com forma geométrica de um cilindro infinito. O estudo apresenta soluções para a equação de difusão com condição de contorno do terceiro tipo. Foram desenvolvidas ferramentas numéricas para a descrição da difusão de calor e massa em produtos com as formas mencionadas. Para as soluções numéricas propostas, a equação de difusão unidimensional foi discretizada usando o método dos volumes finitos, com uma formulação totalmente implícita, usando coordenadas cilíndricas. Para a solução numérica em coordenadas cilíndricas, foram desenvolvidos dois softwares na plataforma Windows, um para a migração de massa e outro para a propagação de calor, utilizando a linguagem de programação Fortran, opção Quick Win Application. O software foi validado usando-se soluções conhecidas para cilindros tanto com parâmetros termofísicos constantes quanto variáveis. Pode-se concluir que as ferramentas desenvolvidas foram adequadas para o estudo de problemas difusivos em geral. As ferramentas desenvolvidas foram usadas para descrever o processo de secagem de bananas inteiras. Na descrição, foram considerados volume e difusividade de calor e massa variáveis. Pode-se concluir que o modelo proposto para descrever o processo apresentou excelentes indicadores estatísticos na descrição da cinética de transferência de calor e massa. Pode-se concluir, também, que a exclusão do aquecimento do vapor nos cálculos efetuados não altera de forma significativa os resultados e que o uso do calor latente da água livre ao invés desta propriedade no produto não produz efeitos significativos. Por outro lado, o desprezo do calor latente de vaporização e a consideração da densidade e do calor específico do produto como propriedades constantes devem ser evitados. / This work aims to study the phenomenon of the transient diffusion of heat and mass in solids with geometric form of an infinite cylinder. The study presents solutions for the diffusion equation with boundary condition of the third kind. Numerical tools for describing the diffusion of heat and mass in the ways mentioned were developed. For the numerical solutions proposed, the one-dimensional diffusion equation was discretized using the finite volume method with a fully implicit formulation, using cylindrical coordinates. For the numerical solution in cylindrical coordinates, two software have been developed on the Windows platform, one for mass migration and one for the heat transfer, using the Fortran programming language, Quick Win Application option. The software was validated using solutions known for cylinders with both constant and variable thermophysical parameters. It can be concluded that the developed tools were appropriate for the study of diffusion problems in general. The above tools were used to describe the process of drying whole banana. In the description, we considered the volume and diffusivities with variables values. It can be concluded that the proposed model to describe the process showed excellent statistical indicators to describe the kinetics of heat and mass transfer. One can also conclude that the exclusion of the vapor heating in the calculations performed does not significantly alter the results. In addition that using the latent heat of free water instead of this property in the product does not produce significant effects. On the other hand, discard the latent heat of vaporization and the consideration of density and specific heat of the product as constant properties should be avoided.