Finite element methods for integrated aerodynamic heating analysis

This report gives a description of the work which has been undertaken during the second year of a three year research program. The objectives of the program are to produce finite element based procedures for the solution of the large scale practical problems which are of interest to the Aerothermal Loads Branch (ALB) at NASA Langley Research Establishment. The problems of interest range from Euler simulations of full three dimensional vehicle configurations to local analyses of three dimensional viscous laminar flow. Adaptive meshes produced for both steady state and transient problems are to be considered. An important feature of the work is the provision of specialized techniques which can be used at ALB for the development of an integrated fluid/thermal/structural modeling capability

Robust topology optimization of three-dimensional photonic-crystal band-gap structures

We perform full 3D topology optimization (in which "every voxel" of the unit cell is a degree of freedom) of photonic-crystal structures in order to find optimal omnidirectional band gaps for various symmetry groups, including fcc (including diamond), bcc, and simple-cubic lattices. Even without imposing the constraints of any fabrication process, the resulting optimal gaps are only slightly larger than previous hand designs, suggesting that current photonic crystals are nearly optimal in this respect. However, optimization can discover new structures, e.g. a new fcc structure with the same symmetry but slightly larger gap than the well known inverse opal, which may offer new degrees of freedom to future fabrication technologies. Furthermore, our band-gap optimization is an illustration of a computational approach to 3D dispersion engineering which is applicable to many other problems in optics, based on a novel semidefinite-program formulation for nonconvex eigenvalue optimization combined with other techniques such as a simple approach to impose symmetry constraints. We also demonstrate a technique for \emph{robust} topology optimization, in which some uncertainty is included in each voxel and we optimize the worst-case gap, and we show that the resulting band gaps have increased robustness to systematic fabrication errors.Comment: 17 pages, 9 figures, submitted to Optics Expres

A discontinuous Galerkin front tracking method for two-phase flows with surface tension

A Discontinuous Galerkin method for solving hyperbolic systems of conservation laws involving interfaces is presented. The interfaces are represented by a collection of element boundaries and their position is updated using an arbitrary Lagrangian-Eulerian method. The motion of the interfaces and the numerical fluxes are obtained by solving a Riemann problem. As the interface is propagated, a simple and effective remeshing technique based on distance functions regenerates the grid to preserve its quality. Compared to other interface capturing techniques, the proposed approach avoids smearing of the jumps across the interface which leads to an improvement in accuracy. Numerical results are presented for several typical two-dimensional interface problems, including flows with surface tension

A new methodology for anisotropic mesh refinement based upon error gradients

We introduce a new strategy for controlling the use of anisotropic mesh refinement based upon the gradients of an a posteriori approximation of the error in a computed finite element solution. The efficiency of this strategy is demonstrated using a simple anisotropic mesh adaption algorithm and the quality of a number of potential a posteriori error estimates is considered

Effect of TNF-α genetic variants and CCR5Δ32 on the vulnerability to HIV-1 infection and disease progression in Caucasian Spaniards

<p>Abstract</p> <p>Background</p> <p>Tumor necrosis factor alpha (TNF-α) is thought to be involved in the various immunogenetic events that influence HIV-1 infection.</p> <p>Methods</p> <p>We aimed to determine whether carriage of the <it>TNF-α-238G>A, -308G>A </it>and <it>-863 C>A </it>gene promoter single nucleotide polymorphisms (SNP) and the <it>CCR5Δ32 </it>variant allele influence the risk of HIV-1 infection and disease progression in Caucasian Spaniards. The study group consisted of 423 individuals. Of these, 239 were uninfected (36 heavily exposed but uninfected [EU] and 203 healthy controls [HC]) and 184 were HIV-1-infected (109 typical progressors [TP] and 75 long-term nonprogressors [LTNP] of over 16 years' duration). <it>TNF-α </it>SNP and the <it>CCR5Δ32 </it>allele were assessed using PCR-RFLP and automatic sequencing analysis methods on white blood cell DNA. Genotype and allele frequencies were compared using the χ 2 test and the Fisher exact test. Haplotypes were compared by logistic regression analysis.</p> <p>Results</p> <p>The distribution of <it>TNF-α-238G>A, -308G>A </it>and <it>-863 C>A </it>genetic variants was non-significantly different in HIV-1-infected patients compared with uninfected individuals: <it>-238G>A</it>, p = 0.7 and p = 0.3; <it>-308G>A</it>, p = 0.05 and p = 0.07; <it>-863 C>A</it>, p = 0.7 and p = 0.4, for genotype and allele comparisons, respectively. Haplotype analyses, however, indicated that carriers of the haplotype H3 were significantly more common among uninfected subjects (p = 0.04). Among the infected patients, the distribution of the three <it>TNF-α </it>genetic variants assessed was non-significantly different between TP and LTNP: <it>-238G>A</it>, p = 0.35 and p = 0.7; <it>-308G>A</it>, p = 0.7 and p = 0.6: <it>-863 C>A</it>, p = 0.2 and p = 0.2, for genotype and allele comparisons, respectively. Haplotype analyses also indicated non-significant associations. Subanalyses in the LTNP subset indicated that the <it>TNF-α-238A </it>variant allele was significantly overrepresented in patients who spontaneously controlled plasma viremia compared with those who had a detectable plasma viral load (genotype comparisons, p = 0.02; allele comparisons, p = 0.03). The <it>CCR5Δ32 </it>distribution was non-significantly different in HIV-1-infected patients with respect to the uninfected population (p = 0.15 and p = 0.2 for genotype and allele comparisons, respectively) and in LTNP vs TP (p = 0.4 and p = 0.5 for genotype and allele comparisons, respectively).</p> <p>Conclusions</p> <p>In our cohort of Caucasian Spaniards, <it>TNF-α </it>genetic variants could be involved in the vulnerability to HIV-1 infection. <it>TNF-α </it>genetic variants were unrelated to disease progression in infected subjects. The <it>-238G>A </it>SNP may modulate the control of viremia in LTNP. Carriage of the <it>CCR5Δ32 </it>variant allele had no effect on the risk of infection and disease progression.</p

EM modelling of arbitrary shaped anisotropic dielectric objects using an efficient 3D leapfrog scheme on unstructured meshes

The standard Yee algorithm is widely used in computational electromagnetics because of its simplicity and divergence free nature. A generalization of the classical Yee scheme to 3D unstructured meshes is adopted, based on the use of a Delaunay primal mesh and its high quality Voronoi dual. This allows the problem of accuracy losses, which are normally associated with the use of the standard Yee scheme and a staircased representation of curved material interfaces, to be circumvented. The 3D dual mesh leapfrog-scheme which is presented has the ability to model both electric and magnetic anisotropic lossy materials. This approach enables the modelling of problems, of current practical interest, involving structured composites and metamaterials

Unstructured Mesh Generation Including Directional Refinement For Aerodynamic Flow Simulation

A method for generating directionally refined unstructured tetrahedral meshes is presented. The motivation is the need to efficiently mesh the complex computational domains which are frequently encountered in aerospace applications. In this context, the method is applicable to the construction of meshes suitable for the simulation of inviscid and viscous aerodynamic flows. Minimum user intervention is required and the user specified stretching distribution is achieved by locally modifying an existing mesh. The interfacing of the mesh generation procedure with available geometry modellers is discussed and certain implementation issues are also addressed. Examples are presented which illustrate the performance of the proposed methodology. Introduction In the area of computational aerodynamics, unstructured mesh technology has developed considerably over the past ten years. In particular, significant progress has been made in the area of automatic tetrahedral mesh generation, so that the..