基于介电泳的电极阵列电场仿真研究

介电泳方法被广泛地应用于微纳颗粒的分离和操纵中,实现介电泳操作的关键是设计满足所需电场分布的电极阵列。针对目前在微电极阵列设计中尚缺乏简单有效的电场解析方法的现状,提出一种基于格林公式的电极阵列电场的解析方法。首先介绍了传统介电泳和行波介电泳的概念和计算模型,分析了介电泳过程与电极上所施加的交变电压的频率和幅度的关系,然后在确立电极电势的边界条件的基础上,采用基于格林公式的电场解析方法,建立了非均匀电场的解析模型,得出不同条件下的电极阵列电场分布的仿真结果,最后利用FEMLAB有限元仿真软件对解析模型进行了对比仿真,验证了该解析模型的可行性。基于格林公式的电场解析求解方法能够有效地提高电极阵列设计中的针对性以及缩短电极设计的时间。

消除采集脚印的处理方法研究

Seismic technique is in the leading position for discovering oil and gas trap and searching for reserves throughout the course of oil and gas exploration. It needs high quality of seismic processed data, not only required exact spatial position, but also the true information of amplitude and AVO attribute and velocity. Acquisition footprint has an impact on highly precision and best quality of imaging and analysis of AVO attribute and velocity. Acquisition footprint is a new conception of describing seismic noise in 3-D exploration. It is not easy to understand the acquisition footprint. This paper begins with forward modeling seismic data from the simple sound wave model, then processes it and discusses the cause for producing the acquisition footprint. It agreed that the recording geometry is the main cause which leads to the distribution asymmetry of coverage and offset and azimuth in different grid cells. It summarizes the characters and description methods and analysis acquisition footprint’s influence on data geology interpretation and the analysis of seismic attribute and velocity. The data reconstruct based on Fourier transform is the main method at present for non uniform data interpolation and extrapolate, but this method always is an inverse problem with bad condition. Tikhonov regularization strategy which includes a priori information on class of solution in search can reduce the computation difficulty duo to discrete kernel condition disadvantage and scarcity of the number of observations. The method is quiet statistical, which does not require the selection of regularization parameter; and hence it has appropriate inversion coefficient. The result of programming and tentat-ive calculation verifies the acquisition footprint can be removed through prestack data reconstruct. This paper applies migration to the processing method of removing the acquisition footprint. The fundamental principle and algorithms are surveyed, seismic traces are weighted according to the area which occupied by seismic trace in different source-receiver distances. Adopting grid method in stead of accounting the area of Voroni map can reduce difficulty of calculation the weight. The result of processing the model data and actual seismic demonstrate, incorporating a weighting scheme based on the relative area that is associated with each input trace with respect to its neighbors acts to minimize the artifacts caused by irregular acquisition geometry.

南海北部陆缘深水区构造－热演化与烃源岩成熟度史研究

Abstract In order to provide basic data for evaluation of the petroleum potential in the deep water area of the northern margin of the South China Sea (SCS), present-day thermal regime and basin tectonothermal evolution are reconstructed and the maturation history of the Cenozoic major source rocks in the study area is derived. The present-day geothermal regime in the deep water area of the northern margin of SCS is defined according to the geothermal gradient, thermal properties and heat flow data. Tectonic subsidence history is reconstructed based on borehole and seismic data, and accordingly the stretching episodes are determined from the subsidence pattern. Heat flow history in the deep water area of the northern margin of SCS is estimated on a finite time, laterally non-uniform and multi-episode stretching model. Maturation history of the main source rocks in the study area is estimated through EASYRo% kinetic model and thermal history, and the potential of petroleum in the deep water area of the northern margin of SCS is evaluated based on the data above. The results show that the present-day geothermal regime in the deep water area of the northern margin of SCS is characterized by “hot basin” with high geothermal gradient (39.1±7.4℃/km) and high heat flow (77.5±14.8 mW/m2), and that the Qiongdongnan Basin (QDNB) underwent three stretching episodes and consequently suffered three heating episodes (Eocene, Oligocene and Pliocene time) with highest paleo-heat flow of 65～90 mW/m2 at the end of the Pliocene, that the Pearl River Mouth Basin (PRMB) two stretching and two heating episodes (Eocene, Oligocene time) with highest paleo-heat flow of 60～70 mW/m2 at the end of the Oligocene, and that the source rocks matured drastically responding to the heating episodes. There are four hydrocarbon generation kitchens in the deep water area of the northern margin of SCS which are favor of its bright petroleum perspective. Tectonothermal analysis indicates that the present-day geothermal regime which is characterized with “hot basin” in the deep water area of the PRMB resulted mainly from the Cenozoic stretching as well as faulting and magmatic activities during the Neotectonic period, and that the Pliocene heating episode of the QDNB is coupled with the transition from sinistral to dextral gliding of the Red Rive fault, and that the deep water basins in the northern margin of SCS are typical of multiple rifting which caused multi-episode heating process.

复杂介质中地震波场数值模拟研究——基于广义正交多项式迭积微分算子法

In this paper, we propose a new numerical modeling method – Convolutional Forsyte Polynomial Differentiator (CFPD), aimed at simulating seismic wave propagation in complex media with high efficiency and accuracy individually owned by short-scheme finite differentiator and general convolutional polynomial method. By adjusting the operator length and optimizing the operator coefficient, both global and local informations can be easily incorporated into the wavefield which is important to invert the undersurface geological structure. The key issue in this paper is to introduce the convolutional differentiator based on Forsyte generalized orthogonal polynomial in mathematics into the spatial differentiation of the first velocity-stress equation. To match the high accuracy of the spatial differentiator, this method in the time coordinate adopts staggered grid finite difference instead of conventional finite difference to model seismic wave propagation in heterogeneous media. To attenuate the reflection artifacts caused by artificial boundary, Perfectly Matched Layer (PML) absorbing boundary is also being considered in the method to deal with boundary problem due to its advantage of automatically handling large-angle emission. The PML formula for acoustic equation and first-order velocity-stress equation are also derived in this paper. There is little difference to implement the PML boundary condition in all kind of wave equations, but in Biot media, special attenuation factors should be taken. Numerical results demonstrate that the PML boundary condition is better than Cerjan absorbing boundary condition which makes it more suitable to hand the artificial boundary reflection. Based on the theories of anisotropy, Biot two-phase media and viscous-elasticity, this paper constructs the constitutive relationship for viscous-elastic and two-phase media, and further derives the first-order velocity-stress equation for 3D viscous-elastic and two-phase media. Numerical modeling using CFPD method is carried out in the above-mentioned media. The results modeled in the viscous-elastic media and the anisotropic pore elastic media can better explain wave phenomena of the true earth media, and can also prove that CFPD is a useful numerical tool to study the wave propagation in complex media.

粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

复杂裂隙岩体统计断裂本构关系与变形参数研究

Rock mass is widely recognized as a kind of geologic body which consists of rock blocks and discontinuities. The deformation and failure of rock mass is not only determined by rock block，but also by discontinuity which is virtually more important. Mutual cutting and combination of discontinuities controlled mechanical property of rock mass. The complex cutting of discontinuities determine the intense anisotropy on mechanical property of rock mass，especially under the effect of ground stress. Engineering practice has show that the brittle failure of hard rock always occurs when its working stress is far lower than the yield strength and compressive strength，the failure always directly related to the fracture propagation of discontinuities. Fracture propagation of discontinuities is the virtue of hard rock’s failure. We can research the rock mass discontinuous mechanical properties precisely by the methods of statistical analysis of discontinuities and Fracture Mechanics. According to Superposition Principle in Fracture Mechanics，A Problem or C Problem could be chosen to research. Problem A mainly calculates the crack-tip stress field and displacement field on internal discontinuities by numerical method. Problem C calculate the crack-tip stress field and displacement field under the assumption of that the mainly rock mass stress field has been known. So the Problem C avoid the complex mutual interference of stress fields of discontinuities，which is called crack system problem in Fracture Mechanics. To solve Problem C, field test on stress field in the rock mass is needed. The linear Superposition of discontinuities strain energies are Scientific and Rational. The difference of Fracture Mechanics between rock mass and other materials can mostly expression as：other materials Fracture Mechanics mostly face the problem A，and can’t avoid multi-crack puzzle, while the Rock mass Fracture Mechanics answer to the Problem C. Problem C can avoid multi-discontinuities mutual interference puzzle via the ground stress test. On the basis of Problem C, Fracture Mechanics could be used conveniently in rock mass. The rock mass statistics fracture constitutive relations, which introduced in this article, are based on the Problem C and the Discontinuity Strain Energy linear superposition. This constitutive relation has several merits: first, it is physical constitutive relation rather than empirical; second, it is very fit to describe the rock mass anisotropy properties; third, it elaborates the exogenous factors such as ground stress. The rock mass statistics fracture constitutive relation is the available approach to answer to the physical, anisotropic and ground stress impacted rock mass problems. This article stand on the foundation of predecessor’s statistics fractures constitutive relation, and improved the discontinuity distributive function. This article had derived the limitation of negative exponential distribution in the course of regression analysis, and advocated to using the two parameter negative exponential distribution for instead. In order to solve the problems of two-dimension stability on engineering key cross-sectional view in rock mass, this article derived the rock mass planar flexibility tensor, and established rock mass two-dimension penetrate statistics fracture constitutive relation on the basis of penetrate fracture mechanics. Based on the crack tip plasticity research production of penetrate fracture, for example the Irwin plasticity equifinality crack, this article established the way to deal with the discontinuity stress singularity and plastic yielding problem at discontinuity tip. The research on deformation parameters is always the high light region of rock mass mechanics field. After the dam foundation excavation of XiaoWan hydroelectric power station, dam foundation rock mass upgrowthed a great deal of unload cracks, rock mass mechanical property gotten intricacy and strong anisotropy. The dam foundation rock mass mostly upgrowthed three group discontinuities： the decantation discontinuity, the steep pitch discontinuity, and the schistosity plane. Most of the discontinuities have got partial unload looseness. In accordance with ground stress field data, the dam foundation stress field greatly non-uniform, which felled under the great impaction of tectonic stress field, self-weight stress field, excavation geometric boundary condition, and excavation, unload. The discontinuity complexity and stress field heterogeneity, created the rock mass mechanical property of dam foundation intricacy and levity. The research on the rock mass mechanics, if not take every respected influencing factor into consideration as best as we can, major errors likely to be created. This article calculated the rock mass elastic modulus that after Xiao Wan hydroelectric power station dam foundation gutter excavation finished. The calculation region covered possession monolith of Xiao Wan concrete double-curvature arch dam. Different monolith were adopted the penetrate fracture statistics constitutive relation or bury fracture statistics constitutive relation selectively. Statistics fracture constitutive relation is fit for the intensity anisotropy and heterogeneity rock mass of Xiao Wan hydroelectric power station dam foundation. This article had contrastive analysis the statistics fracture constitutive relation result with the inclined plane load test actual measurement elastic modulus and RMR method estimated elastic modulus, and find that the three methods elastic modulus have got greatly comparability. So, the statistics fracture constitutive relations are qualified for trust. Generally speaking，this article had finished following works based on predecessors job: “Argumentation the C Problems of superposition principle in Fracture Mechanics, establish two-dimension penetrate statistics fracture constitutive relation of rock mass, argue the negative exponential distribution limitation and improve it, improve of the three-dimension berry statistics fracture constitutive relation of rock mass, discontinuity-tip plastic zone isoeffect calculation, calculate the rock mass elastic modulus on two-dimension cross-sectional view”. The whole research clue of this article inherited from the “statistics rock mass mechanics” of Wu Faquan（1992）.

各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.

网络成瘾者心率变异性频谱特征及人格相关因素的研究

The aim of this research is to explore heart rate variability frequency field characteristics and personality influential factors of internet addicts by experimentation and questionnaire. Two studies were carried out: study 1 was to explore the pathological and mental mechanism of internet addicts by heart rate variability physiological index and EPQ and internet addiction scale. Study 2 was to compare the personality and mental characteristics between internet addicts and Non-addicts. The testees were 30 internet addict schoolboys who were in-patient and 43 schoolboys who accorded with qualification from grade 2 in senior high school and sophomore. It is found that: 1、 Internet addicts have obviously lower HFNU than Non-addicts， but have obviously higher LFNU and LF/HF ratio than Non-addicts. Internet addicts have dysfunction in their sympathetic and parasympathetic system; 2、 Internet addicts and Non-addicts have no significant difference in their EPQ except their lying point, but the nervous characteristics of EPQ of internet addicts have influence on their equipoise of parasympathetic system, that is to say when the score of nervous characteristics of EPQ become higher, equipoise of parasympathetic system become worse and worse. However the EPQ personality characteristics of non-addicts have no influence on their sympathetic and parasympathetic system. 3、 The extent of internet addiction of the internet addicts is independent of their equipoise of parasympathetic system, but non-addicts use internet more time, their equipoise of parasympathetic system become lower, there is the significant difference in the neurophysiology between internet addicts and non-addicts. That is to say internet environment is safe to most adolescent, they can make use of internet environment accurately, but internet addicts, their physiological and psychological level has changed to a certain, need to be treated.； 4、 Serious internet game addicts have the metal characteristics of low social support, low purpose in life and low adventure; 5、 The objective support、support utilization of social support questionnaire and sensation seeking characteristics have prognosticative function for internet addiction degree; 6、 Serious internet game addicts have the metal characteristics of low social support, low lying particularity, they are inclined to self abandonment. Serious internet game addicts of low N characteristic have high sensation seeking characteristic, and at the same E personality foundation, they also seek new stimulus with higher intension. There is a prompt that we should pay attention to internet addicts’ personality so that obtain better curative effect for internet addiction therapy.

Amorphous Placement and Informed Diffusion for Timely Monitoring by Autonomous, Resource-Constrained, Mobile Sensors

Personal communication devices are increasingly equipped with sensors for passive monitoring of encounters and surroundings. We envision the emergence of services that enable a community of mobile users carrying such resource-limited devices to query such information at remote locations in the field in which they collectively roam. One approach to implement such a service is directed placement and retrieval (DPR), whereby readings/queries about a specific location are routed to a node responsible for that location. In a mobile, potentially sparse setting, where end-to-end paths are unavailable, DPR is not an attractive solution as it would require the use of delay-tolerant (flooding-based store-carry-forward) routing of both readings and queries, which is inappropriate for applications with data freshness constraints, and which is incompatible with stringent device power/memory constraints. Alternatively, we propose the use of amorphous placement and retrieval (APR), in which routing and field monitoring are integrated through the use of a cache management scheme coupled with an informed exchange of cached samples to diffuse sensory data throughout the network, in such a way that a query answer is likely to be found close to the query origin. We argue that knowledge of the distribution of query targets could be used effectively by an informed cache management policy to maximize the utility of collective storage of all devices. Using a simple analytical model, we show that the use of informed cache management is particularly important when the mobility model results in a non-uniform distribution of users over the field. We present results from extensive simulations which show that in sparsely-connected networks, APR is more cost-effective than DPR, that it provides extra resilience to node failure and packet losses, and that its use of informed cache management yields superior performance.

An automatable generic strategy for dynamic load balancing in parallel structured mesh CFD codes

Parallel computing is now widely used in numerical simulation, particularly for application codes based on finite difference and finite element methods. A popular and successful technique employed to parallelize such codes onto large distributed memory systems is to partition the mesh into sub-domains that are then allocated to processors. The code then executes in parallel, using the SPMD methodology, with message passing for inter-processor interactions. In order to improve the parallel efficiency of an imbalanced structured mesh CFD code, a new dynamic load balancing (DLB) strategy has been developed in which the processor partition range limits of just one of the partitioned dimensions uses non-coincidental limits, as opposed to coincidental limits. The ‘local’ partition limit change allows greater flexibility in obtaining a balanced load distribution, as the workload increase, or decrease, on a processor is no longer restricted by the ‘global’ (coincidental) limit change. The automatic implementation of this generic DLB strategy within an existing parallel code is presented in this chapter, along with some preliminary results.

Portability, predictability and performance for parallel computing: BSP in practice

We report on practical experience using the Oxford BSP Library to parallelize a large electromagnetic code, the British Aerospace finite-difference time-domain code EMMA T:FD3D. The Oxford BS Library is one of the first realizations of the Bulk Synchronous Parallel computational model to be targeted at numerically intensive scientific (typically Fortran) computing. The BAe EMMA code is one of the first large-scale applications to be parallelized using this library, and it is an important demonstration of the cost effectiveness of the BSP approach. We illustrate how BSP cost-modelling techniques can be used to predict and optimize performance for single-source programs across different parallel platforms. We provide predicted and observed performance figures for an industrial-strength, single-source parallel code for a variety of real parallel architectures: shared memory multiprocessors, workstation clusters and massively parallel platforms.

Multiphysics simulation of microwave curing in micro-electronics packaging applications

Purpose – This paper aims to present an open-ended microwave curing system for microelectronics components and a numerical analysis framework for virtual testing and prototyping of the system, enabling design of physical prototypes to be optimized, expediting the development process. Design/methodology/approach – An open-ended microwave oven system able to enhance the cure process for thermosetting polymer materials utilised in microelectronics applications is presented. The system is designed to be mounted on a precision placement machine enabling curing of individual components on a circuit board. The design of the system allows the heating pattern and heating rate to be carefully controlled optimising cure rate and cure quality. A multi-physics analysis approach has been adopted to form a numerical model capable of capturing the complex coupling that exists between physical processes. Electromagnetic analysis has been performed using a Yee finite-difference time-domain scheme, while an unstructured finite volume method has been utilized to perform thermophysical analysis. The two solvers are coupled using a sampling-based cross-mapping algorithm. Findings – The numerical results obtained demonstrate that the numerical model is able to obtain solutions for distribution of temperature, rate of cure, degree of cure and thermally induced stresses within an idealised polymer load heated by the proposed microwave system. Research limitations/implications – The work is limited by the absence of experimentally derived material property data and comparative experimental results. However, the model demonstrates that the proposed microwave system would seem to be a feasible method of expediting the cure rate of polymer materials. Originality/value – The findings of this paper will help to provide an understanding of the behaviour of thermosetting polymer materials during microwave cure processing.

Fracture mechanics analysis of cracks in solder joint intermetallic compounds

The trend towards miniaturization of electronic products leads to the need for very small sized solder joints. Therefore, there is a higher reliability risk that too large a fraction of solder joints will transform into Intermetallic Compounds (IMCs) at the solder interface. In this paper, fracture mechanics study of the IMC layer for SnPb and Pb-free solder joints was carried out using finite element numerical computer modelling method. It is assumed that only one crack is present in the IMC layer. Linear Elastic Fracture Mechanics (LEFM) approach is used for parametric study of the Stress Intensity Factors (SIF, KI and KII), at the predefined crack in the IMC layer of solder butt joint tensile sample. Contrary to intuition, it is revealed that a thicker IMC layer in fact increases the reliability of solder joint for a cracked IMC. Value of KI and KII are found to decrease with the location of the crack further away from the solder interfaces while other parameters are constant. Solder thickness and strain rate were also found to have a significant influence on the SIF values. It has been found that soft solder matrix generates non-uniform plastic deformation across the solder-IMC interface near the crack tip that is responsible to obtain higher KI and KII.

Determining surface tension using DC magnetic levitation

The values of material physical properties are vital for the successful use of numerical simulations for electromagnetic processing of materials. The surface tension of materials can be determined from the experimental measurement of the surface oscillation frequency of liquid droplets. In order for this technique to be used, a positioning field is required that results in a modification to the oscillation frequency. A number of previous analytical models have been developed that mainly focus on electrically conducting droplets positioned using an A.C. electromagnetic field, but due to the turbulent flow resulting from the high electromagnetic fields required to balance gravity, reliable measurements have largely been limited to microgravity. In this work axisymmetric analytical and numerical models are developed, which allow the surface tension of a diamagnetic droplet positioned in a high DC magnetic field to be determined from the surface oscillations. In the case of D.C. levitation there is no internal electric currents with resulting Joule heating, Marangoni flow and other effects that introduce additional physics that complicates the measurement process. The analytical solution uses the linearised Navier-Stokes equations in the inviscid case. The body force from a DC field is potential, in contrast to the AC case, and it can be derived from Maxwell equations giving a solution for the magnetic field in the form of a series expansion of Legendre polynomials. The first few terms in this expansion represent a constant and gradient magnetic field valid close to the origin, which can be used to position the droplet. Initially the mathematical model is verified in microgravity conditions using a numerical model developed to solve the transient electromagnetics, fluid flow and thermodynamic equations. In the numerical model (as in experiment) the magnetic field is obtained using electrical current carrying coils, which provides the confinement force for a liquid droplet. The model incorporates free surface deformation to accurately model the oscillations that result from the interaction between the droplet and the non-uniform external magnetic field. A comparison is made between the analytical perturbation theory and the numerical pseudo spectral approximation solutions for small amplitude oscillations.

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.