Assessing the scalability of multiphysics tools for modeling solidification and melting processes on parallel clusters

A comprehensive solution of solidification/melting processes requires the simultaneous representation of free surface fluid flow, heat transfer, phase change, nonlinear solid mechanics and, possibly, electromagnetics together with their interactions, in what is now known as multiphysics simulation. Such simulations are computationally intensive and the implementation of solution strategies for multiphysics calculations must embed their effective parallelization. For some years, together with our collaborators, we have been involved in the development of numerical software tools for multiphysics modeling on parallel cluster systems. This research has involved a combination of algorithmic procedures, parallel strategies and tools, plus the design of a computational modeling software environment and its deployment in a range of real world applications. One output from this research is the three-dimensional parallel multiphysics code, PHYSICA. In this paper we report on an assessment of its parallel scalability on a range of increasingly complex models drawn from actual industrial problems, on three contemporary parallel cluster systems.

Modeling the structural breakdown of solder paste using the structural kinetic model

Solder paste is the most important strategic bonding material used in the assembly of surface mount devices in electronic industries. It is known to exhibit a thixotropic behavior, which is recognized by the decrease in apparent viscosity of paste material with time when subjected to a constant shear rate. The proper characterization of this time-dependent rheological behavior of solder pastes is crucial for establishing the relationships between the pastes structure and flow behavior; and for correlating the physical parameters with paste printing performance. In this article, we present a novel method which has been developed for characterizing the time-dependent and non-Newtonian rheological behavior of solder pastes and flux mediums as a function of shear rates. We also present results of the study of the rheology of the solder pastes and flux mediums using the structural kinetic modeling approach, which postulates that the network structure of solder pastes breaks down irreversibly under shear, leading to time and shear-dependent changes in the flow properties. Our results show that for the solder pastes used in the study, the rate and extent of thixotropy was generally found to increase with increasing shear rate. The technique demonstrated in this study has wide utility for R&D personnel involved in new paste formulation, for implementing quality control procedures used in solder-paste manufacture and packaging; and for qualifying new flip-chip assembly lines.

Sun glint estimation in marine satellite images: a comparison of results from calculation and radiative transfer modeling

The intensity and location of Sun glint in two Medium Resolution Imaging Spectrometer (MERIS) images was modeled using a radiative transfer model that includes elevation features as well as the slope of the sea surface. The results are compared to estimates made using glint flagging and correction approaches used within standard atmospheric correction processing code. The model estimate gives a glint pattern with a similar width but lower peak level than any current method, or than that estimated by a radiative transfer model with surfaces that include slope but not height. The MERIS third reprocessing recently adopted a new slope statistics model for Sun glint correction; the results show that this model is an outlier with respect to both the elevation model and other slope statistics models and we recommend that its adoption should be reviewed.

Real-time sea surface air temperature retrieval over Southern Ocean using AMSR-E data

The AMSR-E satellite data and in-situ data were applied to retrieve sea surface air temperature (Ta) over the Southern Ocean. The in-situ data were obtained from the 24~(th) -26~(th) Chinese Antarctic Expeditions during 2008-2010. First, Ta was used to analyze the relativity with the bright temperature (Tb) from the twelve channels of AMSR-E, and no high relativity was found between Ta and Tb from any of the channels. The highest relativity was 0.38 (with 23.8 GHz). The dataset for the modeling was obtained by using in-situ data to match up with Tb, and two methods were applied to build the retrieval model. In multi-parameters regression method, the Tbs from 12 channels were used to the model and the region was divided into two parts according to the latitude of 50°S. The retrieval results were compared with the in-situ data. The Root Mean Square Error (RMS) and relativity of high latitude zone were 0.96℃and 0.93, respectively. And those of low latitude zone were 1.29 ℃ and 0.96, respectively. Artificial neural network (ANN) method was applied to retrieve Ta.The RMS and relativity were 1.26 ℃ and 0.98, respectively.

Chemistry and Release of Gases from the Surface Ocean

The sea-surface layer is the very upper part of the sea surface where reduced mixing leads to strong gradients in physical, chemical and biological properties1. This surface layer is naturally reactive, containing a complex chemistry of inorganic components and dissolved organic matter (DOM), the latter including amino acids, proteins, fatty acids, carbohydrates, and humic-type components,2 with a high proportion of functional groups such as carbonyls, carboxylic acids and aromatic moieties.3 The different physical and chemical properties of the surface of the ocean compared with bulk seawater, and its function as a gateway for molecules to enter the atmosphere or ocean phase, make this an interesting and important region for study. A number of chemical reactions are believed to occur on and in the surface ocean; these may be important or even dominant sources or sinks of climatically-active marine trace gases. However the sea surface, especially the top 1um to 1mm known as the sea surface microlayer (ssm), is critically under-sampled, so to date much of the evidence for such chemistry comes from laboratory and/or modeling studies. This review discusses the chemical and physical structure of the sea surface, mechanisms for gas transfer across it, and explains the current understanding of trace gas formation at this critical interface between the ocean and atmosphere.

A Virtual Inspection Technique for Assessing the Dimensional Accuracy of Forged Compressor Blades Using FE Modeling and CMM Inspection

This paper presents research for developing a virtual inspection system that evaluates the dimensional tolerance of forged aerofoil blades formed using the finite element (FE) method. Conventional algorithms adopted by modern coordinate measurement processes have been incorporated with the latest free-form surface evaluation techniques to provide a robust framework for the dimensional inspection of FE aerofoil models. The accuracy of the approach had been verified with a strong correlation obtained between the virtual inspection data and coordinate measurement data from corresponding aerofoil components.

Two Arginine-Glutamate Ionic Locks Near the Extracellular Surface of FFAR1 Gate Receptor Activation

Activation of a number of class A G protein-coupled receptors (GPCRs) is thought to involve two molecular switches, a rotamer toggle switch within the transmembrane domain and an ionic lock at the cytoplasmic surface of the receptor; however, the mechanism by which agonist binding changes these molecular interactions is not understood. Importantly, 80% of GPCRs including free fatty acid receptor 1 (FFAR1) lack the complement of amino acid residues implicated in either or both of these two switches; the mechanism of activation of these GPCRs is therefore less clear. By homology modeling, we identified two Glu residues (Glu-145 and Glu-172) in the second extracellular loop of FFAR1 that form putative interactions individually with two transmembrane Arg residues (Arg-183(5.39) and Arg-258(7.35)) to create two ionic locks. Molecular dynamics simulations showed that binding of agonists to FFAR1 leads to breakage of these Glu-Arg interactions. In mutagenesis experiments, breakage of these two putative interactions by substituting Ala for Glu-145 and Glu-172 caused constitutive receptor activation. Our results therefore reveal a molecular switch for receptor activation present on the extracellular surface of FFAR1 that is broken by agonist binding. Similar ionic locks between the transmembrane domains and the extracellular loops may constitute a mechanism common to other class A GPCRs also.

Modeling Frequency-Dependent Boundaries as Digital Impedance Filters in FDTD and K-DWM Room Acoustics Simulations

A new method for modeling-frequency-dependent boundaries in finite-difference time-domain (FDTD) and Kirchhoff variable digital waveguide mesh (K-DWM) room acoustics simulations is presented. The proposed approach allows the direct incorporation of a digital impedance filter (DIF) in the Multidimensional (2D or 3D) FDTD boundary model of a locally reacting surface. An explicit boundary update equation is obtained by carefully constructing a Suitable recursive formulation. The method is analyzed in terms of pressure wave reflectance for different wall impedance filters and angles of incidence. Results obtained from numerical experiments confirm the high accuracy of the proposed digital impedance filter boundary model, the reflectance of which matches locally reacting surface (LRS) theory closely. Furthermore a numerical boundary analysis (NBA) formula is provided as a technique for an analytic evaluation of the numerical reflectance of the proposed digital impedance filter boundary formulation.

Wavelength Dependence of Raman Enhancement from Gold Nanorod Arrays: Quantitative Experiment and Modeling of a Hot Spot Dominated System

Surface-enhanced Raman scattering (SERS) spectra from molecules adsorbed on the surface of vertically aligned gold nanorod arrays exhibit a variation in enhancement factor (EF) as a function of excitation wavelength that displays little correlation with the elastic optical properties of the surface. The key to understanding this lack of correlation and to obtaining agreement between experimental and calculated EF spectra lies with consideration of randomly distributed, sub-10 nm gaps between nanorods forming the substrate. Intense fields in these enhancement “hot spots” make a dominant contribution to the Raman scattering and have a very different spectral profile to that of the elastic optical response. Detailed modeling of the electric field enhancement at both excitation and scattering wavelengths was used to quantitatively predict both the spectral profile and the magnitude of the observed EF.

Finite difference time domain modeling of phase grating diffusion

In this paper, a method for modeling diffusion caused by non-smooth boundary surfaces in simulations of room acoustics using finite difference time domain (FDTD) technique is investigated. The proposed approach adopts the well-known theory of phase grating diffusers to efficiently model sound scattering from rough surfaces. The variation of diffuser well-depths is attained by nesting allpass filters within the reflection filters from which the digital impedance filters used in the boundary implementation are obtained. The presented technique is appropriate for modeling diffusion at high frequencies caused by small surface roughness and generally diffusers that have narrow wells and infinitely thin separators. The diffusion coefficient was measured with numerical experiments for a range of fractional Brownian diffusers.

Reaction Mechanisms of Crotonaldehyde Hydrogenation on Pt(111): Density Functional Theory and Microkinetic Modeling

The microkinetics based on density function theory (DFT) calculations is utilized to investigate the reaction mechanism of crotonaldehyde hydrogenation on Pt(111) in the free energy landscape. The dominant reaction channel of each hydrogenation product is identified. Each of them begins with the first surface hydrogenation of the carbonyl oxygen of crotonaldehyde on the surface. A new mechanism, 1,4-addition mechanism generating enols (butenol), which readily tautomerize to saturated aldehydes (butanal), is identified as a primary mechanism to yield saturated aldehydes instead of the 3,4-addition via direct hydrogenation of the ethylenic bond. The calculation results also show that the full hydrogenation product, butylalcohol, mainly stems from the deep hydrogenation of surface open-shell dihydrogenation intermediates. It is found that the apparent barriers of the dominant pathways to yield three final products are similar on P(111), which makes it difficult to achieve a high selectivity to the desired crotyl alcohol (COL).

In vitro modeling of respiratory syncytial virus infection of pediatric bronchial epithelium, the primary target of infection in vivo

Respiratory syncytial virus (RSV) is the major viral cause of severe pulmonary disease in young infants worldwide. However, the mechanisms by which RSV causes disease in humans remain poorly understood. To help bridge this gap, we developed an ex vivo/in vitro model of RSV infection based on well-differentiated primary pediatric bronchial epithelial cells (WD-PBECs), the primary targets of RSV infection in vivo. Our RSV/WD-PBEC model demonstrated remarkable similarities to hallmarks of RSV infection in infant lungs. These hallmarks included restriction of infection to noncontiguous or small clumps of apical ciliated and occasional nonciliated epithelial cells, apoptosis and sloughing of apical epithelial cells, occasional syncytium formation, goblet cell hyperplasia/metaplasia, and mucus hypersecretion. RSV was shed exclusively from the apical surface at titers consistent with those in airway aspirates from hospitalized infants. Furthermore, secretion of proinflammatory chemokines such as CXCL10, CCL5, IL-6, and CXCL8 reflected those chemokines present in airway aspirates. Interestingly, a recent RSV clinical isolate induced more cytopathogenesis than the prototypic A2 strain. Our findings indicate that this RSV/WD-PBEC model provides an authentic surrogate for RSV infection of airway epithelium in vivo. As such, this model may provide insights into RSV pathogenesis in humans that ultimately lead to successful RSV vaccines or therapeutics.

Modeling of subsurface pulsed radar for nondestructive testing of structures

The use of pulsed radar for investigating the integrity of structural elements is gaining popularity and becoming firmly established as a nondestructive test method in civil engineering. Difficulties can often arise in the interpretation of results obtained, particularly where internal details are relatively complex. One approach that can be used to understand and evaluate radar results is through numerical modeling of signal propagation and reflection. By comparing the results of a numerical modeling with those from field measurements, engineers can gain valuable insight into the probable features embedded beneath the surface of a structural element. This paper discusses a series of numerical techniques for modeling subsurface radar and compares the precision of the results with those taken from real field data. It is found that more complex problems require more sophisticated analysis techniques to obtain realistic results, with a consequential increase in the computational resources to carry out the modeling.