首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 377 毫秒
1.
A two-dimensional numerical mesoscale model is used to determine the pressure drag of sinoidal mountains and valleys in a neutral atmosphere. In the first part, pressure distributions and flow patterns for isolated obstacles are considered. For large aspect ratios, the pressure drag exerted by valleys becomes small compared to that of mountains. In the second part, interactions between several obstacles are investigated. For mountains, the drag on downstream obstacles is reduced considerably by the first obstacle when the obstacles are close together. For valleys there is a slight increase of the average drag exerted by each obstacle. In the limit for a large number of obstacles, average drag exerted by one mountain is equal to average drag for one valley. For smaller aspect ratios, this average drag can be entered into the resistence law from the Rossby number similarity theory to yield an effective roughness length.  相似文献   

2.
A previously published technique for using tethered spherical balloons as anemometers for measuring light low-level winds has been further developed. Earlier data on the relationship between the aerodynamic drag coefficient and the Reynolds number of spherical rubber balloons were combined with a large number of new data and re-analysed; and the errors in the relationship were estimated. The results allowed a more accurate calculation of wind speed from the deflection of a tethered balloon from the vertical. When combined with a new technique for calculating the effects of the tether, this enabled light to moderate low-level winds at fixed heights up to 600 m or more to be measured with simple, cheap, and readily mobile equipment; and a slight modification of the technique allowed measurement of winds in and above fog. Wind speeds measured by the ballon technique showed reasonably good agreement with measurements by an anemometer carried beneath the balloon.Glossary of Symbols a, b, c Coefficients in the relationship between lnC d and lnR - A Quantity under square root in solution for lnV whena0 - C d Wind drag coefficient for balloon - C dc Value ofC d given by calibration curve of Table I - D Dynamic wind pressure force on balloon - F Buoyant free lift of balloon with load - Re Reynold's number of balloon (sphere) - R = Re/105 - r Radius of sphere - T Tension in tether - V Wind speed - 83() =(lnC dc -lnC d ) when 83° , or 0 for other - Error in lnC d - Elevation of tether where attached to balloon - Elevation of balloon from ground tether point - Molecular viscosity of air - Ratio of circumference to diameter of circle - Density of air  相似文献   

3.
Summary Idealized numerical simulations using the Weather and Research Forecast (WRF) model indicate that three flow regimes, based on the moist Froude number, can be identified for a conditionally unstable, rotational, horizontally homogeneous, uniformly stratified flow over an idealized, three-dimensional, mesoscale mountain stretched spanwise to the impinging flow: (I) a quasi-stationary upslope convective system and an upstream-propagating convective system, (II) a quasi-stationary upslope convective system, and (III) a stationary upslope convective system and a quasi-stationary downstream convective system. Several major differences from a similar type of flow with no rotation over a two-dimensional mountain range are found. One important finding is that relatively strong mean flow produces a quasi-stationary mesoscale convective system (MCS) and maximum rainfall on the windward slope (upslope rain), instead of on the mountain peak or over the lee side.We found that the Coriolis force helps produce heavy upslope rainfall by making transition from flow-around the eastern part of the upslope to flow-over the western part of the upslope (transits to a higher flow regime) by deflecting the incident southerly flow to become east–southeasterly barrier winds. We found that the addition of the western flank of the arc-shaped mountain helps slow down the barrier wind from east and causes the maximum rainfall to move east of the windward slope. A lower-Froude number flow tends to produce a rainfall maximum near the concave region.Several other important facts can also be found in this study. The ratio of the maximum grid scale rainfall to the sub-grid scale rainfall increases when the moist Froude number increases. When the CAPE decreases, it is found that the upstream moist flow tends to shift to a higher Froude-number regime. Therefore, the Froude number cannot solely be used to define a moist flow regime when different CAPEs are considered. In another word, other parameters, such as CAPE, might play an important role in determining moist flow regimes.  相似文献   

4.
Surface pressure distribution and pressure drag on mountains   总被引：1，自引：0，他引：1  
Summary A mesoscale numerical model is used to compute the different components of the pressure drag on mountains, i.e.: form drag, wave drag, hydrostatic drag, and total pressure drag, for stable stratification. The Froude number is chosen so that non-breaking lee-waves evolve.The paper explains how the different parts of the drag are computed from the numerical results and how they form part of the horizontal momentum budget.For a single mountain the drag from the evolving stationary solution is compared to the wave drag from linear inviscid theory. Wave drag turns out to be about one third of the value expected from linear theory, and nonlinear interaction between wave and form drag is found. Wave drag is responsible for about 75% of the total drag if blocking is negligible.For two obstacles with varying distance the wave drag in the stationary solution varies between 5% and 30% of the value from linear theory due to partial cancellation between the lee-waves from the two mountains.Finally in an instationary simulation the passage of a cold air mass over a mountain and the respective drag components have been computed. 1500 m above the crestline of the obstacle wave drag is only 10% to 30% of the total drag.From the present results it seems realistic that wave drag from ALPEX experimental data was only a few percent of the value expected from the surface pressure distribution and from linear theory.With 7 Figures  相似文献   

5.
Functional forms of the universal similarity functions A, B (for wind components parallel and normal to the surface stress), and C (for potential temperature difference) are determined based on the generalized theory of the resistance laws for the Planetary Boundary Layer (PBL). The similarity-profile functions for the surface layer are matched with the velocity and temperature-defect profiles that are assumed to have shapes modified by certain powers of nondimensional height z/h, where h is the PBL height. The powers of the outer-layer profile functions are determined, so that the functions become negligible in the surface layer. To close the temperature defect law, an assumption that the temperature gradient across the top of the PBL is continuous with the stratification of the overlying atmosphere is used. The result of this assumption is that nondimensional momentum and temperature profiles in the PBL can be described in terms of four basic ratios: (1) roughness ratio = /h (2) scale-height ratio =|f|h/u*, (3) ambient stratification parameter =h/*, and (4) stability parameter =h/L, where L is the Monin-Obukhov length, z0 is the surface roughness, is the upper-air stratification, u * is the friction velocity, and * is the temperature scale at the surface. For stable conditions, the scale-height ratio can be related to the atmospheric stability and the upperair stratification, and the generalized similarity and Rossby number similarity theories become identical. Under appropriate boundary conditions, function A is explicitly dependent on the stability parameter , while B is a function of scale-height ratio , which in turn depends on the stability. Function C is shown to be dependent on the stability and the upper-air stratification, due to the closure assumption used for the temperature profile.The suggested functional forms are compared with other empirical approximations by several authors. The general framework used to determine the functional forms needs to be tested against good boundary-layer measurements.  相似文献   

6.
The extension of Lagrangian similarity theory of diffusion to stratified flow is examined, to improve its prediction of the vertical spread of a passive substance. In the basic equation, where is the average height of a cluster of particles,u * is the friction velocity andL is Monin-Obukhov length. It is shown theoretically, under the assumption of an equivalence between the diffusivities of heat and matter, that the unspecified function is the reciprocal of a more familiar meteorological parameter n , the dimensionless temperature gradient. The universal constantb is found to be approximately equal to von Karman's constant for various stability conditions. The predicted effect of stability on vertical spread shows excellent agreement with that of the published data from the O'Neill experiments.  相似文献   

7.
The structure of atmospheric turbulence in the surface layer over the open ocean is examined under conditions of local free convection. The raw data consist of profile and fluctuation measurements of wind and temperature as obtained from a meteorological buoy. For near neutral conditions and for waves running approximately along the wind direction, wave-induced wind fluctuations can be described by a simplified linear theory based on Miles (1957). In this case, the spectrum of wind velocity is given as the sum of two parts; for the turbulent part, the parameterization as obtained by Kaimal et al. (1972) applies, while the wave-induced part is parameterized using a simplification of Miles' linear theory. For cases of local free convection, the measurements of the vertical component of the wind velocity are well described by similarity theory; as expected, w /(-uw)1/2 is proportional to (- z/L)1/3. In order to scale the longitudinal wind velocity component, it seems to be reasonable to extend the list of relevant parameters by the height of the mixed layer z i. We obtain u /(- uw)1/2 (z/z i)1/3(- z/L)1/3 with only a poor correlation coefficient of r = 0.6. Overall, the results of local free convection scaling obtained from direct measurements show good agreement with those obtained from profile measurements. A comparison between direct and indirect determination of turbulent fluxes of momentum shows an unexplained difference of about 20%. This discrepancy is mainly due to a gap in the uw-cospectrum at the swell frequency.  相似文献   

8.
A Comparative Analysis of Transpiration and Bare Soil Evaporation   总被引：4，自引：0，他引：4  
Transpiration Ev and bare soil evaporation Eb processes are comparatively analysed assuming homogeneous and inhomogeneous areal distributions of volumetric soil moisture content . For a homogeneous areal distribution of we use a deterministic model, while for inhomogeneous distributions a statistical-deterministic diagnostic surface energy balance model is applied. The areal variations of are simulated by Monte-Carlo runs assuming normal distributions of .The numerical experiments are performed for loam. In the experiments we used different parameterizations for vegetation and bare soil surface resistances and strong atmospheric forcing. According to the results theEv()-Eb() differences are great, especially in dry conditions. In spite of this, the available energy flux curves of vegetation Av() and bare soil Ab() surfaces differ much less than the Ev() and Eb() curves. The results suggest that Ev is much more non-linearly related to environmental conditions than Eb. Both Ev and Eb depend on the distribution of , the wetness regime and the parameterization used. With the parameterizations, Eb showed greater variations than Ev. These results are valid when there are no advective effects or mesoscale circulation patterns and the stratification is unstable.  相似文献   

9.
The effect of changes in zonal and meridional atmospheric moisture transports on Atlantic overturning is investigated. Zonal transports are considered in terms of net moisture export from the Atlantic sector. Meridional transports are related to the vigour of the global hydrological cycle. The equilibrium thermohaline circulation (THC) simulated with an efficient climate model is strongly dependent on two key parameters that control these transports: an anomaly in the specified Atlantic–Pacific moisture flux (Fa) and atmospheric moisture diffusivity (Kq). In a large ensemble of spinup experiments, the values of Fa and Kq are varied by small increments across wide ranges, to identify sharp transitions of equilibrium THC strength in a 2-parameter space (between Conveyor On and Off states). Final states from this ensemble of simulations are then used as the initial states for further such ensembles. Large differences in THC strength between ensembles, for identical combinations of Fa and Kq, reveal the co-existence of two stable THC states (Conveyor On and Off)—i.e. a bistable regime. In further sensitivity experiments, the model is forced with small, temporary freshwater perturbations to the mid-latitude North Atlantic, to establish the minimum perturbation necessary for irreversible THC collapse in this bistable regime. A threshold is identified in terms of the forcing duration required. The model THC, in a Conveyor On state, irreversibly collapses to a Conveyor Off state under additional freshwater forcing of just 0.1 Sv applied for around 100 years. The irreversible collapse is primarily due to a positive feedback associated with suppressed convection and reduced surface heat loss in the sinking region. Increased atmosphere-to-ocean freshwater flux, under a collapsed Conveyor, plays a secondary role.  相似文献   

10.
Using a previous treatment of drag and drag partition on rough surfaces, simple analytic expressions are derived for the roughness length (z 0) and zero-plane displacement (d) of vegetated surfaces, as functions of canopy height (h) and area index (). The resulting expressions provide a good fit to numerous field and wind tunnel data, and are suitable for applications such as surface parameterisations in atmospheric models.  相似文献   

11.
Numerical simulation of turbulent convective flow over wavy terrain   总被引：1，自引：1，他引：0  
By means of a large-eddy simulation, the convective boundary layer is investigated for flows over wavy terrain. The lower surface varies sinusoidally in the downstream direction while remaining constant in the other. Several cases are considered with amplitude up to 0.15H and wavelength ofH to 8H, whereH is the mean fluid-layer height. At the lower surface, the vertical heat flux is prescribed to be constant and the momentum flux is determined locally from the Monin-Obukhov relationship with a roughness lengthz o=10–4 H. The mean wind is varied between zero and 5w *, wherew * is the convective velocity scale. After rather long times, the flow structure shows horizontal scales up to 4H, with a pattern similar to that over flat surfaces at corresponding shear friction. Weak mean wind destroys regular spatial structures induced by the surface undulation at zero mean wind. The surface heating suppresses mean-flow recirculation-regions even for steep surface waves. Short surface waves cause strong drag due to hydrostatic and dynamic pressure forces in addition to frictional drag. The pressure drag increases slowly with the mean velocity, and strongly with /H. The turbulence variances increase mainly in the lower half of the mixed layer forU/w *>2.  相似文献   

12.
Equilibrium evaporation beneath a growing convective boundary layer   总被引：1，自引：1，他引：0  
Expressions for the equilibrium surface Bowen ratio ( s ) and equilibrium evaporation are derived for a growing convective boundary layer (CBL) in terms of the Bowen ratio at the top of the mixed layer i and the entrainment parameter A R . If AR is put equal to zero, the solution for s becomes-that previously obtained for the zero entrainment or closed box model. The Priestley-Taylor parameter is also calculated and plotted in terms ofA R and i . Realistic combinations of the atmospheric parameters give values of in the range 1.1 to 1.4.  相似文献   

13.
Zusammenfassung In der vorliegenden Arbeit werden neue Klassifikationsprinzipien für Großwetterlagen entwickelt. Bisher wurde bei Wetterlagenklassifikationen das Druckfeld zugrunde gelegt, wobei quasistationären Druckzentren eine nicht berechtigte Vorrangstellung eingeräumt wurde. In der hier versuchten Klassifikation wird vom Strömungsfeld ausgegangen, das in elementare Formen zerlegt wird. Eine zu diesem Zweck durchgeführte statistische Untersuchung ergab, daß alle im Strömungsfeld auftretenden Zirkulationstypen auf drei Grundformen zurückgeführt werden können. Diese Grundformen sind:Driften, Wellen undWirbel.Die Untersuchung ergab im einzelnen, daß in mittleren Breiten der nördlichen Hemisphäre bei 49% aller untersuchten Fälle Driften, bei 23% Wellen und bei 28% Wirbel auftraten.In der hier durchgeführten Klassifikation wird das Druckfeld durch das Strömungsfeld und der Begriff Großwetterlage durch den umfassenderen Begriff des Zirkulationstyps ersetzt. Damit wird der unberechtigte Vorrang der Druckformen bei der Wetterlagenklassifikation aufgegeben. Die Klassifizierung der Zirkulationstypen ergibt sich schließlich durch Kombination der drei Zirkulationselemente: Drift, Welle und Wirbel.
Summary New principles of classification for large-scale weather situations are outlined in this paper. Hitherto the pressure-field has been taken as a basis for such classifications by conceding a precedence of an unjustified position to semi-permanent centres of pressure. The new classification starts from the field of large-scale motions, which is dissected in elementary models. A statistical test yielded the possibility to reduce all types of atmospheric circulations in the following three elementary models:drifts, waves andeddies.In detail it was found out, that drifts occur in 49%, waves in 23% and eddies in 28% of all cases investigated.In the new classification the term pressure-field is substituted by field of motion and the expression large-scale weather situation by the more comprehensive conception type of circulation. By that the unjustified priority of pressure-centers in classifying weather situations is abolished. At last the classification of the types of circulation follows from a combination of the three elementary models: drift, wave and eddy.
résumé La présente étude développe de nouveaux principes de classification des situations météorologiques. Alors que jusqu'ici on s'est fondé sur le champ de pression ce qui conduisait à attribuer aux centres d'action quasi stationnaires un rôle trop important, l'auteur part ici du champ de mouvement décomposé en formes élémentaires. Un examen statistique lui a montré que tous les types de circulation peuvent se ramener à trois formes fondamentales:courants, ondulations ettourbillons.Aux latitudes moyennes de l'hémisphère Nord les courants représentent le 49%, les ondulations le 23% et les tourbillons le 28%.Au champ de pression se substitue donc le champ de courant, et les situations météorologiques se groupent en types de circulation ce qui supprime le rôle prépondérant des formes isobariques. Le classement final des types de circulation résulte de la combinaison des trois types mentionnés:courants, ondulations ettourbillons.

Mit 15 Textabbildungen  相似文献   

14.
Summary The discrepancy between measured (or estimated) visual range (V g) and the one calculated from the drop size distribution (V b), stated byTrappenberg on the basis of the measurements ofRittberger, is discussed. It is assumed that in the measurement of the drop size distribution the small droplets are discriminated. If the distributions are extrapolated to drop diameters of 2 or 1 ,V b turns out to be in agreement withV g.
Diskussionsbeitrag zur Frage der Sicht in Wolken
Zusammenfassung Die vonTrappenberg auf Grund der Messungen vonRittberger festgestellte Diskrepanz zwischen gemessener (oder geschätzter) SichtweiteV g und der aus der Tropfengrößenverteilung berechnetenV b wird diskutiert. Es wird angenommen, daß bei der Messung der Tropfengrößenverteilung die kleinen Tropfen nicht genügend zur Geltung kommen. Extrapoliert man die Verteilungen in zwangloser Weise bis zu Tropfendurchmessern von 2 oder 1 , so ergibt sich eine gute Übereinstimmung vonV b mitV g.
Résumé L'auteur discute la discrépance existant entre la visibilité mesurée (ou estimée)Vg et la visibilité calculéeV b au moyen de la répartition de la grosseur des gouttes, discrépance établie parTrappenberg en partant des mesures deRittberger. On admet alors que, lors de la mesure de la répartition de la grosseur des gouttes, les plus petites d'entre elles n'ont pas suffisament de poids. Si l'on extrapole cette répartition jusqu'à un diamètre des gouttes de 2 ou 1 , on constate alors une bonne concordance entreV b etV g.

Mit 1 Textabbildung  相似文献   

15.
This paper considers the near-field dispersion of an ensemble of tracer particles released instantaneously from an elevated source into an adiabatic surface layer. By modelling the Lagrangian vertical velocity as a Markov process which obeys the Langevin equation, we show analytically that the mean vertical drift velocity w(t) is w()=bu *(1–e –(1+)), where is time since release (nondimensionalized with the Lagrangian time scale at the source), b Batchelor's constant, and u *, the friction velocity. Hence, the mean height and mean depth of the ensemble are calculated. Although the derivation is formally valid only when 1, the predictions for w, mean height and mean depth are consistent in the downstream limit ( 1) with surface-layer Lagrangian similarity theory and with the diffusion equation. By comparing the analytical predictions with numerical, randomflight solutions of the Langevin equation, the analytical predictions are shown to be good approximations at all times, both near-field and far-field.  相似文献   

16.
In this paper we analyse diabatic wind profiles observed at the 213 m meteorological tower at Cabauw, the Netherlands. It is shown that the wind speed profiles agree with the well-known similarity functions of the atmospheric surface layer, when we substitute an effective roughness length. For very unstable conditions, the agreement is good up to at least 200 m or z/L–7(z is height, L is Obukhov length scale). For stable conditions, the agreement is good up to z/L1. For stronger stability, a semi-empirical extension is given of the log-linear profile, which gives acceptable estimates up to ~ 100 m. A scheme is used for the derivation of the Obukhov length scale from single wind speed, total cloud cover and air temperature. With the latter scheme and the similarity functions, wind speed profiles can be estimated from near-surface weather data only. The results for wind speed depend on height and stability. Up to 80 m, the rms difference with observations is on average 1.1 m s–1. At 200 m, 0.8 m s–1 for very unstable conditions increasing to 2.1 m s–1 for very stable conditions. The proposed methods simulate the diurnal variation of the 80 m wind speed very well. Also the simulated frequency distribution of the 80 m wind speed agrees well with the observed one. It is concluded that the proposed methods are applicable up to at least 100 m in generally level terrain.  相似文献   

17.
Summary ¶This study presents a numerical simulation of the bora wind as it occurs in form of a severe wind blowing down coastal mountains and over the Adriatic Sea. A typical cyclonic bora event, occurring during the period January 3–6, 1995, is simulated using a nested limited area model. An integration, with horizontal resolution of about 14km, and a nested one, with higher resolution, about 5.5km, are presented. The 1997 version of the Eta Model is used for both the lower resolution and the higher resolution runs. Numerous details of the simulation are found to be in good agreement with the understanding as well as the observational knowledge of the bora, thus supporting confidence in the realism of the results. In particular, features of the simulated flow are seen strongly indicative of some basic characteristics of the hydraulic model of the phenomenon, such as the mountain wave breaking and the upstream flow acceleration. Moreover, the increase in horizontal resolution, in combination with an improvement of the coastal SST information, led to a still improved realism of the low-level wind representation over the Adriatic Sea.Received December 31, 2001; revised March 25, 2002; accepted July 19, 2002 Published online: February 20, 2003  相似文献   

18.
Summary The investigations during the German Frontal Experiment 1987 show the modification of structure and movement of cold fronts in the vicinity of the Alps and in an Alpine valley. The two cold fronts that moved nearly parallel to the mountains were modified at the mountain barrier by the channelling effect. In those cases the postfrontal air flowed into the Inn valley from the valley's exit up-valley. The upper and lower parts of the fronts became detached, and the structure of the lower part changes on its way from the foreland to the mountain barrier.In contrast the fourth cold front, which moved from north-west to south-east, crossed the mountains without appreciable deformation although it suffered some degree of blocking and retardation.With 7 Figures  相似文献   

19.
TheConvectiveDiffusionObserved byRemoteSensors (CONDORS) field experiment conducted at the Boulder Atmospheric Observatory used innovative techniques to obtain three-dimensional mappings of plume concentration fields, /Q, of oil fog detected by lidar and chaff detected by Doppler radar. It included extensive meteorological measurements and, in 1983, tracer gases measured at a single sampling arc. Final results from ten hours of elevated and surface release data are summarized here. Many intercomparisons were made. Oil fog /Q measured 40m above the arc are mostly in good agreement withSF 6 values, except in a few instances with large spacial inhomogeneities over short distances. After a correction scheme was applied to compensate for the effect of its settling speed, chaff dy/Q agreed well with those of oil except in two cases of oil fog hot spots. Mass or frequency distribution vs. azimuth or elevation angle comparisons were made for chaff, oil, and wind, with mostly good agreements. Spacial standard deviations, y and z, of chaff and oil agree overall and are consistent at short range with velocity standard deviations vand w 0.6w* (the convective scale velocity), as measured atz>100m. Surface release y is enhanced up to 60% at smallx, consistent with the Prairie Grass measurements and with larger v and reduced wind speed measured near the surface. Decreased y at small dimensionless average times is also noted. Finally, convectively scaled dy, C y, were plotted versus dimensionlessx andz for oil, chaff, and corrected chaff for each 30–60 min period. Aggregated CONDORSC y fields compare well with laboratory tank and LES numerical simulations; surface-released oil fog compares expecially well with the tank experiments. However, large deviations from the norm occurred in individual averaging periods; these deviations correlated strongly with anomalies in measured distributions.On assignment to the US Environmental Protection Agency, Atmospheric Research and Exposure Assessment Laboratory, RTP, NC.  相似文献   

20.
We have devised a partial differential equation for the prediction of dust concentration in a thin layer near the ground. In this equation, erosion (detachment), transport, deposition and source are parameterised in terms of known quantities. The interaction between a wind prediction model in the boundary layer and this equation affects the evolution of the dust concentration at the top of the surface layer. Numerical integrations are carried out for various values of source strength, ambient wind and particle size. Comparison with available data shows that the results appear very reasonable and that the model should be subjected to further development and testing.Notation (x, y, z, t) space co-ordinates and time (cm,t) - u, v components of horizontal wind speed (cm s–1) - u g, vg components of the geostrophic wind (cm s–1) - V=(u2+v2)1/2 (cm s–1) - (û v)= 1/(h – k) k h(u, v)dz(cm s–1) - V * friction velocity (cm s–1) - z 0 roughness length (cm) - k 1 von Karman constant =0.4 - V d deposition velocity (cm s–1) - V g gravitational settling velocity (cm s–1) - h height of inversion (cm) - k height of surface layer (cm) - potential temperature (°K) - gr potential temperature at ground (°K) - K potential temperature at top of surface layer (°K) - P pressure (mb) - P 0 sfc pressure (mb) - C p/Cv - (t)= /z lapse rate of potential temperature (°K cm–1) - A(z) variation of wind with height in transition layer - B(z) variation of wind with height in transition layer - Cd drag coefficient - C HO transfer coefficient for sensible heat - C dust concentration (g m–3) - C K dust concentration at top of surface layer (g m–3) - D(z) variation with height of dust concentration - u, v, w turbulent fluctuations of the three velocity components (cm s–1) - A 1 constant coefficient of proportionality for heat flux =0.2 - Ri Richardson number - g gravitational acceleration =980 cm s–2 - Re Reynolds number = - D s thickness of laminar sub-layer (cm) - v molecular kinematic viscosity of air - coefficient of proportionality in source term - dummy variable - t time step (sec) - n time index in numerical equations On sabbatical leave at University of Aberdeen, Department of Engineering, September 1989–February 1990.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号