920 resultados para Higher order interior point method
Resumo:
This thesis consists of two parts. In Part I, we develop a multipole moment formalism in general relativity and use it to analyze the motion and precession of compact bodies. More specifically, the generic, vacuum, dynamical gravitational field of the exterior universe in the vicinity of a freely moving body is expanded in positive powers of the distance r away from the body's spatial origin (i.e., in the distance r from its timelike-geodesic world line). The expansion coefficients, called "external multipole moments,'' are defined covariantly in terms of the Riemann curvature tensor and its spatial derivatives evaluated on the body's central world line. In a carefully chosen class of de Donder coordinates, the expansion of the external field involves only integral powers of r ; no logarithmic terms occur. The expansion is used to derive higher-order corrections to previously known laws of motion and precession for black holes and other bodies. The resulting laws of motion and precession are expressed in terms of couplings of the time derivatives of the body's quadrupole and octopole moments to the external moments, i.e., to the external curvature and its gradient.
In part II, we study the interaction of magnetohydrodynamic (MHD) waves in a black-hole magnetosphere with the "dragging of inertial frames" effect of the hole's rotation - i.e., with the hole's "gravitomagnetic field." More specifically: we first rewrite the laws of perfect general relativistic magnetohydrodynamics (GRMHD) in 3+1 language in a general spacetime, in terms of quantities (magnetic field, flow velocity, ...) that would be measured by the ''fiducial observers” whose world lines are orthogonal to (arbitrarily chosen) hypersurfaces of constant time. We then specialize to a stationary spacetime and MHD flow with one arbitrary spatial symmetry (e.g., the stationary magnetosphere of a Kerr black hole); and for this spacetime we reduce the GRMHD equations to a set of algebraic equations. The general features of the resulting stationary, symmetric GRMHD magnetospheric solutions are discussed, including the Blandford-Znajek effect in which the gravitomagnetic field interacts with the magnetosphere to produce an outflowing jet. Then in a specific model spacetime with two spatial symmetries, which captures the key features of the Kerr geometry, we derive the GRMHD equations which govern weak, linealized perturbations of a stationary magnetosphere with outflowing jet. These perturbation equations are then Fourier analyzed in time t and in the symmetry coordinate x, and subsequently solved numerically. The numerical solutions describe the interaction of MHD waves with the gravitomagnetic field. It is found that, among other features, when an oscillatory external force is applied to the region of the magnetosphere where plasma (e+e-) is being created, the magnetosphere responds especially strongly at a particular, resonant, driving frequency. The resonant frequency is that for which the perturbations appear to be stationary (time independent) in the common rest frame of the freshly created plasma and the rotating magnetic field lines. The magnetosphere of a rotating black hole, when buffeted by nonaxisymmetric magnetic fields anchored in a surrounding accretion disk, might exhibit an analogous resonance. If so then the hole's outflowing jet might be modulated at resonant frequencies ω=(m/2) ΩH where m is an integer and ΩH is the hole's angular velocity.
Resumo:
把合成孔径激光成像雷达的目标衍射区分为三个区域,提出采用离焦或者附加空间相位调制板的光学接收望远镜补偿回波像差。当目标处于菲涅耳衍射区时可采用离焦或偏置望远镜补偿回波二次项离焦像差并产生用于孔径合成的二次项相位历程;目标处于夫琅和费衍射区时可以采用离焦或偏置望远镜补偿回波二次项离焦像差但不产生相位历程;目标处于瑞利-索末菲衍射区域时不可能补偿回波高阶像差。
Resumo:
We propose a united theory that describes the two-center recording system by taking scattering noise into account. The temporal evolution of the signal-to-noise ratio in doubly doped photorefractive crystals is described based on jointly solving material equations and coupled-wave equations with the fourth-order Runge-Kutta method. Roles of microcosmic optical parameters of dopants on the signal-to-noise ratio are discussed in detail. The theoretical results can confirm and predict experimental results. (c) 2005 Elsevier GmbH. All rights reserved.
Resumo:
O objetivo do presente trabalho foi estudar o comportamento dos potenciais superficiais e do perfil de potencial atraves da membrana de eritr ocito em func ao da forca i onica e das cargas superficiais, usando um modelo que leva em conta as cargas el etricas do glicoc alix e das proteınas citoplasm aticas, al em das cargas superficiais da bicamada lipıdica e os efeitos dos eletr olitos divalentes. Programas especıficos em linguagem C foram elaborados para o c alculo desses potenciais, tomando como dados num ericos resultados experimentais de medidas de mobilidade eletrofor etica de eritr ocitos para diferentes valores de forca i onica. Neste c alculo, o metodo para tratamento dos dados eletrofor eticos indicado por Hsu et al.[57] foi incluıdo em nosso modelo. A equac ao de Poisson-Boltzmann nao linear foi resolvida por computac ao num erica, usando o metodo de Runge-Kutta de quarta ordem, obtendo-se os perfis de potencial. Os resultados mostraram que a estimativa da densidade de carga el etrica na superfıcie de c elulas usando a equac ao cl assica de Helmholtz-Smoluchowski conduz a valores que nao conseguem refletir as forcas que regem o comportamento eletrofor etico das mesmas. O presente modelo gerou valores de potenciais superficiais e perfis de potencial para a membrana do eritr ocito bem distintos daqueles obtidos anteriormente para um modelo descrito por uma equac ao de Poisson-Boltzmann linear. Nossos resultados confirmam que a avaliac ao de parametros el etricos superficiais da membrana de eritr ocito, envolvendo dados oriundos de eletroforese, deve incluir c alculos hidrodin amicos al em de eletroest aticos, como sugerido por Hsu et al. [57].
Resumo:
The objective of this investigation has been a theoretical and experimental understanding of ferromagnetic resonance phenomena in ferromagnetic thin films, and a consequent understanding of several important physical properties of these films. Significant results have been obtained by ferromagnetic resonance, hysteresis, torque magnetometer, He ion backscattering, and X-ray fluorescence measurements for nickel-iron alloy films.
Taking into account all relevant magnetic fields, including the applied, demagnetizing, effective anisotropy and exchange fields, the spin wave resonance condition applicable to the thin film geometry is presented. On the basis of the simple exchange interaction model it is concluded that the normal resonance modes of an ideal film are expected to be unpinned. The possibility of nonideality near the surface of a real film was considered by means of surface anisotropy field, inhomogeneity in demagnetizing field and inhomogeneity of magnetization models. Numerical results obtained for reasonable parameters in all cases show that they negligibly perturb the resonance fields and the higher order mode shapes from those of the unpinned modes of ideal films for thicknesses greater than 1000 Å. On the other hand for films thinner than 1000 Å the resonance field deviations can be significant even though the modes are very nearly unpinned. A previously unnoticed but important feature of all three models is that the interpretation of the first resonance mode as the uniform mode of an ideal film allows an accurate measurement of the average effective demagnetizing field over the film volume. Furthermore, it is demonstrated that it is possible to choose parameters which give indistinguishable predictions for all three models, making it difficult to uniquely ascertain the source of spin pinning in real films from resonance measurements alone.
Spin wave resonance measurements of 81% Ni-19% Fe coevaporated films 30 to 9000 Å thick, at frequencies from 1 to 8 GHz, at room temperature, and with the static magnetic field parallel and perpendicular to the film plane have been performed. A self-consistent analysis of the results for films thicker than 1000 Å, in which multiple excitations can be observed, shows for the first time that a unique value of exchange constant A can only be obtained by the use of unpinned mode assignments. This evidence and the resonance behavior of films thinner than 1000 Å strongly imply that the magnetization at the surfaces of permalloy films is very weakly pinned. However, resonance measurements alone cannot determine whether this pinning is due to a surface anisotropy, an inhomogeneous demagnetizing field or an inhomogeneous magnetization. The above analysis yields a value of 4πM=10,100 Oe and A = (1.03 ± .05) x 10-6 erg/cm for this alloy. The ability to obtain a unique value of A suggests that spin wave resonance can be used to accurately characterize the exchange interaction in a ferromagnet.
In an effort to resolve the ambiguity of the source of pinning of the magnetization, a correlation of the ratio of magnetic moment and X-ray film thickness with the value of effective demagnetizing field 4πNM as determined from resonance, for films 45 to 300 Å has been performed. The remarkable agreement of both quantities and a comparison with the predictions of five distinct models, strongly imply that the thickness dependence of both quantities is related to a thickness dependent average saturation magnetization, which is far below 10,100 Oe for very thin films. However, a series of complementary experiments shows that this large decrease of average saturation magnetization cannot be simply explained by either oxidation or interdiffusion processes. It can only be satisfactorily explained by an intrinsic decrease of the average saturation magnetization for very thin films, an effect which cannot be justified by any simple physical considerations.
Recognizing that this decrease of average saturation magnetization could be due to an oxidation process, a correlation of resonance measurements, He ion backscattering, X-ray fluorescence and torque magnetometer measurements, for films 40 to 3500 Å thick has been performed. On basis of these measurements it is unambiguously established that the oxide layer on the surface of purposefully oxidized 81% Ni-19% Fe evaporated films is predominantly Fe-oxide, and that in the oxidation process Fe atoms are removed from the bulk of the film to depths of thousands of angstroms. Extrapolation of results for pure Fe films indicates that the oxide is most likely α-Fe2O3. These conclusions are in agreement with results from old metallurgical studies of high temperature oxidation of bulk Fe and Ni-Fe alloys. However, X-ray fluorescence results for films oxidized at room temperature, show that although the preferential oxidation of Fe also takes place in these films, the extent of this process is by far too small to explain the large variation of their average saturation magnetization with film thickness.
Resumo:
Organic contaminants are readily bioaccumulated by aquatic organisms. Exposure to and toxic effects of contaminants can be measured in terms of the biochemical responses of the organisms (i.e. molecular biomarkers). The hepatic biotransformation enzyme cytochrome P4501A (CYP1A) in vertebrates is specifically induced by organic contaminants such as aromatic hydrocarbons, PCBs and dioxins, and is involved in chemical carcinogenesis via catalysis of the covalent binding of organic contaminants to DNA (DNA-adducts). Hepatic CYP1A induction has been used extensively and successfully as a biomarker of organic contaminant exposure in fish. Fewer but equally encouraging studies in fish have used hepatic bulky, hydrophobic DNA-adducts as biomarkers of organic contaminant damage. Much less is known of the situation in marine invertebrates, but a CYPlA-like enzyme with limited inducibility and some potential for biomarker application is indicated. Stimulation of reactive oxygen species (ROS) production is another potential mechanism of organic contaminant-mediated DNA and other damage in aquatic organisms. A combination of antioxidant (enzymes, scavengers) and pro-oxidant (oxidised DNA bases, lipid peroxidation) measurements may have potential as a biomarker of organic contaminant exposure (particularly those chemicals which do not induce CYP1A) and/or oxidative stress, but more studies are required. Both CYP1A- and ROS-mediated toxicity are indicated to result in higher order deleterious effects, including cancer and other aspects of animal fitness.
Resumo:
O presente trabalho aborda um problema inverso associado a difus~ao de calor em uma barra unidimensional. Esse fen^omeno e modelado por meio da equac~ao diferencial par- cial parabolica ut = uxx, conhecida como equac~ao de difus~ao do calor. O problema classico (problema direto) envolve essa equac~ao e um conjunto de restric~oes { as condic~oes inicial e de contorno {, o que permite garantir a exist^encia de uma soluc~ao unica. No problema inverso que estudamos, o valor da temperatura em um dos extremos da barra n~ao esta disponvel. Entretanto, conhecemos o valor da temperatura em um ponto x0 xo no interior da barra. Para aproximar o valor da temperatura no intervalo a direita de x0, propomos e testamos tr^es algoritmos de diferencas nitas: diferencas regressivas, leap-frog e diferencas regressivas maquiadas.
Resumo:
In this paper, a refractive index pro. le design enabling us to obtain a. at modal field around the fibre centre is investigated. The theoretical approach for designing such multilayer large flattened mode (LFM) optical fibres is presented. A comparison is made between the properties of a three-layer LFM structure and a standard step-index pro. le with the same core size. The obtained results indicate that the effective area of the LFM fibre is about twice as large as that of the standard step-index fibre, but the LFM fibre has less effective ability to filter out the higher order modes than the standard step-index fibre with the same bending radius.
Resumo:
缠绕技术通常被用于大模场面积多模光纤激光器的横模控制。将光纤弯曲成不同半径的圆环,多模光纤激光器的高阶横模逐个被抑制并在缠绕半径为20mm时,获得15.4W的单模激光输出。实验研究表明,缠绕半径减小时,由于激光器高阶模式被抑制,其光束质量变好,同时其斜率效率降低。实验测量了多模光纤激光器在不同缠绕半径下的输出性能,并理论计算了各阶模式在不同缠绕半径下的损耗特性。实验测量结果与理论计算结果符合得较好。
Resumo:
O presente trabalho trata do escoamento bifásico em meios porosos heterogêneos de natureza fractal, onde os fluidos são considerados imiscíveis. Os meios porosos são modelados pela equação de Kozeny-Carman Generalizada (KCG), a qual relaciona a porosidade com a permeabilidade do meio através de uma nova lei de potência. Esta equação proposta por nós é capaz de generalizar diferentes modelos existentes na literatura e, portanto, é de uso mais geral. O simulador numérico desenvolvido aqui emprega métodos de diferenças finitas. A evolução temporal é baseada em um esquema de separação de operadores que segue a estratégia clássica chamada de IMPES. Assim, o campo de pressão é calculado implicitamente, enquanto que a equação da saturação da fase molhante é resolvida explicitamente em cada nível de tempo. O método de otimização denominado de DFSANE é utilizado para resolver a equação da pressão. Enfatizamos que o DFSANE nunca foi usado antes no contexto de simulação de reservatórios. Portanto, o seu uso aqui é sem precedentes. Para minimizar difusões numéricas, a equação da saturação é discretizada por um esquema do tipo "upwind", comumente empregado em simuladores numéricos para a recuperação de petróleo, o qual é resolvido explicitamente pelo método Runge-Kutta de quarta ordem. Os resultados das simulações são bastante satisfatórios. De fato, tais resultados mostram que o modelo KCG é capaz de gerar meios porosos heterogêneos, cujas características permitem a captura de fenômenos físicos que, geralmente, são de difícil acesso para muitos simuladores em diferenças finitas clássicas, como o chamado fenômeno de dedilhamento, que ocorre quando a razão de mobilidade (entre as fases fluidas) assume valores adversos. Em todas as simulações apresentadas aqui, consideramos que o problema imiscível é bidimensional, sendo, portanto, o meio poroso caracterizado por campos de permeabilidade e de porosidade definidos em regiões Euclideanas. No entanto, a teoria abordada neste trabalho não impõe restrições para sua aplicação aos problemas tridimensionais.
Resumo:
A tese descreve a ingestão de nutrientes segundo variáveis demográficas e socioeconômicas em adultos brasileiros, com base nos dados da primeira avaliação nacional do consumo alimentar individual, o Inquérito Nacional de Alimentação (INA), realizado entre 2008 e 2009. Um total de 34.003 indivíduos com pelo menos 10 anos de idade participaram do estudo. O presente estudo incluiu 21.003 indivíduos adultos, de 20 a 59 anos de idade, com exceção das mulheres gestantes e lactantes (n=1.065). O consumo alimentar individual foi estimado utilizando dois dias de registros alimentares não consecutivos. O consumo usual de nutrientes foi estimado pelo método do National Cancer Institute que permitiu a correção da variabilidade intraindividual. As prevalências de ingestão inadequada de nutrientes foram estimadas segundo o sexo e faixas etárias utilizando o método da necessidade média estimada como ponte de corte. A inadequação de sódio foi avaliada pelo consumo acima do nível de ingestão máximo tolerável. Os resultados são apresentados na forma de dois artigos. No primeiro artigo, estimaram-se as prevalências de inadequação segundo as cinco grandes regiões (Norte, Nordeste, Sudeste, Sul e Centro-Oeste) e a situação do domicílio (urbano e rural). Observaram-se prevalências de inadequação maiores ou iguais a 70% para cálcio entre os homens e magnésio, vitamina A, sódio em ambos os sexos. Prevalências maiores ou iguais a 90% foram encontradas para cálcio entre as mulheres e vitaminas D e E em ambos os sexos. No geral, os grupos com maior risco de inadequação de micronutrientes foram as mulheres e os que residem na área rural e na região Nordeste. No segundo artigo, estimaram-se as prevalências de inadequação do consumo segundo renda e escolaridade. A renda foi caracterizada pela renda mensal familiar per capita e a escolaridade definida pelo número de anos completos de estudo. Ambas variáveis foram categorizadas em quartis. Modelos de regressão linear simples e mutuamente ajustados foram estimados para verificar a associação independente entre o consumo de nutrientes e as variáveis socioeconômicas. Foram testadas as interações entre renda e escolaridade. Verificou-se que a inadequação da maioria dos nutrientes diminuiu com o aumento da renda e escolaridade; porém, o consumo excessivo de gordura saturada e o baixo consumo de fibra aumentaram com ambas variáveis. Grande parte dos nutrientes foi independentemente associada à renda e escolaridade, contudo, o consumo de ferro, vitamina B12 e sódio entre mulheres foi associado somente com a educação. Observou-se interação entre renda e escolaridade na associação com o consumo de sódio em homens, fósforo em mulheres e cálcio em ambos os sexos. Os achados indicam que melhorar a educação é um passo importante na melhoria do consumo de nutrientes no Brasil, além da necessidade de formulação de estratégias econômicas que permitam que indivíduos de baixa renda adotem uma dieta saudável. Nossos resultados mostram também um grande desafio das ações de saúde pública na área de nutrição, com importantes inadequações de consumo em toda população adulta brasileira e particularmente em grupos populacionais e regiões mais vulneráveis do país.
Resumo:
An advanced beam propagation model was developed to show that the far field narrows with good suppression of higher order modes for an appropriate temperature rise, without significant power penalty. To verify the accuracy of the model, the dependence of far field pattern on bias conditions were assessed both experimentally and theoretically, initially under pulsed conditions to reduce thermal effects. The results highlight the optimum taper angle and the role of local heating effects.
Resumo:
Lab-on-a-chip (LOC) is one of the most important microsystem applications with promise for use in microanalysis, drug development, diagnosis of illness and diseases etc. LOC typically consists of two main components: microfluidics and sensors. Integration of microfluidics and sensors on a single chip can greatly enhance the efficiency of biochemical reactions and the sensitivity of detection, increase the reaction/detection speed, and reduce the potential cross-contamination, fabrication time and cost etc. However, the mechanisms generally used for microfluidics and sensors are different, making the integration of the two main components complicated and increases the cost of the systems. A lab-on-a-chip system based on a single surface acoustic wave (SAW) actuation mechanism is proposed. SAW devices were fabricated on nanocrystalline ZnO thin films deposited on Si substrates using sputtering. Coupling of acoustic waves into a liquid induces acoustic streaming and motion of droplets. A streaming velocity up to ∼ 5cm/s and droplet pumping speeds of ∼lcm/s were obtained. It was also found that a higher order mode wave, the Sezawa wave is more effective in streaming and transportation of microdroplets. The ZnO SAW sensor has been used for prostate antigen/antibody biorecognition systems, demonstrated the feasibility of using a single actuation mechanism for lab-on-a-chip applications. © 2010 Materials Research Society.
Resumo:
Many types of oceanic physical phenomena have a wide range in both space and time. In general, simplified models, such as shallow water model, are used to describe these oceanic motions. The shallow water equations are widely applied in various oceanic and atmospheric extents. By using the two-layer shallow water equations, the stratification effects can be considered too. In this research, the sixth-order combined compact method is investigated and numerically implemented as a high-order method to solve the two-layer shallow water equations. The second-order centered, fourth-order compact and sixth-order super compact finite difference methods are also used to spatial differencing of the equations. The first part of the present work is devoted to accuracy assessment of the sixth-order super compact finite difference method (SCFDM) and the sixth-order combined compact finite difference method (CCFDM) for spatial differencing of the linearized two-layer shallow water equations on the Arakawa's A-E and Randall's Z numerical grids. Two general discrete dispersion relations on different numerical grids, for inertia-gravity and Rossby waves, are derived. These general relations can be used for evaluation of the performance of any desired numerical scheme. For both inertia-gravity and Rossby waves, minimum error generally occurs on Z grid using either the sixth-order SCFDM or CCFDM methods. For the Randall's Z grid, the sixth-order CCFDM exhibits a substantial improvement , for the frequency of the barotropic and baroclinic modes of the linear inertia-gravity waves of the two layer shallow water model, over the sixth-order SCFDM. For the Rossby waves, the sixth-order SCFDM shows improvement, for the barotropic and baroclinic modes, over the sixth-order CCFDM method except on Arakawa's C grid. In the second part of the present work, the sixth-order CCFDM method is used to solve the one-layer and two-layer shallow water equations in their nonlinear form. In one-layer model with periodic boundaries, the performance of the methods for mass conservation is compared. The results show high accuracy of the sixth-order CCFDM method to simulate a complex flow field. Furthermore, to evaluate the performance of the method in a non-periodic domain the sixth-order CCFDM is applied to spatial differencing of vorticity-divergence-mass representation of one-layer shallow water equations to solve a wind-driven current problem with no-slip boundary conditions. The results show good agreement with published works. Finally, the performance of different schemes for spatial differencing of two-layer shallow water equations on Z grid with periodic boundaries is investigated. Results illustrate the high accuracy of combined compact method.
Resumo:
Computational Fluid Dynamics CFD can be used as a powerful tool supporting engineers throughout the steps of the design. The combination of CFD with response surface methodology can play an important role in such cases. During the conceptual engineering design phase, a quick response is always a matter of urgency. During this phase even a sketch of the geometrical model is rare. Therefore, the utilisation of typical response surface developed for congested and confined environment rather than CFD can be an important tool to help the decision making process, when the geometrical model is not available, provided that similarities can be considered when taking into account the characteristic of the geometry in which the response surface was developed. The present work investigates how three different types of response surfaces behave when predicting overpressure in accidental scenarios based on CFD input. First order, partial second order and complete second order polynomial expressions are investigated. The predicted results are compared with CFD findings for a classical offshore experiment conducted by British Gas on behalf of Mobil and good agreement is observed for higher order response surfaces. The higher order response surface calculations are also compared with CFD calculations for a typical offshore module and good agreement is also observed. © 2011 Elsevier Ltd.