Modelling technology to predict flip-chip assembly

This paper describes modelling technology and its use in providing data governing the assembly of flip-chip components. Details are given on the reflow and curing stages as well as the prediction of solder joint shapes. The reflow process involves the attachment of a die to a board via solder joints. After a reflow process, underfill material is placed between the die and the substrate where it is heated and cured. Upon cooling the thermal mismatch between the die, underfill, solder bumps, and substrate will result in a nonuniform deformation profile across the assembly and hence stress. Shape predictions then thermal solidification and stress prediction are undertaken on solder joints during the reflow process. Both thermal and stress calculations are undertaken to predict phenomena occurring during the curing of the underfill material. These stresses may result in delamination between the underfill and its surrounding materials leading to a subsequent reduction in component performance and lifetime. Comparisons between simulations and experiments for die curvature will be given for the reflow and curing process

CFD approach for solid mechanics analysis

A cell-centred finite volume(CC-FV) solid mechanics formulation, based on a computational fluid dynamics(CFD) procedure, is presented. A CFD code is modified such that the velocity variable is used as to the displacement variable. Displacement and pressure fields are considered as unknown variables. The results are validated with finite element(FE) and cell-vertex finite volume(CV-FV) predictions based on discretisation of the equilibrium equations. The developed formulation is applicable for both compressible and incompressible solids behaviour. The method is general and can be extended for the simultaneous analysis of problems involving flow-thermal and stress effects.

Finite element modelling of crack detection tests

Four non-destructive tests for determining the length of fatigue cracks within the solder joints of a 2512 surface mount resistor are investigated. The sensitivity of the tests is obtained using finite element analysis with some experimental validation. Three of the tests are mechanically based and one is thermally based. The mechanical tests all operate by applying different loads to the PCB and monitoring the strain response at the top of the resistor. The thermal test operates by applying a heat source underneath the PCB, and monitoring the temperature response at the top of the resistor. From the modelling work done, two of these tests have shown to be sensitive to cracks. Some experimental results are presented but further work is required to fully validate the simulation results.

Mathematical modelling: a laser soldering process for an optoelectronics butterfly package

For sensitive optoelectronic components, traditional soldering techniques cannot be used because of their inherent sensitivity to thermal stresses. One such component is the Optoelectronic Butterfly Package which houses a laser diode chip aligned to a fibre-optic cable. Even sub-micron misalignment of the fibre optic and laser diode chip can significantly reduce the performance of the device. The high cost of each unit requires that the number of damaged components, via the laser soldering process, are kept to a minimum. Mathematical modelling is undertaken to better understand the laser soldering process and to optimize operational parameters such as solder paste volume, copper pad dimensions, laser solder times for each joint, laser intensity and absorption coefficient. Validation of the model against experimental data will be completed, and will lead to an optimization of the assembly process, through an iterative modelling cycle. This will ultimately reduce costs, improve the process development time and increase consistency in the laser soldering process.

Microwave cure of conductive adhesives for flip-chip & microsystems applications

The curing of conductive adhesives and underfills can save considerable time and offer cost benefits for the microsystems and electronics packaging industry. In contrast to conventional ovens, curing by microwave energy generates heat internally within each individual component of an assembly. The rate at which heat is generated is different for each of the components and depends on the material properties as well as the oven power and frequency. This leads to a very complex and transient thermal state, which is extremely difficult to measure experimentally. Conductive adhesives need to be raised to a minimum temperature to initiate the cross-linking of the resin polymers, whilst some advanced packaging materials currently under investigation impose a maximum temperature constraint to avoid damage. Thermal imagery equipment integrated with the microwave oven can offer some information on the thermal state but such data is based on the surface temperatures. This paper describes computational models that can simulate the internal temperatures within each component of an assembly including the critical region between the chip and substrate. The results obtained demonstrate that due to the small mass of adhesive used in the joints, the temperatures reached are highly dependent on the material properties of the adjacent chip and substrate.

Computer modelling of the reliability of flip chips with metal column bumping

Recently, research has been carried out to test a novel bumping method which omits the under bump metallurgy forming process by bonding copper columns directly onto the Al pads of the silicon dies. This bumping method could be adopted to simplify the flip chip manufacturing process, increase the productivity and achieve a higher I/O count. This paper describes an investigation of the solder joint reliability of flip-chips based on this new bumping process. Computer modelling methods are used to predict the shape of solder joints and response of flip chips to thermal cyclic loading. The accumulated plastic strain energy at the comer solder joints is used as the damage indicator. Models with a range of design parameters have been compared for their reliability. The parameters that have been investigated are the copper column height, radius and solder volume. The ranking of the relative importance of these parameters is given. For most of the results presented in the paper, the solder material has been assumed to be the lead-free 96.5Sn3.5Ag alloy but some results for 60Sn40Pb solder joints have also been presented.

Dynamic model for metal cleanness evaluation by melting in a cold crucible

Melting of metallic samples in a cold crucible causes inclusions to concentrate on the surface owing to the action of the electromagnetic force in the skin layer. This process is dynamic, involving the melting stage, then quasi-stationary particle separation, and finally the solidification in the cold crucible. The proposed modeling technique is based on the pseudospectral solution method for coupled turbulent fluid flow, thermal and electromagnetic fields within the time varying fluid volume contained by the free surface, and partially the solid crucible wall. The model uses two methods for particle tracking: (1) a direct Lagrangian particle path computation and (2) a drifting concentration model. Lagrangian tracking is implemented for arbitrary unsteady flow. A specific numerical time integration scheme is implemented using implicit advancement that permits relatively large time-steps in the Lagrangian model. The drifting concentration model is based on a local equilibrium drift velocity assumption. Both methods are compared and demonstrated to give qualitatively similar results for stationary flow situations. The particular results presented are obtained for iron alloys. Small size particles of the order of 1 μm are shown to be less prone to separation by electromagnetic field action. In contrast, larger particles, 10 to 100 μm, are easily “trapped” by the electromagnetic field and stay on the sample surface at predetermined locations depending on their size and properties. The model allows optimization for melting power, geometry, and solidification rate.

Thermoelectric MHD effects on equiaxed crystal morphology

The effects of a constant uniform magnetic field on a growing equiaxed crystal are investigated using a 3-dimensional enthalpy based numerical model. Two cases are considered: The first case looks at unconstrained growth, where the current density is generated through the thermo-electric effect and the current circulates between the tips and roots of the dendrite, the second represents an imposed potential difference across the domain. A jump in the electrical conductivity between the liquid and solid causes the current density to be non uniform. In both cases the resulting Lorentz force drives fluid flow in the liquid phase, this in turn causes advection of the thermal and solute field altering the free energy close to the interface and changing the morphology of the dendrite. In the first case the flow field is complex comprising of many circulations, the morphological changes are modelled using a 2D model with a quasi 3D approximation. The second case is comparable to classic problems involving a constant velocity boundary.

Seasonal dynamics of the carbonate system in the Western English Channel

We present over 900 carbonate system observations collected over four years (2007–2010) in the Western English Channel (WEC). We determined CO2 partial pressure (pCO2), Total Alkalinity (TA) and Dissolved Inorganic Carbon (DIC) along a series of 40 km transects, including two oceanographic stations (L4 and E1) within a sustained coastal observatory. Our data follow a seasonal pattern of CO2 undersaturation from January to August, followed by supersaturation in September–October and a return to near-equilibrium thereafter. This pattern is explained by the interplay of thermal and biological sinks in winter and spring–summer, respectively, followed by the breakdown of stratification and mixing with deeper, high-CO2 water in autumn. The drawdown of DIC and inorganic N between March and June with a C:N ratio of 8.7–9.5 was consistent with carbon over-consumption during phytoplankton growth. Monthly mean surface pCO2 was strongly correlated with depth integrated chlorophyll a highlighting the importance of subsurface chlorophyll a maxima in controlling C-fluxes in shelf seas. Mixing of seawater with riverine freshwater in near-shore samples caused a reduction in TA and the saturation state of calcite minerals, particularly in winter. Our data show that the L4 and E1 oceanographic stations were small, net sinks for atmospheric CO2 over an annual cycle (−0.52±0.66 mol C m−2 y−1 and −0.62±0.49 mol C m−2 y−1, respectively).

The Impact of Offshore Wind Farms on Marine Ecosystems: A Review Taking an Ecosystem Services Perspective

The rapid increase in renewable energy generation from wind has increased concerns about the impacts that wind arrays have on the marine environment and what these impacts mean for society. One method for identifying the impacts of offshore wind farms (OWFs) on human welfare is through the assessment and valuation of ecosystem services. Using an ecosystem services approach, this paper reviews the impacts of OWFs on the ecosystem services delivered by marine environments. During the construction phase, supporting services such as reduced energy capture and nutrient cycling are changed due to the introduction of hard substrate and the reduction in soft sediment habitat at turbine bases. This may lead to changes in all other ecosystem services, both negative and positive. Quantifying these changes, however, is a challenge partly due to data limitations and a lack of clear understanding of the impacts of OWFs on the marine ecosystems. Scientific effort needs to quantitatively explore the impacts of OWFs on ecosystem functionality and the gathering of data that enables the assessment of changes to ecosystem services. Data needed to better quantify and value the impacts of OWFs on ecosystem services are suggested. The development of methods which integrate socioeconomic valuation of ecosystem services into the evaluation of renewable energy devices compliments efforts in assessing the environmental impacts and should enable a holistic assessment of the impact of renewable energy production and greenhouse gas mitigation technologies on the U. K. carbon footprint.

Multiple environmental controls on phytoplankton growth strategies determine adaptive responses of the N : P ratio

The controls on the 'Redfield' N:P stoichiometry of marine phytoplankton and hence the N:P ratio of the deep ocean remain incompletely understood. Here, we use a model for phytoplankton ecophysiology and growth, based on functional traits and resource-allocation trade-offs, to show how environmental filtering, biotic interactions, and element cycling in a global ecosystem model determine phytoplankton biogeography, growth strategies and macromolecular composition. Emergent growth strategies capture major observed patterns in marine biomes. Using a new synthesis of experimental RNA and protein measurements to constrain per-ribosome translation rates, we determine a spatially variable lower limit on adaptive rRNA:protein allocation and hence on the relationship between the largest cellular P and N pools. Comparison with the lowest observed phytoplankton N:P ratios and N:P export fluxes in the Southern Ocean suggests that additional contributions from phospholipid and phosphorus storage compounds play a fundamental role in determining the marine biogeochemical cycling of these elements.

Mesoscale fronts as foraging habitats: composite front mapping reveals oceanographic drivers of habitat use for a pelagic seabird

The oceanographic drivers of marine vertebrate habitat use are poorly understood yet fundamental to our knowledge of marine ecosystem functioning. Here, we use composite front mapping and high-resolution GPS tracking to determine the significance of mesoscale oceanographic fronts as physical drivers of foraging habitat selection in northern gannets Morus bassanus. We tracked 66 breeding gannets from a Celtic Sea colony over 2 years and used residence time to identify area-restricted search (ARS) behaviour. Composite front maps identified thermal and chlorophyll-a mesoscale fronts at two different temporal scales—(i) contemporaneous fronts and (ii) seasonally persistent frontal zones. Using generalized additive models (GAMs), with generalized estimating equations (GEE-GAMs) to account for serial autocorrelation in tracking data, we found that gannets do not adjust their behaviour in response to contemporaneous fronts. However, ARS was more likely to occur within spatially predictable, seasonally persistent frontal zones (GAMs). Our results provide proof of concept that composite front mapping is a useful tool for studying the influence of oceanographic features on animal movements. Moreover, we highlight that frontal persistence is a crucial element of the formation of pelagic foraging hotspots for mobile marine vertebrates.

XoxF encoding an alternative methanol dehydrogenase is widespread in coastal marine environments

The xoxF gene, encoding a pyrroloquinoline quinone-dependent methanol dehydrogenase, is found in all known proteobacterial methylotrophs. In several newly discovered methylotrophs, XoxF is the active methanol dehydrogenase, catalysing the oxidation of methanol to formaldehyde. Apart from that, its potential role in methylotrophy and carbon cycling is unknown. So far, the diversity of xoxF in the environment has received little attention. We designed PCR primer sets targeting clades of the xoxF gene, and used 454 pyrosequencing of PCR amplicons obtained from DNA of four coastal marine environments for a unique assessment of the diversity of xoxF in these habitats. Phylogenetic analysis of the data obtained revealed a high diversity of xoxF genes from two of the investigated clades, and substantial differences in sequence composition between environments. Sequences were classified as being related to a wide range of both methylotrophs and non-methylotrophs from Alpha-, Beta- and Gammaproteobacteria. The most prominent sequences detected were related to the family Rhodobacteraceae, the genus Methylotenera and the OM43 clade of Methylophilales, and are thus related to organisms that employ XoxF for methanol oxidation. Furthermore, our analyses revealed a high degree of so far undescribed sequences, suggesting a high number of unknown species in these habitats.

Conceptual Ecological Modelling of Sublittoral Rock Habitats to Inform Indicator Selection

The purpose of this study is to produce a series of Conceptual Ecological Models (CEMs) that represent sublittoral rock habitats in the UK. CEMs are diagrammatic representations of the influences and processes that occur within an ecosystem. They can be used to identify critical aspects of an ecosystem that may be studied further, or serve as the basis for the selection of indicators for environmental monitoring purposes. The models produced by this project are control diagrams, representing the unimpacted state of the environment free from anthropogenic pressures. It is intended that the models produced by this project will be used to guide indicator selection for the monitoring of this habitat in UK waters. CEMs may eventually be produced for a range of habitat types defined under the UK Marine Biodiversity Monitoring R&D Programme (UKMBMP), which, along with stressor models, are designed to show the interactions within impacted habitats, would form the basis of a robust method for indicator selection. This project builds on the work to develop CEMs for shallow sublittoral coarse sediment habitats (Alexander et al 2014). The project scope included those habitats defined as ‘sublittoral rock’. This definition includes those habitats that fall into the EUNIS Level 3 classifications A3.1 Atlantic and Mediterranean high energy infralittoral rock, A3.2 Atlantic and Mediterranean moderate energy infralittoral rock, A3.3 Atlantic and Mediterranean low energy infralittoral rock, A4.1 Atlantic and Mediterranean high energy circalittoral rock, A4.2 Atlantic and Mediterranean moderate energy circalittoral rock, and A4.3 Atlantic and Mediterranean low energy circalittoral rock as well as the constituent Level 4 and 5 biotopes that are relevant to UK waters. A species list of characterising fauna to be included within the scope of the models was identified using an iterative process to refine the full list of species found within the relevant Level 5 biotopes. A literature review was conducted using a pragmatic and iterative approach to gather evidence regarding species traits and information that would be used to inform the models and characterise the interactions that occur within the sublittoral rock habitat. All information gathered during the literature review was entered into a data logging pro-forma spreadsheet that accompanies this report. Wherever possible, attempts were made to collect information from UK-specific peer-reviewed studies, although other sources were used where necessary. All data gathered was subject to a detailed confidence assessment. Expert judgement by the project team was utilised to provide information for aspects of the models for which references could not be sourced within the project timeframe. A multivariate analysis approach was adopted to assess ecologically similar groups (based on ecological and life history traits) of fauna from the identified species to form the basis of the models. A model hierarchy was developed based on these ecological groups. One general control model was produced that indicated the high-level drivers, inputs, biological assemblages, ecosystem processes and outputs that occur in sublittoral rock habitats. In addition to this, seven detailed sub-models were produced, which each focussed on a particular ecological group of fauna within the habitat: ‘macroalgae’, ‘temporarily or permanently attached active filter feeders’, ‘temporarily or permanently attached passive filter feeders’, ‘bivalves, brachiopods and other encrusting filter feeders’, ‘tube building fauna’, ‘scavengers and predatory fauna’, and ‘non-predatory mobile fauna’. Each sub-model is accompanied by an associated confidence model that presents confidence in the links between each model component. The models are split into seven levels and take spatial and temporal scale into account through their design, as well as magnitude and direction of influence. The seven levels include regional to global drivers, water column processes, local inputs/processes at the seabed, habitat and biological assemblage, output processes, local ecosystem functions, and regional to global ecosystem functions. The models indicate that whilst the high level drivers that affect each ecological group are largely similar, the output processes performed by the biota and the resulting ecosystem functions vary both in number and importance between groups. Confidence within the models as a whole is generally high, reflecting the level of information gathered during the literature review. Physical drivers which influence the ecosystem were found to be of high importance for the sublittoral rock habitat, with factors such as wave exposure, water depth and water currents noted to be crucial in defining the biological assemblages. Other important factors such as recruitment/propagule supply, and those which affect primary production, such as suspended sediments, light attenuation and water chemistry and temperature, were also noted to be key and act to influence the food sources consumed by the biological assemblages of the habitat, and the biological assemblages themselves. Output processes performed by the biological assemblages are variable between ecological groups depending on the specific flora and fauna present and the role they perform within the ecosystem. Of particular importance are the outputs performed by the macroalgae group, which are diverse in nature and exert influence over other ecological groups in the habitat. Important output processes from the habitat as a whole include primary and secondary production, bioengineering, biodeposition (in mixed sediment habitats) and the supply of propagules; these in turn influence ecosystem functions at the local scale such as nutrient and biogeochemical cycling, supply of food resources, sediment stability (in mixed sediment habitats), habitat provision and population and algae control. The export of biodiversity and organic matter, biodiversity enhancement and biotope stability are the resulting ecosystem functions that occur at the regional to global scale. Features within the models that are most useful for monitoring habitat status and change due to natural variation have been identified, as have those that may be useful for monitoring to identify anthropogenic causes of change within the ecosystem. Biological, physical and chemical features of the ecosystem have been identified as potential indicators to monitor natural variation, whereas biological factors and those physical /chemical factors most likely to affect primary production have predominantly been identified as most likely to indicate change due to anthropogenic pressures.

Control estructural de agrietamientos poligonales en la Pedriza del Manzanares (Madrid)

Los agrietamientos poligonales representan una forma común de modelado granítico cuyo origen y evolución continúa en fase de estudio, no existiendo una sistematización de estas estructuras diversas. Algunos autores explican su origen por procesos geodinámicos internos, relacionándolo con movimientos de planos de fractura en estados tardíos de consolidación magmática. Otros autores atribuyen su formación y desarrollo a factores externos relacionados con el régimen climático. La gran variedad de agrietamientos poligonales requiere la utilización de un número mayor de variables para definir los distintos orígenes, y las posibles interrelaciones entre los factores externos e internos, así como para explicar la evolución de dichas estructuras y avanzar en la clasificación de los patrones concretos. Este trabajo pretende contribuir a sistematizar los mecanismos que intervienen en el desarrollo de agrietamientos poligonales. Para ello se estudian únicamente agrietamientos poligonales formados sobre planos de fractura verticales o subverticales. En particular se establecen relaciones entre la presencia de agrietamientos poligonales y la red de fracturación, la altura de aparición, la orientación e inclinación de la pared, la morfología de las placas y la profundidad de incisión de las grietas perimetrales. Por otra parte, establece relaciones entre procesos geodinámicos internos y procesos de meteorización externos.