一个改进的随机动力学水平衡模型及其应用研究Ⅰ.模型

本文在Eagleson随机动力学水平衡模型的基础上 ,改进了原模型中蒸腾、蒸发量等项的计算方法 ;并根据黄土区降水分布特点 ,将全年分为 1 6个时间长度不等的降水时段 ,分别统计降水特征参数 ,利用质量守衡定律和逐步校正演算法 ,可计算出水文生态系统平衡体系水均衡要素的定量分配和活动层土壤含水量的年动态变化 ,代表了原模型中确定土壤含水量的简单统计法 ,从而可提高模型计算精度和实用性 .改进后的随机动力学水平衡模型更适用地黄土高原的特定水文生态系统

SOES3010/6005 - Scaling analysis (Week 1, end of lecture)

This is the end of the scaling analysis we saw in class on Friday. In class, we managed to scale the mass conservation equation and the x-momentum equation, but we didn't finish scaling the z-momentum equation in order to arrive at the hydrostatic approximation.

The influence of hilly terrain on canopy-atmosphere carbon dioxide exchange

Topography influences many aspects of forest-atmosphere carbon exchange; yet only a small number of studies have considered the role of topography on the structure of turbulence within and above vegetation and its effect on canopy photosynthesis and the measurement of net ecosystem exchange of CO2 (N-ee) using flux towers. Here, we focus on the interplay between radiative transfer, flow dynamics for neutral stratification, and ecophysiological controls on CO2 sources and sinks within a canopy on a gentle cosine hill. We examine how topography alters the forest-atmosphere CO2 exchange rate when compared to uniform flat terrain using a newly developed first-order closure model that explicitly accounts for the flow dynamics, radiative transfer, and nonlinear eco physiological processes within a plant canopy. We show that variation in radiation and airflow due to topography causes only a minor departure in horizontally averaged and vertically integrated photosynthesis from their flat terrain values. However, topography perturbs the airflow and concentration fields in and above plant canopies, leading to significant horizontal and vertical advection of CO2. Advection terms in the conservation equation may be neglected in flow over homogeneous, flat terrain, and then N-ee = F-c, the vertical turbulent flux of CO2. Model results suggest that vertical and horizontal advection terms are generally of opposite sign and of the same order as the biological sources and sinks. We show that, close to the hilltop, F-c departs by a factor of three compared to its flat terrain counterpart and that the horizontally averaged F-c-at canopy top differs by more than 20% compared to the flat-terrain case.

Salinity changes in the world ocean since 1950 in relation to changing surface freshwater fluxes

Global hydrographic and air–sea freshwater flux datasets are used to investigate ocean salinity changes over 1950–2010 in relation to surface freshwater flux. On multi-decadal timescales, surface salinity increases (decreases) in evaporation (precipitation) dominated regions, the Atlantic–Pacific salinity contrast increases, and the upper thermocline salinity maximum increases while the salinity minimum of intermediate waters decreases. Potential trends in E–P are examined for 1950–2010 (using two reanalyses) and 1979–2010 (using four reanalyses and two blended products). Large differences in the 1950–2010 E–P trend patterns are evident in several regions, particularly the North Atlantic. For 1979–2010 some coherency in the spatial change patterns is evident but there is still a large spread in trend magnitude and sign between the six E–P products. However, a robust pattern of increased E–P in the southern hemisphere subtropical gyres is seen in all products. There is also some evidence in the tropical Pacific for a link between the spatial change patterns of salinity and E–P associated with ENSO. The water cycle amplification rate over specific regions is subsequently inferred from the observed 3-D salinity change field using a salt conservation equation in variable isopycnal volumes, implicitly accounting for the migration of isopycnal surfaces. Inferred global changes of E–P over 1950–2010 amount to an increase of 1 ± 0.6 % in net evaporation across the subtropics and an increase of 4.2 ± 2 % in net precipitation across subpolar latitudes. Amplification rates are approximately doubled over 1979–2010, consistent with accelerated broad-scale warming but also coincident with much improved salinity sampling over the latter period.

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.

Transferência de calor combinando radiação e convecção no interior de dutos de geradores de vapor fumotubulares

Esta dissertação de mestrado considera a transferência de calor combinando convecção e radiação térmica no escoamento de gases participantes em dutos de seção circular. Partindo de uma metodologia geral, o trabalho enfoca principalmente os casos típicos de aplicação em geradores de vapor fumotubulares de pequeno e médio porte, em que gases em alta temperatura escoam através de um tubo mantido em temperatura uniforme. O escoamento é turbulento e o perfil de velocidade é plenamente desenvolvido desde a entrada do duto. A temperatura do gás, contudo, é uniforme na entrada, considerando-se a região de desenvolvimento térmico. Duas misturas de gases são tratadas, ambas constituídas por dióxido de carbono, vapor d’água e nitrogênio, correspondendo a produtos típicos da combustão estequiométrica de óleo combustível e metano. As propriedades físicas dos gases são admitidas uniformes em todo o duto e calculadas na temperatura de mistura média, enquanto que as propriedades radiantes são modeladas pela soma-ponderada-de-gases-cinzas. O campo de temperatura do gás é obtido a partir da solução da equação bidimensional da conservação da energia, sendo os termos advectivos discretizados através do método de volumes de controle com a função de interpolação Flux-Spline; as trocas de energia radiantes são avaliadas por meio do método das zonas, onde cada zona de radiação corresponde a um volume de controle. Em um primeiro passo, a metodologia é verificada pela comparação com resultados apresentados na literatura para a transferência de calor envolvendo apenas convecção e combinando convecção com radiação. Em seguida, discutem-se alguns efeitos da inclusão da radiação térmica, por exemplo, no número de Nusselt convectivo e na temperatura de mistura do gás. Finalmente, são propostas correlações para o número de Nusselt total, que leva em conta tanto a radiação quanto a convecção. Essa etapa exige inicialmente uma análise dos grupos adimensionais que governam o processo radiante para redução do número elevado de parâmetros independentes. As correlações, aplicáveis a situações encontradas em geradores de vapor fumotubulares de pequeno e médio porte, são validadas estatisticamente pela comparação com os resultados obtidos pela solução numérica.

Modelos analíticos para rádio fontes extragalácticas

The great majority of analytical models for extragalactic radio sources suppose self-similarity and can be classified into three types: I, II and III. We have developed a model that represents a generalization of most models found in the literature and showed that these three types are particular cases. The model assumes that the area of the head of the jet varies with the jet size according to a power law and the jet luminosity is a function of time. As it is usually done, the basic hypothesis is that there is an equilibrium between the pressure exerted both by the head of the jet and the cocoon walls and the ram pressure of the ambient medium. The equilibrium equations and energy conservation equation allow us to express the size and width of the source and the pressure in the cocoon as a power law and find the respective exponents. All these assumptions can be used to calculate the evolution of the source size, width and radio luminosity. This can then be compared with the observed width-size relation for radio lobes and the power-size (P-D) diagram of both compact (GPS and CSS) and extended sources from the 3CR catalogue. In this work we introduce two important improvement as compared with a previous work: (1)We have put together a larger sample of both compact and extended radio sources

Análise do escoamento de HFC-134a em tubos capilares usando o modelo de dois fluidos

This work presents a numerical model to simulate refrigerant flow through capillary tubes, commonly used as expansion devices in refrigeration systems. The capillary tube is considered straight and horizontal. The flow is taken as one-dimensional and adiabatic. Steady state and thermodynamic equilibrium conditions are assumed. The two-fluid model, involving four conservation equations and considering the hidrodynamic nonequilibrium between the liquid and vapor phases is applied to the flow region. The pressure profiles and the mass flow rates given by the model are compared with experimental data.

Object oriented modeling of piston engines

This work reports a conception phase of a piston engine global model. The model objective is forecast the motor performance (power, torque and specific consumption as a function of rotation and environmental conditions). Global model or Zero-dimensional is based on flux balance through each engine component. The resulting differential equations represents a compressive unsteady flow, in which, all dimensional variables are areas or volumes. A review is presented first. The ordinary differential equation system is presented and a Runge-Kutta method is proposed to solve it numerically. The model includes the momentum conservation equation to link the gas dynamics with the engine moving parts rigid body mechanics. As an oriented to objects model the documentation follows the UML standard. A discussion about the class diagrams is presented, relating the classes with physical model related. The OOP approach allows evolution from simple models to most complex ones without total code rewrite. Copyright © 2001 Society of Automotive Engineers, Inc.

Análise numérica de condensadores do tipo arame-sobre-tubo usados em refrigeradores domésticos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Modelagem do escoamento ao longo de evaporadores de serpentina com tubos aletados

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Análise teórica e experimental de uma ventaneira de alto forno, visando seu comportamento térmico e fluidodinâmico

In this work, it was studied the thermal and fluid dynamic behavior from theorical and experimental point of view of a blast furnace tuyere. The tuyere is responsible for the hot air insuflation into the blast furnace. The parameter used in the comparison was the difference between the cooling water inlet and outlet temperatures. There were used forced convection correlations inside circular sections with adequations for non circular sections. Based on operations dates of flux and thermal loads it was possible to model numerically the tuyere, and, since it was obtained the wall temperatures, estimate the conduction and convection resistances and the heat flux through the walls in contact with the water. The total heat fluxes from wall to water were applied to the energy conservation equation where could be estimated the theorical temperature variation. Compared to the real value, the theorical value presented a difference of 0,2 °C. Considering that the boundary conditions around the tuyere are transitories and that your channels have some rugosity irregularities we can consider the estimation method for cooling system coherent with the real operational parameters

Modelo distribuído aplicado à análise de evaporadores do tipo tubo aletado

Pós-graduação em Engenharia Mecânica - FEIS

(Table 1) Horizontal and vertical ice flow velocities for seven markers around Vostok Station

Combined geodetic, geophysical and glaciological in situ measurements are interpreted regarding surface height changes over subglacial Lake Vostok and the local mass balance of the ice sheet at Vostok station. Repeated GPS observations spanning 5 years and long-term surface accumulation data show that the height of the lake surface has not changed over the observation period. The application of the mass conservation equation to purely observational data yields an ice mass balance for Vostok station close to equilibrium.

Three-dimensional kinetic simulations of active suspensions: effect of chemotaxis and steric interaction

In this work, the effects of chemotaxis and steric interactions in active suspensions are analyzed by extending the kinetic model proposed by Saintillan and Shelley [1, 2]. In this model, a conservation equation for the active particle configuration is coupled to the Stokes equation for the flow arising from the force dipole exerted by the particles on the fluid. The fluid flow equations are solved spectrally and the conservation equation is solved by second-order finite differencing in space and second-order Adams-Bashforth time marching. First, the dynamics in suspensions of oxytactic run-and-tumble bacteria confined in thin liquid films surrounded by air is investigated. These bacteria modify their tumbling behavior by making temporal comparisons of the oxygen concentration, and, on average, swim towards high concentrations of oxygen. The kinetic model proposed by Saintillan and Shelley [1, 2] is modified to include run-and-tumble effects and oxygentaxis. The spatio-temporal dynamics of the oxygen and bacterial concentration are analyzed. For small film thicknesses, there is a weak migration of bacteria to the boundaries, and the oxygen concentration is high inside the film as a result of diffusion; both bacterial and oxygen concentrations quickly reach steady states. Above a critical film thickness (approximately 200 micron), a transition to chaotic dynamics is observed and is characterized by turbulent-like 3D motion, the formation of bacterial plumes, enhanced oxygen mixing and transport into the film, and hydrodynamic velocities of magnitudes up to 7 times the single bacterial swimming speed. The simulations demonstrate that the combined effects of hydrodynamic interactions and oxygentaxis create collective three-dimensional instabilities which enhances oxygen availability for the bacteria. Our simulation results are consistent with the experimental findings of Sokolov et al. [3], who also observed a similar transition with increasing film thickness. Next, the dynamics in concentrated suspensions of active self-propelled particles in a 3D periodic domain are analyzed. We modify the kinetic model of Saintillan and Shelley [1, 2] by including an additional nematic alignment torque proportional to the local concentration in the equation for the rotational velocity of the particles, causing them to align locally with their neighbors (Doi and Edwards [4]). Large-scale three- dimensional simulations show that, in the presence of such a torque both pusher and puller suspensions are unstable to random fluctuations and are characterized by highly nematic structures. Detailed measures are defined to quantify the degree and direction of alignment, and the effects of steric interactions on pattern formation will be presented. Our analysis shows that steric interactions have a destabilizing effect in active suspensions.