潮波模型的优化开边界方法研究

数值模式是潮波研究的一种有利手段，但在研究中会面临各种具体问题，包括开边界条件的确定、底摩擦系数和耗散系数的选取等。数据同化是解决这些问题的一种途径，即利用有限数量的潮汐观测资料对潮波进行最优估计，其根本目的是迫使模型预报值逼近观测值，使模式不要偏离实际情况太远。本文采用了一种优化开边界方法，沿着数值模型的开边界优化潮汐水位信息，目的是设法使数值解在动力约束的意义下接近观测值，获得研究区域的潮汐结果。边界值由指定优化问题的解来定，以提高模拟区域的潮汐精度，最优问题的解是基于通过开边界的能量通量的变化，处理开边界处的观测值与计算值之差的最小化。这里提供了辐射型边界条件，由Reid 和Bodine（本文简称为RB）推导，我们将采用的优化后的RB方法（称为ORB）是优化开边界的特殊情况。 本文对理想矩形海域（ E- E， N- N， 分辨率 ）进行了潮波模拟，有东部开边界，模式采用ECOM3D模式。对数据结果的误差分析采用，振幅平均偏差，平均绝对偏差，平均相对误差和均方根偏差四个值来衡量模拟结果的好坏程度。 需要优化入开边界的解析潮汐值本文采用的解析解由方国洪《海湾的潮汐与潮流》（1966年）方法提供，为验证本文所做的解析解和方文的一致，本文做了其第一个例子的关键值a,b,z,结果与其结果吻合的相当好。但略有差别，分析的可能原因是两法在具体迭代方案和计算机保留小数上有区别造成微小误差。另外，我们取m=20，得到更精确的数值，我们发现对前十项的各项参数值，取m=10,m=20各项参数略有改进。当然我们可以获得m更大的各项参数值。 同时为了检验解析解的正确性讨论m和l变化对边界值的影响，结果指出，增大m，m=20时，u的模最大在本身u1或u2的模的6%；m=100时，u的模最大在本身u1或u2的模的4%；m再增大，m=1000时，u的模最大在本身u1或u2的模的4%，改变不大。当l<1时， =0处u的模最大为2。当l=1时， =0处u的模最大为0.1，当l>1时，l越大，u的模越小，当l=10时，u的模最大为0.001，可以认为为0。 为检验该优化方法的应用情况，我们对理想矩形区域进行模拟，首先将本文所采用的优化开边界方法应用于30m的情况，在开边界优化入开边界得出模式解，所得模拟结果与解析解吻合得相当好，该模式解和解析解在整个区域上，振幅平均绝对偏差为9.9cm,相位平均绝对偏差只有4.0 ，均方根偏差只有13.3cm，说明该优化方法在潮波模型中有效。 为验证该优化方法在各种条件下的模拟结果情况，在下面我们做了三类敏感性试验： 第一类试验：为证明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，我们来比较ORB条件与RB条件的优劣，我们模拟用了两个不同的摩擦系数，k分别为：0,0.00006。 结果显示，针对不同摩擦系数，显示在开边界上使用ORB条件的解比使用RB条件的解无论是振幅还是相位都有显著改善，两个试验均方根偏差优化程度分别为84.3%，83.7%。说明在开边界上使用优化方法相比于没有采用优化方法的模拟解更接近于解析解，大大提高了模拟水平。上述的两个试验得出， k=0.00006优化结果比k=0的好。 第二类试验，使用ORB条件确定优化开边界情况下，在东西边界加入出入流的情况，流考虑线性和非线性情况，结果显示，加入流的情况，潮汐模拟的效果降低不少，流为1Sv的情况要比5Sv的情况均方根偏差相差20cm,而不加流的情况只有0.2cm。线性流和非线性流情况两者模式解相差不大，振幅，相位各项指数都相近, 说明流的线性与否对结果影响不大。 第三类试验，不仅在开边界使用ORB条件，在模式内部也使用ORB条件，比较了内部优化和不优化情况与解析解的偏差。结果显示，选用不同的k,振幅都能得到很好的模拟，而相位相对较差。另外，在内部优化的情况下，考虑不同的k的模式解， 我们选用了与解析解相近的6个模式解的k,结果显示，不同的k,振幅都能得到很好的模拟，而相位较差。 总之，在开边界使用ORB条件比使用RB条件好，振幅相位都有大幅度改进，在加入出入流情况下，流的大小对模拟结果有影响，但线形流和非线性流差别不大。内部优化的结果显示，模式采用不同的k都能很好模拟解析解的振幅。

Magnetic levitation fluid dynamics

This work is concerned with the accurate computation of flow in a rapidly deforming liquid metal droplet, suspended in an AC magnetic field. Intense flow motion due to the induced electromagnetic force distorts dynamically the droplet envelope, which is initially spherical. The relative positional change between the liquid metal surface and the surrounding coil means that fluid flow and magnetic field computations need to be closely coupled. A spectral technique is used to solve this problem, which is assumed axisymmetric. The computed results are compared against a physical experiment and "ideal sphere" analytic solutions. A comparison between the "magnetic pressure" approximation and the full electromagnetic force solutions, shows fundamental differences; the full electromagnetic force solution is necessary for accurate results in most practical applications of this technique. The physical reason for the fundamental discrepancy is the difference in the electromagnetic force representation: only the gradient part of the full force is accounted for in the "magnetic pressure" approximation. Figs 9, Refs 13.

Dynamics of magnetically suspended fluid

Magnetic suspension is a technique for processing pure or reactive materials without contact to walls. This work is concerned with the flow in the rapidly deforming liquid volume, suspended in an AC magnetic field. Intense flow motion due to the induced electromagnetic force distorts dynamically the droplet envelope. The relative positional change between the liquid surface and the surrounding coil means that fluid flow and magnetic field computations need to be closely coupled. The computed results are compared against a physical experiment and nearly spherical analytic solutions. A comparison between the "magetic pressure" approximation and the full electromagnetic force solutions shows fundamental differences; the full electromagnetic force is necessary for accurate results in most practical applications of this technique. The physical reason for the fundamental discrepancy is the difference in the electromagnetic force representation: only the gradient part of the full force is accounted for in the "magnetic pressure" approximation.

Superluminous supernovae from PESSTO

We present optical spectra and light curves for three hydrogen-poor superluminous supernovae followed by the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). Time series spectroscopy from a fewdays aftermaximum light to 100 d later shows them to be fairly typical of this class, with spectra dominated by Ca II, MgII, FeII, and Si II, which evolve slowly over most of the post-peak photospheric phase. We determine bolometric light curves and apply simple fitting tools, based on the diffusion of energy input by magnetar spin-down, Ni-56 decay, and collision of the ejecta with an opaque circumstellar shell. We investigate how the heterogeneous light curves of our sample (combined with others from the literature) can help to constrain the possible mechanisms behind these events. We have followed these events to beyond 100-200 d after peak, to disentangle host galaxy light from fading supernova flux and to differentiate between the models, which predict diverse behaviour at this phase. Models powered by radioactivity require unrealistic parameters to reproduce the observed light curves, as found by previous studies. Both magnetar heating and circumstellar interaction still appear to be viable candidates. A large diversity is emerging in observed tail-phase luminosities, with magnetar models failing in some cases to predict the rapid drop in flux. This would suggest either that magnetars are not responsible, or that the X-ray flux from the magnetar wind is not fully trapped. The light curve of one object shows a distinct rebrightening at around 100 d after maximum light. We argue that this could result either from multiple shells of circumstellar material, or from a magnetar ionization front breaking out of the ejecta.

Conventional and Reciprocal Approaches to the Forward and Inverse Problems of Electroencephalography

Le problème inverse en électroencéphalographie (EEG) est la localisation de sources de courant dans le cerveau utilisant les potentiels de surface sur le cuir chevelu générés par ces sources. Une solution inverse implique typiquement de multiples calculs de potentiels de surface sur le cuir chevelu, soit le problème direct en EEG. Pour résoudre le problème direct, des modèles sont requis à la fois pour la configuration de source sous-jacente, soit le modèle de source, et pour les tissues environnants, soit le modèle de la tête. Cette thèse traite deux approches bien distinctes pour la résolution du problème direct et inverse en EEG en utilisant la méthode des éléments de frontières (BEM): l’approche conventionnelle et l’approche réciproque. L’approche conventionnelle pour le problème direct comporte le calcul des potentiels de surface en partant de sources de courant dipolaires. D’un autre côté, l’approche réciproque détermine d’abord le champ électrique aux sites des sources dipolaires quand les électrodes de surfaces sont utilisées pour injecter et retirer un courant unitaire. Le produit scalaire de ce champ électrique avec les sources dipolaires donne ensuite les potentiels de surface. L’approche réciproque promet un nombre d’avantages par rapport à l’approche conventionnelle dont la possibilité d’augmenter la précision des potentiels de surface et de réduire les exigences informatiques pour les solutions inverses. Dans cette thèse, les équations BEM pour les approches conventionnelle et réciproque sont développées en utilisant une formulation courante, la méthode des résidus pondérés. La réalisation numérique des deux approches pour le problème direct est décrite pour un seul modèle de source dipolaire. Un modèle de tête de trois sphères concentriques pour lequel des solutions analytiques sont disponibles est utilisé. Les potentiels de surfaces sont calculés aux centroïdes ou aux sommets des éléments de discrétisation BEM utilisés. La performance des approches conventionnelle et réciproque pour le problème direct est évaluée pour des dipôles radiaux et tangentiels d’excentricité variable et deux valeurs très différentes pour la conductivité du crâne. On détermine ensuite si les avantages potentiels de l’approche réciproquesuggérés par les simulations du problème direct peuvent êtres exploités pour donner des solutions inverses plus précises. Des solutions inverses à un seul dipôle sont obtenues en utilisant la minimisation par méthode du simplexe pour à la fois l’approche conventionnelle et réciproque, chacun avec des versions aux centroïdes et aux sommets. Encore une fois, les simulations numériques sont effectuées sur un modèle à trois sphères concentriques pour des dipôles radiaux et tangentiels d’excentricité variable. La précision des solutions inverses des deux approches est comparée pour les deux conductivités différentes du crâne, et leurs sensibilités relatives aux erreurs de conductivité du crâne et au bruit sont évaluées. Tandis que l’approche conventionnelle aux sommets donne les solutions directes les plus précises pour une conductivité du crâne supposément plus réaliste, les deux approches, conventionnelle et réciproque, produisent de grandes erreurs dans les potentiels du cuir chevelu pour des dipôles très excentriques. Les approches réciproques produisent le moins de variations en précision des solutions directes pour différentes valeurs de conductivité du crâne. En termes de solutions inverses pour un seul dipôle, les approches conventionnelle et réciproque sont de précision semblable. Les erreurs de localisation sont petites, même pour des dipôles très excentriques qui produisent des grandes erreurs dans les potentiels du cuir chevelu, à cause de la nature non linéaire des solutions inverses pour un dipôle. Les deux approches se sont démontrées également robustes aux erreurs de conductivité du crâne quand du bruit est présent. Finalement, un modèle plus réaliste de la tête est obtenu en utilisant des images par resonace magnétique (IRM) à partir desquelles les surfaces du cuir chevelu, du crâne et du cerveau/liquide céphalorachidien (LCR) sont extraites. Les deux approches sont validées sur ce type de modèle en utilisant des véritables potentiels évoqués somatosensoriels enregistrés à la suite de stimulation du nerf médian chez des sujets sains. La précision des solutions inverses pour les approches conventionnelle et réciproque et leurs variantes, en les comparant à des sites anatomiques connus sur IRM, est encore une fois évaluée pour les deux conductivités différentes du crâne. Leurs avantages et inconvénients incluant leurs exigences informatiques sont également évalués. Encore une fois, les approches conventionnelle et réciproque produisent des petites erreurs de position dipolaire. En effet, les erreurs de position pour des solutions inverses à un seul dipôle sont robustes de manière inhérente au manque de précision dans les solutions directes, mais dépendent de l’activité superposée d’autres sources neurales. Contrairement aux attentes, les approches réciproques n’améliorent pas la précision des positions dipolaires comparativement aux approches conventionnelles. Cependant, des exigences informatiques réduites en temps et en espace sont les avantages principaux des approches réciproques. Ce type de localisation est potentiellement utile dans la planification d’interventions neurochirurgicales, par exemple, chez des patients souffrant d’épilepsie focale réfractaire qui ont souvent déjà fait un EEG et IRM.

Planning rigid body motions using elastic curves

This paper tackles the problem of computing smooth, optimal trajectories on the Euclidean group of motions SE(3). The problem is formulated as an optimal control problem where the cost function to be minimized is equal to the integral of the classical curvature squared. This problem is analogous to the elastic problem from differential geometry and thus the resulting rigid body motions will trace elastic curves. An application of the Maximum Principle to this optimal control problem shifts the emphasis to the language of symplectic geometry and to the associated Hamiltonian formalism. This results in a system of first order differential equations that yield coordinate free necessary conditions for optimality for these curves. From these necessary conditions we identify an integrable case and these particular set of curves are solved analytically. These analytic solutions provide interpolating curves between an initial given position and orientation and a desired position and orientation that would be useful in motion planning for systems such as robotic manipulators and autonomous-oriented vehicles.

A simple analysis of thermodynamic properties for classical plasmas: I. theory

By eliminating the short range negative divergence of the Debye–Hückel pair distribution function, but retaining the exponential charge screening known to operate at large interparticle separation, the thermodynamic properties of one-component plasmas of point ions or charged hard spheres can be well represented even in the strong coupling regime. Predicted electrostatic free energies agree within 5% of simulation data for typical Coulomb interactions up to a factor of 10 times the average kinetic energy. Here, this idea is extended to the general case of a uniform ionic mixture, comprising an arbitrary number of components, embedded in a rigid neutralizing background. The new theory is implemented in two ways: (i) by an unambiguous iterative algorithm that requires numerical methods and breaks the symmetry of cross correlation functions; and (ii) by invoking generalized matrix inverses that maintain symmetry and yield completely analytic solutions, but which are not uniquely determined. The extreme computational simplicity of the theory is attractive when considering applications to complex inhomogeneous fluids of charged particles.

Analysis of wave phenomena in a morphogenetic mechanochemical model and an application to post-fertilization waves on eggs

Wave solutions to a mechanochemical model for cytoskeletal activity are studied and the results applied to the waves of chemical and mechanical activity that sweep over an egg shortly after fertilization. The model takes into account the calcium-controlled presence of actively contractile units in the cytoplasm, and consists of a viscoelastic force equilibrium equation and a conservation equation for calcium. Using piecewise linear caricatures, we obtain analytic solutions for travelling waves on a strip and demonstrate uiat the full nonlinear system behaves as predicted by the analytic solutions. The equations are solved on a sphere and the numerical results are similar to the analytic solutions. We indicate how the speed of the waves can be used as a diagnostic tool with which the chemical reactivity of the egg surface can be measured.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.

Stability of numerical schemes on staggered grids

This paper considers the stability of explicit, implicit and Crank-Nicolson schemes for the one-dimensional heat equation on a staggered grid. Furthemore, we consider the cases when both explicit and implicit approximations of the boundary conditions arc employed. Why we choose to do this is clearly motivated and arises front solving fluid flow equations with free surfaces when the Reynolds number can be very small. in at least parts of the spatial domain. A comprehensive stability analysis is supplied: a novel result is the precise stability restriction on the Crank-Nicolson method when the boundary conditions are approximated explicitly, that is, at t =n delta t rather than t = (n + 1)delta t. The two-dimensional Navier-Stokes equations were then solved by a marker and cell approach for two simple problems that had analytic solutions. It was found that the stability results provided in this paper were qualitatively very similar. thereby providing insight as to why a Crank-Nicolson approximation of the momentum equations is only conditionally, stable. Copyright (C) 2008 John Wiley & Sons, Ltd.

Modelo Estocástico para bloqueio de poros e redução de permeabilidade

Modeling transport of particulate suspensions in porous media is essential for understanding various processes of industrial and scientific interest. During these processes, particles are retained due to mechanisms like size exclusion (straining), adsorption, sedimentation and diffusion. In this thesis, a mathematical model is proposed and analytical solutions are obtained. The obtained analytic solutions for the proposed model, which takes pore and particle size distributions into account, were applied to predict the particle retention, pore blocking and permeability reduction during dead-end microfiltration in membranes. Various scenarios, considering different particle and pore size distributions were studied. The obtained results showed that pore blocking and permeability reduction are highly influenced by the initial pore and particle size distributions. This feature was observed even when different initial pore and particle size distributions with the same average pore size and injected particle size were considered. Finally, a mathematical model for predicting equivalent permeability in porous media during particle retention (and pore blocking) is proposed and the obtained solutions were applied to study permeability decline in different scenarios

Modelagem matemática para o transporte de partículas sujeitas a múltiplos mecanismos de retenção

Discrepancies between classical model predictions and experimental data for deep bed filtration have been reported by various authors. In order to understand these discrepancies, an analytic continuum model for deep bed filtration is proposed. In this model, a filter coefficient is attributed to each distinct retention mechanism (straining, diffusion, gravity interception, etc.). It was shown that these coefficients generally cannot be merged into an effective filter coefficient, as considered in the classical model. Furthermore, the derived analytic solutions for the proposed model were applied for fitting experimental data, and a very good agreement between experimental data and proposed model predictions were obtained. Comparison of the obtained results with empirical correlations allowed identifying the dominant retention mechanisms. In addition, it was shown that the larger the ratio of particle to pore sizes, the more intensive the straining mechanism and the larger the discrepancies between experimental data and classical model predictions. The classical model and proposed model were compared via statistical analysis. The obtained p values allow concluding that the proposed model should be preferred especially when straining plays an important role. In addition, deep bed filtration with finite retention capacity was studied. This work also involves the study of filtration of particles through porous media with a finite capacity of filtration. It was observed, in this case, that is necessary to consider changes in the boundary conditions through time evolution. It was obtained a solution for such a model using different functions of filtration coefficients. Besides that, it was shown how to build a solution for any filtration coefficient. It was seen that, even considering the same filtration coefficient, the classic model and the one here propposed, show different predictions for the concentration of particles retained in the porous media and for the suspended particles at the exit of the media

A marker-and-cell approach to free surface 2-D multiphase flows

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Noise and ultraviolet divergences in the dynamics of the chiral condensate in QCD

The time evolution of the matter produced in high energy heavy-ion collisions seems to be well described by relativistic viscous hydrodynamics. In addition to the hydrodynamic degrees of freedom related to energy-momentum conservation, degrees of freedom associated with order parameters of broken continuous symmetries must be considered because they are all coupled to each other. of particular interest is the coupling of degrees of freedom associated with the chiral symmetry of QCD. Quantum and thermal fluctuations of the chiral fields act as noise sources in the classical equations of motion, turning them into stochastic differential equations in the form of Ginzburg-Landau-Langevin (GLL) equations. Analytic solutions of GLL equations are attainable only in very special circumstances and extensive numerical simulations are necessary, usually by discretizing the equations on a spatial lattice. However, a not much appreciated issue in the numerical simulations of GLL equations is that ultraviolet divergences in the form of lattice-spacing dependence plague the solutions. The divergences are related to the well-known Rayleigh-Jeans catastrophe in classical field theory. In the present communication we present a systematic lattice renormalization method to control the catastrophe. We discuss the implementation of the method for a GLL equation derived in the context of a model for the QCD chiral phase transition and consider the nonequilibrium evolution of the chiral condensate during the hydrodynamic flow of the quark-gluon plasma.

High-order upwinding and the hydraulic jump

This work is concerned with the computation of incompressible axisymmetric and fall three-dimensional free-surface flows. In particular, the circular-hydraulic jump is simulated and compared with approximate analytic solutions. However, the principal thrust of this paper is to provide a real problem as a test bed for comparing the many existing convective approximations. Their performance is compared; SMART, HLPA and VONOS emerge as acceptable upwinding methods for this problem. Copyright (C) 2002 John Wiley Sons, Ltd.