Weakly Enforced Boundary Conditions for the NURBS-Based Finite Cell Method

In this paper, we present a variationally consistent formulation for the weak enforcement of essential boundary conditions as an extension to the finite cell method, a fictitious domain method of higher order. The absence of boundary fitted elements in fictitious domain or immersed boundary methods significantly restricts a strong enforcement of essential boundary conditions to models where the boundary of the solution domain coincides with the embedding analysis domain. Penalty methods and Lagrange multiplier methods are adequate means to overcome this limitation but often suffer from various drawbacks with severe consequences for a stable and accurate solution of the governing system of equations. In this contribution, we follow the idea of NITSCHE [29] who developed a stable scheme for the solution of the Laplace problem taking weak boundary conditions into account. An extension to problems from linear elasticity shows an appropriate behavior with regard to numerical stability, accuracy and an adequate convergence behavior. NURBS are chosen as a high-order approximation basis to benefit from their smoothness and flexibility in the process of uniform model refinement

Observations of the J = 10 manifold of the pure rotational band of phosphine on Saturn

Saturn was observed in the vicinity of the J = 10 manifold of the pure rotational band of phosphine on 1984 July 10 and 12 from NASA's Kuiper Airborne Observatory with the facility far-infrared cooled grating spectrometer. On each night observations of the full disk plus rings were made at 4 to 6 discrete wavelengths which selectively sampled the manifold and the adjacent continuum. The previously reported detection of this manifold is confirmed. After subtraction of the flux due to the rings, the data are compared with disk-averaged models of Saturn. It is found that PH3 must be strongly depleted above the thermal inversion (approx. 70 mbar). The best fitting models consistent with other observational constaints indicate that PH3 is significantly depleted at even deeper atmospheric levels ( or = 500 mbar), implying an eddy diffusion coefficient for Saturn of 10 to the 4 cm sq/sec

Bayesian multiscale deconvolution applied to gamma-ray spectroscopy

A common task in gamma-ray astronomy is to extract spectral information, such as model constraints and incident photon spectrum estimates, given the measured energy deposited in a detector and the detector response. This is the classic problem of spectral “deconvolution” or spectral inversion. The methods of forward folding (i.e., parameter fitting) and maximum entropy “deconvolution” (i.e., estimating independent input photon rates for each individual energy bin) have been used successfully for gamma-ray solar flares (e.g., Rank, 1997; Share and Murphy, 1995). These methods have worked well under certain conditions but there are situations were they don’t apply. These are: 1) when no reasonable model (e.g., fewer parameters than data bins) is yet known, for forward folding; 2) when one expects a mixture of broad and narrow features (e.g., solar flares), for the maximum entropy method; and 3) low count rates and low signal-to-noise, for both. Low count rates are a problem because these methods (as they have been implemented) assume Gaussian statistics but Poisson are applicable. Background subtraction techniques often lead to negative count rates. For Poisson data the Maximum Likelihood Estimator (MLE) with a Poisson likelihood is appropriate. Without a regularization term, trying to estimate the “true” individual input photon rates per bin can be an ill-posed problem, even without including both broad and narrow features in the spectrum (i.e., amultiscale approach). One way to implement this regularization is through the use of a suitable Bayesian prior. Nowak and Kolaczyk (1999) have developed a fast, robust, technique using a Bayesian multiscale framework that addresses these problems with added algorithmic advantages. We outline this new approach and demonstrate its use with time resolved solar flare gamma-ray spectroscopy

On the natural stabilization of convection dominated problems using high order Bubnov–Galerkin finite elements

In the case of dominating convection, standard Bubnov–Galerkin finite elements are known to deliver oscillating discrete solutions for the convection–diffusion equation. This paper demonstrates that increasing the polynomial degree (p-extension) limits these artificial numerical oscillations. This is contrary to a widespread notion that an increase of the polynomial degree destabilizes the discrete solution. This treatise also provides explicit expressions as to which polynomial degree is sufficiently high to obtain stable solutions for a given Péclet number at the nodes of a mesh

Predicting Fracture in the Proximal Humerus using Phase Field Models

Proximal humerus impacted fractures are of clinical concern in the elderly population. Prediction of such fractures by CT-based finite element methods encounters several major obstacles such as heterogeneous mechanical properties and fracture due to compressive strains. We herein propose to investigate a variation of the phase field method (PFM) embedded into the finite cell method (FCM) to simulate impacted humeral fractures in fresh frozen human humeri. The force-strain response, failure loads and the fracture path are compared to experimental observations for validation purposes. The PFM (by means of the regularization parameter $l_0$) is first calibrated by one experiment and thereafter used for the prediction of the mechanical response of two other human fresh frozen humeri. All humeri are fractured at the surgical neck and strains are monitored by Digital Image Correlation (DIC). Experimental strains in the elastic regime are reproduced with good agreement ($R^2 = 0.726$), similarly to the validated finite element method [9]. The failure pattern and fracture evolution at the surgical neck predicted by the PFM mimic extremely well the experimental observations for all three humeri. The maximum relative error in the computed failure loads is $3.8\%$. To the best of our knowledge this is the first method that can predict well the experimental compressive failure pattern as well as the force-strain relationship in proximal humerus fractures

Are There Hidden Supernovae?

Ames Research Center and UCSC have been working on the development of a Mid IR Camera for the KAO in order to search for extra galactic supernovae. The development of the camera and its associated data reduction software have been successfully completed. Spectral Imaging of the Orion Bar at 6.2 and 7.8 microns demonstrates the derotation and data reduction software which was developed